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PREFACE 


The object of this book is to consider from an elementary 
standpoint as many different types of wave motion as 
possible. In almost every case the fundamental problem 
is the same, since it consists in solving the standard equation 
of wave motion ; the various applications differ chiefly 
in the conditions imposed upon these solutions. For tins 
reason it is desirable that the subject of waves should bo 
treated as one whole, rather than in several distinct parts ; 

the present tendency is in this direction. 

It is presupposed that the reader is familiar with the 
elements of vector analysis, the simpler results of which 
are freely quoted. In a sense this present volume may 
be regarded as a sequel to Rutherford’s Vector Methods, 
published in this series. 

In a volume of this size, it is not possible to deal 
thoroughly with any one branch of the subject : nor 
indeed is this desirable in a book which is intended as 
an introduction to the more specialised and elaborate 
treatises necessary to the specialist. This book is intended 
for University students covering a general course of Appli<"d 
Mathematics or Natural Philosophy in the final year of 
their honours degree. A few topics, such as elastic waves 
in continuous media, or at the common boundary of two 
media, and radiation from aerials, have unavoidably had 
to be omitted for lack of space. The reader is referred to 
any of the standard works on elasticity and wireless for 
a discussion of these problems. 
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This book would not be complete without a reference 
of gratitude to my friends Dr D. E. Rutherford and 
Prof. G. S. Rushbrooke, who have read the proofs, checked 
most of the examples and contributed in no small way 
to the clarity of my arguments. My thanks are also 
offered to my wife for her share in the preparation of 
the manuscript. 

C. A. C. 

January 1941. 
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OHAPTER I 


THE EQUATION OF WAVE MOTION 


§ 1. We are all familiar with the idea of a wave ; thus, 
when a pebble is dropped into a pond, water waves travel 
radially outwards ; when a piano is played, the wires 
vibrate and sound waves spread through the room ; when 
a wireless station is transmitting, electric waves move 
through the ether. These are all examples of wave motion, 
and they have two important properties in common : 
firstly, energy is propagated to distant points ; and 
secondly, the disturbance travels tlirough the medium 
without giving the medium as a whole any permanent 
displacement. Thus the ripples spread outwards. ver a 
pond carrying energy with them, but as we can see by 
watching the motion of a small floating body, the water 
of the pond itself does not move with the waves. In the 
following chapters we shall find that whatever the nature 
of the medium which transmits the waves, whether it be 
air, a stretched string, a liquid, an electric cable or the 
ether, these two properties which arc common to all these 
types of wave motion, will enable us to relate them 
together. They are all governed by a certain differentia 
equation, the Equation of Wave Motion (sec §5), and 
the mathematical part of each separate problem merely 
consists in solving this equation with the right boundary 
conditions, and then interpreting the solution appropriately. 

8 2. Consider a disturbance <f> which is propagated 

along the x axis with velocity c. There is no need to 

A 

1 
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state explicitly what <f> refers to ; it may be the elevation 
of a water wave or the magnitude of a fluctuating electric 
field. Then, since the disturbance is moving, <f) will depend 
on x and t. When l = 0, <f> will be some function of x 
which we may call f(x). f(x) is the wave profile, since 
if we plot the disturbance (f> against x, and “ photograph ” 
the wave at t = 0, the curve obtained will be <f> = f(x). 
If we suppose that the wave is propagated without change 
of shape, then a photograph taken at a later time l will 
bo identical with that at t == 0, except that the wave 
profile has moved a distance cl in the positive direction 
of the x axis. If we took a new origin at the point x = ct, 
and let distances measured from this origin be called X, 
so that x = A'-f ct, then the equation of the wave profile 
referred to this new origin would be 

<t>=f(X)- 

Referred to the original fixed origin, this means that 

.... ( 1 ) 

I his equation is the most general expression of a wave 
moving with constant velocity c and without change of 
shape, along the positivo direction of x. If the wave is 
travelling in the negative direction its form is given by 
(1) with the sign of c changed, i.e. 

4>=f(x+ct) (2) 

g 3. The simplest example of a wave of this kind is the 
harmonic wave, in which the wave profile is a sine or 
cosme curve. Thus if the wave profile at t = 0 is 

('£)<= o = a cos mx, 

then at time t, the displacement, or disturbance, is 

<f> = a cos m(x-d) ... (3) 

The maximum value of the disturbance, viz. a, is called 
the amplitude. I he wave profile repeats itself at regular 
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distances 277 /m. This is known as the wavelength A. 
Equation (3) could therefore be written 

2tt 

<f> = a cos —{x— cl) .... (4) 


The time taken for one complete wave to pass any point 
is called the period r of the wave. It follows from (4) that 


^(x—cl) must pass through a complete cycle of values 
A 


as t is increased by r. Thus 


277CT 

“T 



i.e. 




The frequency n of the wave is the number of waves 
passing a fixed observer in unit time. Clearly 

n = 1/r . . . . (0) 

so that c = riA (7) 


and equation (4) may be written in either of the equivalent 
forms, 




( 8 ) 

(9) 


Sometimes it is usoful to introduce the wave number k, 
which is the number of wavos in unit distance. Then 

k = 1/A (10) 


and we may write equation (9) 

<f> = a cos 2ir(kx—nl) 


. ( 11 ) 
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If we compare two similar waves 

<f>i~ a cos 2ir(k‘x—nt), 

<f> 2 = a cos{2tt(Icx— nl) -f-c}, 

we see that <f> 2 is the same as <f> l except that it is displaced 
a distance e!'2irlc, i.e. eXj'ln. € is called the phase of<f> 3 
relative to </>,. If € = 2 tt, 477, ... then the displacement 
is exactly one, two, ... wavelengths, and we say that the 
waves are in phase ; if e = it, 37 r, ... then the two waves 
are exactly out of phase. 

Even if a wave is not a harmonic wave, but the wave 
profile consists of a regularly repeating pattern, the 
definitions of wavelength, period, frequency and wave 
number still apply, and equations (5), (6), (7) and (10) 
are still valid. 

§ 4. It is possible to generalise equation (1) to deal 
with the case of plane waves in three dimensions. A 
plane wave is one in which the disturbance is constant 
over all points of a plane drawn perpendicular to the 
direction of propagation. Such a plane is called a wave- 
front, and the wavefront moves perpendicular to itself 
with the velocity of propagation c. If the direction of 
propagation is x : y : z = l : m : n, where /, m, n are the 
direction cosines of the normal to each wavefront, then the 
equation of the wavefronts is 

Ix+my + m = const., . . . (12) 

and at any moment /, </> is to be constant for all x, y, z 
satisfying (12). It is clear that 

<f> =f(l * +»»»/ + nz-ci) . . . (13) 

is a function which fulfils all these requirements and 
therefore represents a plane wave travelling with velocity 
c in the direction l :m:n without change of form. 
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§5. The expression (13) is a particular solution of the 
equation of wave motion referred to on p. 1. Since 
l, m, n are direction cosines, l 2 +m 2 -\-n 2 = 1, and it is 
easily verified that <f> satisfies the differential equation * 



dty c*4> 

cx 2 8y 2 ' cz 1 


i ^ 

c 2 ct- 



Tliis is the equation of wave motion. | It is one of the 
most important differential equations in the whole of 
mathematics, since it represents all types of wave motion 
in which the velocity is constant. The expressions in 
(1), (2), (£), (9), (11) and (13) are all particular solutions 
of this equation. We shall find, as we investigate different 
types of wave motion in subsequent chapters, that equation 
(14) invariably appears, and it will be our task to select 
the solution that is appropriate to our particular problem. 
There are certain types of solution that occur often, and 
we shall discuss some of them in the rest of this chapter, 
but before doing so, there is one important property of 
the fundamental equation that must be explained. 


§ 6. The equation of wave motion is linear. That is 
to say, <f> and its differential coefficients never occur in 
any form other than that of the first degree. Consequently, 
if <£, and <f> 2 are any two solutions of (14), ai<f> t +a 2 <f> 2 is 
also a solution, a, and a 2 being two arbitrary constants. 
This is an illustration of the principle o t superposition, 
which states that, when all the relevant equations are 
linear, we may superpose any number of individual 
solutions to form new functions which are themselves also 
solutions. We shall often have occasion to do this. 

A particular instance of this superposition, which is 

• This equation hus a dose resemblance to Laplace/ a Equation 
which ifl discussed in Rutherford, Vector Method *, Chapter VII. 
The roforoneee, hero and later, to this book refer to tho first edition. 

t Soinot lines called the wave equation, but wo do not use this 
phnmo to uvoid confusion with inodorn wave mechanics. 
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important in many problems, comes by adding together 
two harmonic waves going in different directions with the 
same amplitude and velocity. Thus, with two waves 
similar to (11) in opposite directions, we obtain 

<f> = a cos 2tt(Icx— nt)-\-a cos 2-n(kx-\-nl) 

= 2 a cos 2ttJcx cos 27ml (15) 


This is known as a stationary wave, to distinguish it from 
the earlier progressive waves. It owes its name to the 
fact that the wave profile does not move forward. In fact, 
<f) always vanishes at the points for which cos 2-nkx = 0, 
13 5 

viz. x = ±j7» ±77. ±77, ... • These points are called the 

nodes, and the intermediate points, where the amplitude 
of (f) (i.e. 2a cos 2ttI?x) is greatest, are called antinodes. 
The distance between successive nodes, or successive 
antinodes, is lj2lc, which, by (10), is half a wavelength. 

Using harmonic wave functions similar to (13), we find 
stationary waves in three dimensions, given by 

. 2t7 2tt 

9 = a cos y (lx \-my+nz-ct ) + a cos — (Ix+my+nz+ci) 


2rr 


2tt 


2a cos — (fr-fmy-f 7i2) cos ct 
A A 


(16) 


In thij case ^ always vanishes on the planes lx-\-my-\-nz 

A 3A 

— ± ± and these are known as nodal planes. 


§ 7. We shall now obtain some special types of solution 
of the equation of wave motion ; wo shall then be able to 
apply them to specific problems in later chapters. Wo 
may divide our solutions into two main types, representing 
stationary and progressive waves. 

W 'j have already dealt with progressive waves in one 
dm, "sion. The equation to be solved is 

_ \ ey 

dx 2 ~ c 2 Tt 2 • 
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Its most general solution may be obtained by a 

method due to D’Alembert. We change to new variables 

u = x—ct, and v = x+ct. Then it is easily verified that 

dJ) 8<h d<b dd> r 8(f> d(f> 

^ - transforms to ^ — transforms to -c — +c — ; 

n rt'ti slti sit S' // Si I) 


dx 

d*<f> 

dc 2 


du dv dt 

, o r-t . &<t> 

du 2 


dudv 8v 


becomes ^ +2 and finally becomes 

Cl “ 


When these chang e are made 

[du 2 dudv dv 2 ) 


p2i 

the equation becomes =0; the most general solution 

dudv 

of this is (f) = f(u) 4-<7(id. 

/ and g being arbitrary functions. In the original variables 
this is <f> = f(r— ct)-\-g{x+ct) . . . (17) 

The harmonic waves of § 2 are special cases of this, in 
which / and g are cosine functions. The waves / and g 
travel with velocity c, in opposite directions. 

In two dimensions the equation of wave motion is 


d 2 <f> d 2 <f> _ 1 d 2 <f> 

dx 2 + By* ~ c 2 dt 2 


(18) 


By a method similar to D’Alembert’s, it can be shown 
that the most general solution involving only plane * waves 
is <f> =f(lx-\-wj— cl)+g(lx-\-my+cl), . (19) 

where, as before, / and g arc arbitrary functions and 
l 2 +m 2 = 1. 

In three dimensions the differential equation is 


d 2 <f> d 2 <f> d 2 <f>_ l efrf) 

dx 2 + dy 2 + dz 2 ~ c 2 dt 2 


( 20 ) 


and the most general solution involving only plane waves is 
(f) = f(lx+my+nz—cl)+g(lj:+iny+nz+cl) . (21) 

in which l 2 -\-m 2 +n 2 = 1. 

• Strictly these should bo called lino waves, since at any 
moment ^ is constant along tho linos lx-\- my «= const. 
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There are, however, other solutions of progressive type, 
not involving plane waves. For suppose tliat we transform 
(20) to spherical polar coordinates r, 9, ip * The equation 
of wave motion becomes 


?V,2ty a *l>\ , i &<f> l g»0 

Sr 2 r dr r- sin 9 c9 \ c9j + r 2 sin* 9 difi 2 c 2 dt 2 ' 

( 22 ) 

If we are interested in solutions possessing spherical 
symmetry (i.e. independent of 9 and 0) we shall have to 
solve the simpler equation 


<^0 s i 

cr 2 r dr = 6 - cl 2 



This may bo written 


a 2 , \ d 2 , 

& (r *> = ? e? "*>■ 

showing (cf. eq. (17)) that it has solutions 

r0 = f[r-cl) -f g(r+ci), 

f and r, ag tin being arbitrary functions. Wo see, therefore, 
that there are progressive type solutions 

<f> = l ~fr-cl) + l -g(r fd) . . ( 24 ) 

Let us now turn to solutions of stationary tvpe. These 

M ‘ a >’ a11 1 ''thned by the method known as the separation 

ot variables. In one dimension we have to solve 

dty _ 1 c '4> 

Cx 2 ~ c 2 J7- ' 

Let us try to find a solution of the form 

6 = XU)T(t), 

* See e.g. Uuthorfonl, p. tiJ, equation 20. 
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X and T being functions of x and t respectively, whose 
form is still to be discovered. Substituting this value of 
<f> in the differential equation and dividing both sides 
by X(x)T(t) we obtain 


1 d-X _ 1 d-T 
X dx 1 ~ c-T dt ! 



The left-hand side is independent of /, being only a function 
of x, and the right-hand side is independent of x. Since the 
two sides are identically equal, this implies that each is 
independent both of x and t, and must therefore be constant. 
Putting this constant equal to — p 2 , wc find 

X'+p 2 X = 0, T'+c 2 p 2 T = 0. . . (26) 


These equations give, apart from arbitrary constants 



cos 

px 

sin 



cos 

sin 




A typical solution therefore is a cos px cos cpl, in which 
p is arbitrary. In this expression we could replace either 
or both of the cosines by sines, and by the principle of 
superposition the complete solution is the sum of any 
number of terms of this kind with different values of p. 

The constant — p 2 which we introduced, is known as 
the separation constant. VVe were able to introduce it in 
(25) because the variables x and t had been completely 
separated from each other and were in fact on opposite sides 
of the equation. There was no reason why the separation 
constant should have had a negative value of —p 2 except 
that this enabled us to obtain harmonic solutions (27). 
If wc had put each side of (25) equal to +p 2 , the solutions 
would have been 

X = e ipz , T = e ±f, ' f . • • (28) 

and our complete solution should therefore include terms 
of both types (27) and (28). The same distinction between 
the harmonic and exponential types of solution will occur 
frequently. 

• Thin is read : X equals cos or sin px. <*U<. 
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This method of separation of variables can be extended 
to any number of dimensions. Thus in two dimensions a 
typical solution of (IS) is 


, cos cos cos 

r ~ ■ P x ■ TJ ■ rcl , 
sin r sin sin 



in which p 2 -\-q 2 = r 2 , p and q being allowed arbitrary 
values. An alternative version of (29), in which one of 
the functions is exponential, is 



cos 

. px e ±9V 
sm 


cos 

sin 


rd. 


(30) 


in which p 2 —q 2 = r 2 . 

It is easy to see that there is a variety of forms similar 
to (30) in which one or more of the functions is altered 
from a harmonic to a hyperbolic or exponential term. 

In three dimensions we have solutions of the same typo, 
two typical examples being 



cos cos cos 

px . q y . 
sin ^ sm sin 



cos 

set , 
sin 


r + ? 2 +r 2 = fi 2 (31) 



cosh 

sinli 


pxe±w 


cos 
. rz 
sin 


cos 

. set . 
sin 


-j) 2 -g 2 +r 2 = 5 2 (32, 


1 1 lore are two other examples of solution in three 
dimensions that we shall discuss. In the first case we put 

r < 0 * ’ ^ S ' n an< * we use r > ^ an, l 2 us cylindrical 
coordinates. The equation of wave motion becomes * 


?££ lev JV lev 

ir ' r cr r 2 c6 2 ' cz 2 ~ c 2 IT 1 ' 


A solution can be found of the form 

<f> = R(r)0(d)Z{z)T(t), . . . (33) 

• See Rutherford, p. 63. 
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where, by the method of separation of variables, It, 0, Z, 
T satisfy the equations 

£? + 1 — - - R bn 2 /? = 0, 

dr 2 r dr r- 


d 2 0 

dd 2 


= —m 2 0, 


d'-Z 2/ 
d? = - q z ' 


d^T 

It 


- = -cY-T , n 2 = p 2 -q 2 • 


(34 


The only difficult equation is the first, and this* is just 
Bessel’s equation of order w», with solutions Jjnr) and 
Yjnr). J m is finite and Y m is infinite when r = 0, so 
that we shall usually require only the J m solutions, iho 
final form of (f> is therefore 

J COS n COS COS /o-v 

<& — '"(nr) . mO . qz . cpt . . (•>•>) 

r — v sin sill * sm 

•* m 

If (f> is to be single valued, in must be an integer ; but n, 
a and p may be arbitrary provided that n 2 = jr-r- 
Hyperbolic modifications of (35) are possible, similar ... 

all respects to (31) and (32). .. 

Our final solution is one in spherical polar coordinates 
r, 0 , I The equation of wave motion (22) has a solution 

/i(rj 0 ( 0 ) V (<J>) T[t), where 

• S«o Jnco, Intiqral ion o/ Ordinary OiJJercntinl Equation,. 121. 
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m, n and p are arbitrary constants, but if 'P(iJj) is to bo 
single valued, m must be integral. The first two of these 
equations present no difficulties. The 0-equation is the 
generalised Legendre’s Equation * with solution 

0(0) = P n m (cos 0), 

and if 0 is to be finite everywhere, n must be a positive 
integer. When m = 0 and n is integral, P n m (cos 0) 
reduces to a polynomial in cos 0 of degree n, known 
as the Legendre's polynomial P n (cos 0). For other 
integral values of m, P n m (cos 0) is defined by the equation 

P n "(cos6) = a in"6j^Lj^{PJcos0)). 


A few values of P„ (cos 0) and P n m (cos 0) are given 
below, for small integral values of n and m. When m>n, 
P n ' n (cos 0) vanishes identically. 

Pq (cos 0) = 1 

/’, (cos 0) = cos 0 

(co.s 0) = i (3 cos 2 0—1) 

/’ 3 (cos 0) = j (3 cos 3 0—3 cos 0) 

P t (cos 0) = l (33 cos* 0-30 cos 2 0+3) 

/V (cos 0) = sin 0 
/V (cos0) = 3 sin 0 cos 0 
/y (cos 0) = j sin 0 (5 cos 2 0—1) 
p? (cos 0) = 3 sin 2 0. 


To solve the /{-equation put Ii(r) = r~^ 2 S(r), and we 
find that the equation for S(r) is just Bessel’s equation 

dr 2 1 rdr \ P r 2 
Therefore S(r) = J nUI2 (pr) or T n+1 2 (pr). 

* See I nee. 1 nt.j, n j Ordinary Differential Equations, p. 119, 
for the cilso m = i) 
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Collecting the various terms, the complete solution, apart 
from hyperbolic modifications, is seen to he 

* = y" +1B {pr)P fl m (cos 0) 008 “ 8 cpt. (36) 

• r n + l/2 Sin r SU1 r ' 

If has axial symmetry, we must only take functions 
with m = 0, and if it has spherical symmetry, terms with 
w — w = o. Now J V2 (z) = y/(2lnz) sin z, and also 
'1/2(2) = — V(2/tt2) cos z, so that this becomes 

, , cos cos 

* = f sin V sin C P L • • (37) 

A solution finite at the origin is obtained by omitting the 
cos pr term. b 

§ 8. We shall now gather together for future reference 
U 80 "f ,ons obtained in the preceding pages. 

Progressive waves 

1 dimension 

_ l <*}> 
dz* ~ c 2 it 2 

<f> =f(x-(U)+g:x-{(U) . . (17) 

2 dimensions 

_ 1 

d,j * ~ c s Hi 

<P -ffa+my-ty+giU + my+ct), l 2 +m*= 1 . (19) 

3 dimensions 

_u*4> 

. te* dy t+ cz* c * i t* 

* - , nix +my + ( nz - ( ay -\g[U +my + nz-\ c t), P \ m 2 + n 2 = I 

3 dimensions, spherical symmetry ’ 

20^ __ 1 ^ j • j 

<** r6r ~ c l I7* ’ < P = ;f^ r cl) + -g(r+cl) . (24) 
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Stationary waves 
1 dimension 


8*$ 
8x 2 


c 2 a a ’ 


, cos cos 

* = Bin ^ sin °P‘ ■ < 27 > 


<f> = e 


±Px P ±rpl 


(28) 


2 dimensions 


8*<f> 8 2 (f> 1 c n -(f> 

d? + W ~ ?Tr- ’ 


, COS cos cos 

^ = sin sin sin rd -^ ^' =r ' 

(29) 

, COS , cos A 

^ = ain P xe < 1 \[n rct ' P—r= r- 

(30) 


3 dimensions 


o uiutuiiaiuiis COS COS COS COS 

< f , = s\nP x sin** sin” .in**' 

C 3 + C 3 + 5 = 15J P 2 +? 2 +r*=« 2 . (31) 

c'j : 2 6y 2 cte 2 c 2 ci 2 ’ . cosh , cos cos 

t = sinli P*' sin r2 sin “*• 


l -? s -?Hr J =i 2 . (32) 

Plane Polar Coordinates (r, 0) 

1 8<ft 1 8 2 tf» _ 1 c 2 tf> J m cos cos 

cr 2 r rcr ‘ r*d$* ~ c 2 a 2 ’ * “ 7 m (nr) sin ,nd sin cnL 

(35a) 

Cylindrical polar coordinates (r, 0, z) 

<?± +l?tj.L e V , _ 1 di <f> 

dr 2 T r cr ' r 2 £0 2 + c: a ~ c 2 t/ 2 ’ 
j _ , . COS COS CoS 

* - Y m f nr > sin sin sin ^ P'-? ! ■ (35b) 


Spherical Polar Coordinates (r, 0, 0) 

^ + ?a + _2_ »/. , i av 

£ r ' r 1 3 0 £0 V 9 80 j + r 2 sin 2 0 80r 


<f> = r 


-1/2 ^n-fl /2 


r "" iI (i.r)P„-(cos«) COS m ^ COS c P i 

■*n + l ( 2 sm 7 Sin ^ 


c 2 £t 2 ’ 
(36) 
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Spherical symmetry 
d 2 (f> 2 d<p _ 1 djf> 

fr* + r ~dr ~ c 2 dt 2 ’ 


<f> = r~' 


cos cos 

pr . cpt . 
sin sin 



It should be noted that there are exponential modifications 
of all the above equations. 

In solving problems, we shall more often require 
progressive type solutions in cases where the variables x, 
y, z are allowed an infinite range of values, and stationary 
type solutions when their allowed range is finite. 


§ 9. There is an important modification of the equation 
of wave motion which arises when friction, or some other 
dissipative force, produces a damping. The damping effect 
is usually allowed for (see e g. Chapter II and elsewhere) 

bv a term of the form k^$, which will arise when the 

damping force is proportional to the velocity of the vibra- 
tions. The revised form of the fundamental equation, 
known as the equation of telegraphy, is 


2 i 1 

V 2 <£= Ill^+ A: 77 


(38) 


^21 

If we omit the term -? this equation is the same ns that 

vl 

occurring in the flow of heat. If we put <f> = ue t " 2 , wo 
obtain an equation for u of the form 


1 ( c 2 u 1..,) 

r • ' 


(39) 


Very often k is so small that wo may neglect k 2 , and then 
(39) is in the standard form which we have discussed in 
§ 8, and the solutions given there will apply. In such a 
case tho presence of the dissipative term is shown by a 
decay factor e"»*. If this is written in the form e"«S 
then < 0 (= 2/k) is called the modulus of decay. When 
the term in k 2 may not bo neglected, wo have to solve 
(38) and tho method of separation of variables usually 
enables a satisfactory solution to bo obtained without 
much difficulty. 
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There is an alternative solution to the equation of 
telegraphy that is sometimes useful. Taking the case of 
one dimension, and supposing that k is so small that 
k- may be neglected, we have shown that the solution of 
(38) may be written in the form 

<f> = e~ m J(x-d) (40) 

where / is any function. Since / is arbitrary, we can put 

f(x-cl) = e rt> g(x-ci), 

and g is now an arbitrary function. Substituting this in 
(40) we get 

( f> = e- kx:Sf g(i- C t). . . . (41) 

This expression resembles (40) except that the exponential 
factor varies with x instead of with t. 


§ 10. Most of the waves with which we shall be 
concerned in later chapters will bo harmonic. This is 
partly because, as we have seen in § 8, harmonic functions 
arise very naturally when we try to solve the equation of 
wave motion ; it is also due to the fact that by means of 
a Fourier analysis, any function may be split into harmonic 
components, and hence by the principle of superposition, 

any wave may be regarded as the resultant of a set of 
harmonic waves. 

M lii'ii dealing with progressive waves of harmonic type 
then' is one simplification that is often useful and which 
is especially important in the electromagnetic theory of 
'glit waves. We have seen in (11) that a progressive 
larmonic wave in one dimension can bo represented bv 
<l> ~a cos 2 t T (kx-nt). If wo allow for a phase e , it 

H ‘ Wntten < f , = a cos {Mkx-nt)+ € }. Now this 
latter function may be regarded as the real part of the 
complex q,„„,,i,v„ j t u most ' convcnient 

or our 8u »< on-nt work il we choose the minus sign and 
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also absorb the phase « and the amplitude a into one 
complex number A. We shall then write 


A = A e iiriw ~ kx) , A = a e~ u . 


. ( 42 ) 


This complex quantity is itself a solution of the equation 
of wave motion, as can easily be seen by substitution, and 
consequently both its real and imaginary parts are also 
solutions. Since all our equations in <f> are linear, it is 
possible to use (42) itself as a solution ol the equation 
of wave motion, instead of its real part. In any equation 
in which 6 appears to the first degree, we can, if we wish, 
use the function (42) and assume that we always refer to 
the real part, or we can just use (42) as it stands, without 
reference to its real or imaginary parts. In such a case the 
apparent amplitude A is usually complex, and since 
A== a e~ if , we can say that \A\ is the true amplitude, 
and -ary A is the true phase. The velocity, of course, 

as given by (7) and (10), is njk. 

We can extend this representation of <f> to cover waves 

travelling in the opposite direction by using in such a case 


<f> = A e 2irilnt 4 tr) • 


. (43) 


There is obviously no reason why we should not extend 
this to two or three dimensions. For instance, in three 

dimensions , 

(f, = A e ^ nt ~ { P z+n ^ r,)) ■ • • (44) 

would represent a harmonic wave with amplitude A 
moving with velocity n/Vd^W+r 2 ) in the direction 
x : y : z = p : q : r. 

§ 11. We shall conclude this chapter with an example. 

Let us find a solution of ^ ^ 8Uch that ^ 

vanishes on tho lines x = 0, x = a, y — 0, y — b. Since 

B 
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the lines x = 0, a, and y = 0, b are nodal lines, our 
solution must be of the stationary type. Referring to § 8, 
equation (29), we see that possible solutions are 


. cos cos cos , „ . „ 

<p = . px . ay . ret . where p 2 +<7 2 = 
T sin^ sin sin r 



Since (f> is identically zero at x = 0, and y = 0, we shall 
have to take the sine rather than the cosine in the first 
two factors. Further, since at x = a, <{> = 0 for all 
values of y, therefore 

sin pa = 0. 

Similarly, sin qb = 0. 

Hence p = ninja, and q = rnrjb , in and n being integers 
A solution satisfying all the conditions is therefore 


whore 


, . WlTlX . 72771/ COS 

ffl = sin sm — - . ret 

r a b srn rct • 

r 2 = 7r 2 (»i 2 /a 2 -|-7i 2 /6 2 ). 


The most general solution is the sum of an arbitrary 
number of such terms, e.g. 


. * m7TX • W7 TV _ 

+ == m % sm v 81,1 T {Cm " cos rcl +D mn sin ret). 


(45) 


At l = 0, this gives 


[<R-o = rC mn sin ^ sin 

b 


a 


= rrcZ) mn sin— sin 

a b 


we can 


By suitable choice of the constants C mn and D mn 

inake and 0 have any chosen form at t = 0. The value 
at any subsequent time is then given by (45) 
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§ 12. Examples 

(1) Show that <f> = /(x cos 0+y sin 0-cl) represents a 
wave in two dimensions, the direction of propagation making 
an angle 6 with the axis of x. 

(2) Show that <f> = a cos (Ix+my—cl) is a wave in two 
dimensions and find its wavelength. 

(3) What is the wavelength and velocity of the system of 
plane waves <f> = a sin (Ax+By+Cz — Dt) 1 

(4) Show that three equivalent harmonic waves with 120° 
phase between each pair have zero sum. 

(5) Show that 4 = r-V 2 cos \9J(r±ct) is a progressive typo 
wave in two dimensions, r and 6 being plane polar coordinates, 
and / being an arbitrary function. By superposing two of 
these waves in which / is a harmonic function, obtain a 
stationary wave, and draw its nodal lines. Note that tine is 
not a single-valued Junction unless we put restrictions upon the 

allowed range oj 6. 

(6) By taking the special case of /(*) = g(x) = sin px in 
equation (24). show that it reduces to the result of equation 
(36) in which m = n = 0. Use the relation 




- V (I) sin *■ 


(7) Find a solution of + ~ x - °. such tl,at * 0 


when t 


dx x c 

co, and <f> = 0 when x = 0. 

d*<f> 

C£ % C 

4-co or J = + co. 


(8) Find a eolation of g - i g .hot * - 0 when 


(8) Solve the equation g - O g given that . is never 

infinite for real values of x and t , and z = 0 when x = 0, or 
when t = 0. 

(10) Solve % = % given that V - 0 whan < - oo and 
' ' cx 1 ot 

when x = 0, and when x => (. 
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( 1 1 ) x, y , z are given in terms of the tliroe quantities 
£, t /, £ by the equations 

x = a sinh £ sin 77 cos £ 
y = a sinh £ sin 77 sin £ 

2 = a cosh £ cos 77 


ci 4X 4 4l . &4> d 2 <f> 1 W. , t 

Show that the equation — -) -| = is of the 

dx* dy* dz* c 1 dl l 


correct form for solution by the method of separation of 
variables, when £, 77, £ are used as the independent variables. 
Write down the subsidiary equations into which the whole 


equation breaks down. 


12. Show that the equation of telegraphy (38) in one 
dimension has solutions of the form 

cos cos 

* = sin " lX sin P lerkt *• 

whoro m and p aro constants satisfying the equation 
P* = mV*- {lc-. 


[Answers: 2. 2.t/(.* + m‘)l ; 3. A = 2ttI(A* + B* + C*)l, 

vol- = AZA/2rr ; 7. A sin nx e ~ cnt ; 8. 9. A sin px 

sin cpt; 10. Ae~ * sin px, p = */l, 2*11. . 11. Show 

that £ — const., rj = coast., £ = coast, form an orthogonal 

system of coordinates, and transform y 1 * in terms of (, q, £ 

88 ,n Rutherford, Vector Methods, §47. The result is 

^ — ( j j) 2(£)7'(/), whero m, p and q are arbitrary con- 
stunts, and 


m 2 


1 (/ ■ , . dX 

Stahl if 8mh 1 7l - 5tai^ X + P* X - 3<X. 


1 


m- 


ST* ^ ^ y = -9»r. 

f Z _ , t7 d-T c'p * 

— — m-Z. — //» i 

cfc* a* * J 


dz 1 


CHAPTER II 


WAVES ON STRINGS 

§ 13. In this chapter we shall discuss the transverse 
vibrations of a heavy string of mass p per unit length. l>y 
transverse vibrations we mean vibrations in which the 
displacement of each particle of the string is in a direction 
perpendicular to the length. When the displacement is 
in the same direction as the string, we call the waves 
longitudinal ; these waves will be discussed in Ghaptei IN • 
We shall neglect the effect of gravity; in practice this 
may be npproximitely achieved by supposing that the 
whole motion takes place on a smooth horizontal pi me. 

In order that a wave may travel along the string, 
it is necessary that the string should bo at least slightly 
extensible ; in our calculations, however, wo shall assume 
that the tension docs not change appreciably from its 
normal value F . The condition for this (see § 14) is that 

the wave disturbance is not too large. 

Let us consider the motion of a small element of the 
string PQ (fig. 1) of length ds. Suppose that in the 
equilibrium state the string lies along the axis of x, and 
that PQ is originally at P 0 Q 0 . Let the displacement of 
PQ from the a; axis be denoted by y. Then wo shall obtain 
an equation for the motion of PQ in terms of the tension 
and density of the string. The forces acting on this 
element, when the string is vibrating, are merely the two 
tensions F acting along the tangents at P and Q as shown 
in the figure ; let 0 and </>+</</> bo the angles made by these 
two tangents with the x axis. We can easily write down 
21 
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the equation of motion of the element PQ in the y direction; 
for the resultant force acting parallel to the y axis is 
F sin (tp-\-d>p)— F sin tp. Neglecting squares of small 
quantities, this is F cos tp dtp. The equation of motion is 
therefore 



Fio. 1 


Now tan tp 
from (1) 


fa* 30 l * iat sec-tpdip = C ^d.c, and so, 


c x y „ chj dx 
P ~=Fco^-.- 

— F cos *tp — 

Y dx* 


. . ( 2 ) 


But cos *tp = 


w! * 


so that if the displacements 
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are small enough for us to neglect 



compared with 


unity, we may write (2) in the standard form for wavo 
motion* (Chapter I, § 5), viz., 


' c l± 

8x 2 


1 '2 

— , where c 2 = F/p 

c 2 ct 2 ’ Ir 



It follows from Chapter I, equation (17) that the general 
solution of this equation may be put in the form 

y = f(x-ct)+g(x+ct), . . . (•*) 


/ and g being arbitrary functions. f(x—ct) represents a 
progressive wave travelling in the positive direction of 
the x axis with velocity c, and g(x+ct) represents a 
progressive wave with the same velocity in the negative 
direction of x. Thus waves of any shape can travel in 
either direction with velocity c=^(F/p) t and without 
change of form. A more complete discussion, in which 
we did not neglect terms of the second order, would show 
us that the velocity was not quite independent of the 
shape, and indeed, that the wave profile would change 
slowly with the time. These corrections are difficult to 
apply, and we shall be content with (4), which is, indeed, 
an excellent approximation except where there is a sudden 

, IW 

“ kink ” in y, in which case we cannot neglect I • 


§ 14. Since the velocity of any point of the string is y, 
wo can soon determine the kinetic energy of vibration. It is 


T=j)tpy 2 dx .... ( 5 ) 

• The student who is interest**! in geometry will bo able to 
prove that the two tensions fit P find Q are together equivalent to a 
singlo force of magnitude Fds/R, where R is the radius of curvature 
of the string. This force acts perpendicularly to PQ. I fitting 


R =» { 1 + (fi) 2 } 3 ~/ ttlld n °e ,octin 8 (“) » wo obuiw 
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Tie potential energy V is found by considering the increase 
of length of the element PQ. This element has increased 
its length from dx to ds. We have therefore done an 
amount of work F(ds—dx). Summing for all the elements 
of the string, we obtain the formula 



F.di i - dr) 




approximately. 



The integrations in (.3) and (6) are both taken over the 
length of the string. 

With a progressive wave y = f(x—cl), these equations 
give 


T =jlpc'-(fr‘dr = IF . . (7) 

V = ( 8 ) 


Thus the kinetic and potential energies are equal. The 
same result applies to the progressive wave y = g(x+ct), 
but it does not, in general, apply to the stationary type 
waves y=f(x-ct)+g(x+ct). 

We can now decide whether our initial assumption is 
correct, that the tension remains effectively constant. 
It the string is elastic, the change in tension will be pro- 
portional to the change in length. Wo have seen in (6) 


that the change in length of an element dx is - l—Ydx. 

2 \dx) 

Thus, provided that c ± is of the first order of small 

quantities, the change of tension is of the second order, 
aiu nu} sa o j e neglected. This assumption is equivalent 
to asserting that the wave profile does not have any largo 
kinks, but has a relatively gradual variation with z. 
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§ 15. The functions / and g of (4) are arbitrary. But 
they may be fixed by a knowledge of the initial conditions. 
Thus, with a string of unlimited length, such that 

y t - o = #*). vt - o = fa) * we raust have - from {4)> 

fa)+9( x ) = fa)< 

-cf'(x)+cg'(x) = fa). 


Integrating this last equation we have 

-fa) +9(*) = (1/c) J fa fa, 


and so 


fa) = \ {fa)~- c \ * {x)dx 

g[x) = ^ ^fa) + ~ e j fa) dx ^j- 

The displacement at any subsequent time l is therefore 

y = \ {fa-ci) + fa+ci) - fa) dx + cj * (x)dj: j 

= i ! fa-ci) + fa+cl) + fafa j 


(9) 


§ 16. The discussion above applies specifically to 
strings of infinite length. Before we discuss strings ot 
finite length, we shall solve two problems of reflection 
of waves from a discontinuity in tlio string, ilw rs is 
when two strings of different densities are joined toget ier, 
and the second is when a mass is concentrated at a I K,in 
of the string. In each case wo shall find that an inci< on 
wave gives rise to a reflected and a transmitted wave 

Consider first, then, the case of two sonii-m finite- 
strings 1 and 2 joined at the origin (fig. 2). bet tlio 

* This function Hz) must he distinguished from the unglo i m 

5 13. 
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densities of the two strings be p l and p 2 . Denote the dis- 
placements in the two strings by y, and y 2 . Let us suppose 
that a train of harmonic waves is incident from the negative 
direction of x. When these waves meet the change of 
wire, they will suffer partial reflection and partial trans- 
mission. If wo choose the exponential functions of § 10 
to represent each of these waves, we may write 

Vl ~ 2/incident 2/ reflected 

J/a = J/traiumltted • • 

where 

2/incident = !<**«-*# 

2/rcflcctcd = B l e- ni{nt +i,x) 

2/iraiisiiiitted = . . ( 11 ) 



is real, but A» and B t may be complex. According to 
§ 10 equation (42), the arguments of A s and B, will give 
tlieir phases relative to the incident wave. All three waves 
in (11) must have the same frequency n, but since the 
velocities in the two wires are different, they will have 
different wavelengths 1/it, and l/it 2 . The reflected wave 
must, of course, have the same wavelength as the incident 
"ave. Since the velocities of the two types of waves are 
and njk 2 (Chapter I, equations (7) and (10)), and we 
have shown in (3) that c* = Flp, therefore 

Pi Pi • 


. ( 12 ) 
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In order to determine A s and B 1 we use what are known 
as the boundary conditions. These are the conditions 
which must hold at the boundary point x=0. Since the two 
strings are continuous, we must have y l = y 2 identically 
for all values of t, and also the two slopes must be the same, 


that = ^? for all t. 
dx dx 


If this latter condition were 


not satisfied, we should have a finite force acting on an 
infinitesimal piece of wire at the common point, thus 
giving it infinite acceleration. Wo shall often meet 
boundary conditions in other parts of this book ; their 
precise form will depend of course upon the particular 
problem under discussion. In the present case, the two 
boundary conditions give 


A x A-B x = Ao, 

2w*(-Mi+*A) = 27Ti(-M«)- 

These equations have a solution 


k,+k 2 ’ A, ~ k x +k 2 ' • {,6) 

Since k lt k 2 and A x are real, this shows that B l and A 2 
are both real. A % is positive for all k x and k 2 , but B x is 
positive if k x >k 2 , and negative if k x <k 2 . Thus the 
transmitted wave is always in phase with the incident 
wave, but the reflected wave is in phase only when the 
incident wave is in the denser medium ; otherwise it is 
exactly out of phase. 

The coefficient of reflection R is defined to be the 


ratio, \BJAtf, i.e. which, by (12), we may write 



Since, from (7) and (8), the energy of a progressive 
wave is proportional to the square of its amplitude, it 
follows that R represents the ratio of reflected energy to 
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incident energy. Similarly, since no energy is wasted, the 
coefficient of transmission T . which gives the ratio of 
transmitted energy to incident energy, is equal to 1— R, 

WpiVpz 


l.e. 


(Vpi+VPz ) 2 


(15) 


§ 17. A similar discussion can be given for the case of 
a mass .1/ concentrated at a point of the string. Let us 
take the equilibrium position of the mass to be the origin 
(fig. 3) and suppose that the string is identical on the two 



sides. Then if the incident wave comes from the negative 
side of the origin, we may write, just as in (10) and (11) : 

l/l = .'/incident d" 1/ reflected 
Vi ~ Z/transmltted 

where 

incident = 

2/ reflected = 

•/transmitted = A 2 e'-" ,,n<-fcr) . . . (16) 


1 be boundary conditions art* that for all values of t 
(*) lyjlr-o = (y 2 )x-o .... 



(17) 

(18) 
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The first equation expresses the continuity of the string 
and the second is the equation of motion of the mass M. 
We can see this as follows : the net force on M is the 
difference of the components of F on either side, so that 
if i/q and < fj 2 are the angles made with the x axis, we have 

M ^ 2 ] =F( sin 0,-sin 0 t ). 

Since 0, and 0, are small, we may put sin i Ji 2 = 

tan 02 = ^ . 8^ 0. = and < 18 > is thon obtilin( ‘ d - 

Substituting from (16) into (17) and (18), and cancelling 
the term er nint } which is common to both sides, we find 

+ 1+^1 = 8 ’ 

2irikF(A i —A l -\-B l ) = 4ir m n m MA t . 

Let us write irr^M/kF = p • • • • ^9) 

A solution of the equations then gives 

_ -*? _ -y*-»g . . (20) 

A x 1+t'p 1+p 2 

d 2 = _!_ = . . . (21) 

A x 1 +ip 1+P 2 

In this problem, unlike the last, B 1 and A % are complex, 
so that there are phase changes. These phases (according 
to §10) are given by the arguments of (20) and (21). 
They are therefore tan -1 (p) and tan -, ( — 1/p) respective!} . 
The coefficient of reflection R is | BjAtf, which equals 
p*/(l+p 2 ), and the coefficient of transmission T is 
1-i?, i.e. 1/(1 +J> 2 )- If wc wrile p = Ian 0, where 
O< 0 < 7 t/ 2 , then we find that the phase changes are U 
and tt/ 2 + 0, and also R = sin 2 0, T = cos 2 0 ■ 


§ 18. The two problems in §§ 16, 17 could In* solved 
quite easily by taking a real form for vach <>! the wuws 
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instead of the complex forms (11) and (16). The student 
is advised to solve these problems in this way, taking, for 
example, in § 17, the forms 

^Incident = a l COS 2tr(nl-kx) 

^reflected = *>i cos {2v[nt+kx) +e} 

2/tranamltUd = a 3 COs{2ir(nt-kx)+7 ) } . (22) 

In most cases of progressive waves, however, the complex 
form is the easier to handle ; the reason for this is that 
exponentials are simpler than harmonic functions, and 
also the amplitude and phase are represented by one 
complex quantity rather than by two separate terms. 


§ 19. So far we have been dealing with strings of 
infinite length. When we deal with strings of finite length 
it is easier to use stationary type waves instead of progres- 
sive type. Let us now consider waves on a string of length 
l, fastened at the ends where x = 0,1. We have to find a 


solution of the equation (3), viz. ^~ = \ ^ subject to 

ox* c- or 

the boundary conditions y = 0, at x = 0, l, for all t. 
Now by Chapter I, § 8, we see that suitable solutions are 
of the typo 


cos cos 
Bin P x sin C P‘- 


It is clear that the cosine term in x will not satisfy the 
boundary condition at x = 0, and we may therefore write 
the solution 

y = sin px (a cos cpt-\-b sin cpt). 


The constants a, b and p are arbitrary, but we have still 
to mako y = 0 at x = l. This implies that sin pi = 0, 
i.o. pi = 77 , 277, 3:7. . . . It follows that the solution is 


rrr.r / net 
sin ~y I o cos — 


j- + b Bin , r= 1,2.3, . . . (23) 



WAVES ON STRINGS 


31 


Each of the solutions (23) in which r may have any positive 
integral value, is known as a normal mode of vibration. 
The most general solution is the sum of any number ot 
terms similar to (23) and may therefore be written 

y = E sin^Jfl,. cos ^ +& r sin -y j • < 24 ) 

The values of a r and b T are determined from the initial 
conditions ; thus, when t = 0, 


Vt - o = ^ a r 8in 


tttx 


_ l 7*7 tC . rrrx 
tit - o = Z b r - y-sin j 


(23) 

( 20 ) 


If we are told the initial velocity and shape of the string, 
then each a, and 6, is found from (25) and (26), and hence 
the full solution is obtained. We shall write down the 

results for reference. If we suppose that when 1- 0. 
y = y = </>(*), then the Fourier analysis represented 

by (25) and (26) gives 


2 f ' . rirx . 

a, | I -j dx 

9 f ' . rirx 

b. = — I </'(*) s,n T d£ 

rnc J o 1 


(27) 


In particular, if the string is released from rest when 
t = 0, every b r = 0. 

§ 20. As an illustration of the theory of the last section 
let us consider the case of a plucked string of length 1 
released from rest when the midpoint is drawn i aside 
through a distance h (fig. 4). In accordance with (25) 
and (26) we can assume that 


y = Za r sin -j 


rrrx rvd 

cos — 
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TTTX 


When t = 0, this reduces to £a r sin — , and the coefficients 

f * 

a T have to be chosen so that this is identical with 


y = 
y = 


V x , O^x^U 

l • 

j (l-x) , 



If we multiply both sides of the equation y = £a r sin 

r 


rnx 

T 


by sin ™ , and integrate from x = 0 to x = /, as in the 

method of Fourier analysis, all the terms except one will 
disappear on the right-hand side, and we shall obtain 


l 



Whence 


r‘-2h . TTTX f' 2h f t . r7T.r 

J 0 Y" a, "T <i ^J„ 2 T || - I > 8m T 


8 h . TTT 

— sm 7 T when r is odd, 
TT-r- 2 


= 0 when r is even. 

So the full solution, giving the value of y at all subsequent 
times, is 


X 


8h „ 
!/ = ->£ 


1 


. (2/i (-I)tt . (2nfl)7rx (2n + l)7TC/ 

ir* „r„ I2n + 1 j* * 2 8,i ‘ j C ° B 7 

(28) 
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Thus the value of y is the result of superposing certam 
normal modes with their appropriate amplitudes. These 
are known as the partial amplitudes. The parted 
amplitude of any selected normal mode (the rth fo 
example), is just the coefficient a,. In tins example, , 
vanishes except when r is odd, and then a is projwrtional 
to 1 1,\ so that the amplitude of the higher modes is 

relatively small. 

S 21. The rth normal mode (23) has a frequency rc/21 

with the finger lightly resting on the P° ,nt i/r ® 

found that this mode of vibration is cxc.tcd. With c vc 
order vibrations (r even) the mid-point is a node, and 

odd order vibrations it is an antmode 

We can find the energy associated with t t » * 

vibration most conveniently by rewriting (- ) 


nr* \ rvet 

y = A sin — cos I — + « 


. . (29) 


Here A is the amplitude and c is the phase. According 
to (5) the kinetic energy is 

_ 1 f' .. , nWp At i . (30) 

T = 2 P )/ dx= ~4T l I / 

Similarly, by (G) the potential energy is 




* cos ! ! 


1 l 


-+e\. (3D 


Now by (3) Flp = c* and 

and (31) are equal. Inc total tn 0 > 

ia therefore 

irVr»p A2 ' m . (32) 

Al 


o 
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The total energy is thus proportional to the square of the 
amplitude and also to the square of the frequency. This 
is a result that we shall often find as wo investigate various 
types of wave motion. 

As a rule, however, there are several normal modes 
present at the same time, and we can then write the dis- 
placement (24) in the more convenient form 



Z A r sin 
r-1 


tttx 
— — cos 





A r is the amplitude, and € r is the phase, of the rth 
normal mode. When we evaluate the kinetic energy 
as in (30) we find that the “ cross-terms ” vanish, since 

I* sin ^ sin ^ dx = 0, if r # 3. Consequently the total 

Jo * } 

kinetic energy is just 


TT’C-p 

~i r 


Zr 2 A? sin* 




and in a precisely similar way the total potential energy is 


77 


-F 


4 1 


Zr 2 A 2 cos 2 


r-nct 

T 


+ C 




By addition we find that the total energy of vibration is 


7T~C~ p 

4f 




This formula is important. It shows that the total energy 
is merely the sum of the energies obtained separately for 
each normal modo. It is duo to this simple fact, which 
arises because there are no cross-terms involving AfA t , 
that the separate modes of vibration are called normal 
modes. It should bo observed that this result holds for 
both the kinetic and potential energies separately as well 
as for their sum. 
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We have already seen that when a string vibrates more 
than one mode is usually excited. The lowest frequency, 
viz. c/21, is called the ground note, or fundamental, am 
the others, with frequencies re/21, are harmonics or over- 
tones. The frequency of the fundamental vanes direct ) 
as the square root of the tension and inversely as the 
length and square root of the density. Tins is known as 
Mersenne’s law. The tone, or quality, of a vibration is 
governed by the proportion of energy in eac o 
harmonics, and it is tins that is character, st,c of each 
musical instrument. The tone must he carefully d.stm- 
guished from the pitch, which is merely the frequency of 

the fundamental. , 

We can use the results of (34) to determine the total 

energy in each normal mode of the vibrating stem® winch 

we discussed in § 20. According tn (-8) and (33) A in 

. ■ 1 « 


, 8 h 

and i4 2n +i — o 


sin 


_ (2n + l)7r Consequently, the 

„2 (2n+l)* 2 

total energy of the even mdm i. aero a„d the energy 
of the (2n + l)th mode is 10c*A 2 p/(2n + l) n i. r 

us that the main part of the energy » “ 

normal modes of low order. We can check these formula 

for the energies in this example quit. *or « 0 ‘ ol “ 

energy of the whole vibration is the sum of the energ.es 

each normal mode separately : i.e. 


IGc'h'p 


rr> 


total energy = ^T n f 0 (2n + 1 ) 2 ' 

It is shown in books on algebra that the sum of tho scries 

l/l 2 + l/3 2 + l/5*+ ... is * 2 / 8 - Hence tho ^ 1 ■ i 8 « llf i 

2c 2 k t pll ie 2 Fh*ll. But the string was drawn aside 

icleasS from rest in the posit, on of fig. 4, and a ha 
moment tho whole energy was m the form of po c nted 
energy. This potential energy is lust F tuna. th e mcre.se 

in length, i...2F{(!*/4 +**>'«-«■ A s.mplc ealculatum 
■hows that if wo neglect powers of h above tho second, 
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as we have already done in our formulation of the equation 
of wave motion, this becomes 2 Fhrjl, thus verifying our 
earlier result. 

This particular example corresponds quite closely to 
the case of a violin string bowed at its mid-point. A listener 
would thus hear not only the fundamental, but also a 
variety of other frequencies, simply related to the funda- 
mental numerically. This would not therefore be a pure 
note, though the small amount of the higher harmonics 
makes it much purer than that of many musical instru- 
ments, particularly a piano. 

If the string had been bowed at some other point than 
its centre, the partial amplitudes would have been different, 
and thus the tone would be changed. By choosing the 
point properly any desired harmonic may bo emphasised 
or diminished, a fact well known to musicians. 

§ 22. We have seen in § 21 that it is most convenient 
to analyse the motion of a string of finite length in terms 
of its normal modes. According t-o (33) the rth mode is 

. kttx rrret , 

V , = A, Sill — cos — + €f . 

We often write this 

y r = <f> r sin ^ . . (35) 

The expressions </> r are known as the normal coordinates 
for t he string. There are an infinite number of these 
coordinates, since there are an infinite number of degrees 
of freedom in a vibrating string. The advantage of using 
these coordinates can be seen from (30) and (31) ; if the 
displacement of the string is 

V- i:4> r sin ~ . . (36) 

r-1 * 
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then 


T = £ plZir 2 

4 1 r 


(37) 


The reason why we call <f> r a normal coordinate is 
now clear; for in mechanics the normal coordinates 
q., q 2 ...q n are suitable combinations of the original 
variables so that the kinetic and potential energies can be 
written in the form 


T = a l q*+a 2 q i *+a 3 q 3 *+ . . • 

V = b l q l i +b 1 q i *+b 3 q 3 t + 

The similarity between (37) and (38) is obvious Further, 
it can be shown, though we shall not reproduce the analyse, 
here, that Lagrange’s equations of motion app y " ‘ 1 
set of coordinates <f>, in just the same way as wi 
coordinates q r in ordinary mechanics. 


§ 23. We shall next discuss the normal modes of a 
string of length l when a mass M is tied to its mi P°' 
(fig. 5). Now we have already seen in § 21 that 



normal vibrations of an unloaded string the normal mil ■ 
of even order have a node at the m.d-point. In such a 
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vibration there is no motion at this point, and it is clearly 
irrelevant whether there is or is not a mass concentrated 
there. Accordingly, the normal modes of even order are 
unaffected by the presence of the mass, and our discussion 
will apply to the odd normal modes. 

Just as in the calculations of §§ 16, 17, in which there 
was a discontinuity in the string, we shall have two 
separate expressions y l and y 2 valid in the regions 0^x^i/2 
and l/2<x^l. It is obvious that the two expressions 
must be such that y is symmetrical about the mid-point 
of the string, y, must vanish at x = 0 and y 3 at x — l. 
Consequently, we may try the solutions 

y, = a sin px cos ( cpt -f e) 

y 2 = a sin p{l—x) cos (cpt-fc) . . (39) 


We have already satisfied the boundary condition y, = y 2 
at x = 1/2. There is still the other boundary condition 
which arises from the motion of J/. Just as in (18) wo 
may write this 


F 






Substituting the values of y, and y a as given by (39) and 
using the relation F = c~p, we find 


pi . pi pi 
2 *“ 2 = J? “ CO " St ' 



The quaitity plj 2 is therefore any one of the roots of the 
equation x tan x = plj M . If we draw the curves y = tan x, 
y = pi Mx, wo can see that these roots lie in the regions 
0 to it / 2 , 7 t to 377/2, 277 to 077/2, etc. If we call the roots 
x i» *2 ••• then the frequencies cp / 277 become cxJttI. If M 
is zero so that the string is unloaded, x r = (r + 1/2)77, 
so tin* presence of M has the effect of decreasing the 
frequencies of odd order. 

It we v. n for the frequency of a normal mode, 
then, sir, 1 -.= cpl'lit, it follows that (40) can bo written 
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in the form of an equation to determine n directly ; viz., 

x tan x = pl/M, where x = (irllc)n . • ( 41 ) 

This equation is called the period equation. Its solutions 
are the various permitted frequencies (and hence periods) 
of the normal modes. Period equations occur very fre- 
quently, especially when we have stationary type waves, 
and we shall often meet them in later chapters, ilus 
particular period equation is a transcendental equation 
with an infinite number of roots. 


§ 24. In the previous paragraphs we have assumed tha 
there was no frictional resistance, so that the vibrations were 
undamped. In practice, however, the air does l ,r0V1< 1 
a resistance to motion ; this is roughly proportional to the 
velocity. Let us therefore discuss the motion of a string 
of length l fixed at its ends but subject to a resistance 
proportional to the velocity. The fundamental equation 
of wave motion (3) has to be supplemented by a term in 


d l 

dt 


and it becomes 


dx 2 




A solution by the method of separation of variables (cf. § 9) 

is easily made, and we find __ 

y = Ae~^ u sin pz cos (V(c 2 p 2 —k 2 l4)t-\-c)- 

Since y is to vanish at the two ends, we niust 
before, sin pi = 0, and hence p = r W » r * 

The normal modes of vibration are there ore 


y = A r e~ ikt sin ^ cos (?<+€,) • • ( 43) 

where 

rVc 5 k- 
? = — “ 4 ' 
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The exponential term represents a decaying amplitude 
with modulus (see § 9) equal to 2/jfc. The frequency 
ql2ir is slightly less than when there is no frictional resist- 
ance. However, k is usually small, so that this decrease 
in frequency is often so small that it may be neglected. 


§ 25. There is another interesting method of obtaining 
the velocity of propagation of waves along a string, which 
we shall now describe and which is known as the method 
of reduction to a steady wave. Suppose that a wave is 
moving from left to right in fig. 6 with velocity c. Then, 



if we superimpose on the whole motion a uniform velocity 
— c the wave profile itself will be reduced to rest, and 
the string will everywhere be moving with velocity c, 
keeping all the time to a fixed curvo (the wave profile). 
We are thus led to a different problem from our original 
•me ; for now the string is moving and the wave profile 
is at rest, whereas originally the wave profile was moving 
and the string as a whole was at rest. Consider the motion 
of the small element PQ of length ds situated at the top 
of the hump of a wave. If R is the radius of curvature at 


this toj) point, and we suppose, as in § 13, that the string 
is almost inextonsible, then the acceleration of the element 
I Q is c 2 h downwards. Consequently, the forces acting on 
it must reduce to (c 2 , 11) p<ls. But these forces are merely the 
two tension F at P and Q, and just as in § 13 (especially 
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note at foot of p. 23 ), they give a resultant Fds’R downwards 
Equating the two expressions, we have 

C ^i* = F% i.e. c‘ = F/p. 

R n 

This is, naturally, the same result as found before. The 
disadvantage of this method is that it does not descr.be 
in detail the propagation of the wave, nor does it deal 
with stationary waves, so that we cannot use it to get 
the equation of wave motion, etc. It is, however very 
useful if we are only concerned with the wave velocity, 
and we shall see later that this simple artifice of reducing the 
wave to rest can bo used in other problems as well. 


§ 26 . 


Examples 


(1) Find the velocity of waves along a string whose density 
is 4 gms. per cm. when stretched to a teasion 90000 dynes. 

(2) A string of unlimited length is pulled into a harmonic 
shape y = a cos fcx, and at time « - 0 it is released. Show 
that if F is the teasion and p the dens.ty of the string its 
shape at any subsequent time * is y = a cos kx cos ket, where 
c* = F/p. Find the mean kinetic and potential energies per 

unit length of string. . 

(3) Find the reflection coefficient for two strings which 

are joined together and whoso deasities are 25 gms. per cm. 

and 9 gms. per cm. . n , tin , 

(4) An infinite string lies along the x axis. At t = 0 that 

part of it between x = ± a is given a transverse velocity 
«* -x*. Describe, with the help of equation (9) the subsequent 
motion of the string, the veloc ity of wave motion i being ! c. 

(5) Investigate the same problem as in question ( ) P 
that the string is finite and of length 2a, fastened at the 

P °'(0) Whafta the total energy of the various normal modes 
in question (5) t Verify, by summation over all the 
modes, that this is equal to the initial kinetic en > g> • 

(7) The iwo ends of a uniform stretched string are fa tci 
to light rings that can slide freely on two fixed parallel wire 
L ring being on o«cl, „,re. the norn.ol 

modes of vibration. 
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(8) A uniform string of length 3 1 fastened at its ends, is 
plucked a distance a at a point of trisection. It is then 
released from rest. Find the energy in each of the normal 
modes and verify that the sum is indeed equal to the work 
done in plucking the string originally. 

(9) Discuss fully the period equation (41) in § 23. Show in 
particular that successive values of x approximate to rn , 
and that a closer approximation is x = rn+pl/Mrir. 

(10) Show that the total energy of vibration (43) is 
\plA*er kl {q*+kq cos ( ql+t r ) sin (qt+tj+^k 1 co>* (qt+t T )}, 

nnd hence prove that the rate of dissipation of energy is 
$k plA*er kt {2q sin (qt+t r )+Je cos (qt+tr)}*. 

(11) Two uniform wires of densities Pl and p t and of equal 
length aro fastened together at one end and the other two ends 
aro tied to two fixed points a distance 21 apart. The tension 
is F. Find the normal periods of vibration. 

(12) Tho density of a stretched string is m/x*. The end- 
points aro at x = a, 2a, and tho tension is F. Verify that the 
normal vibrations are given by the expression 

y = A sin [0 log, (x/a)](?\ ' 009 pt, where 6* = ~ - - . 

\a/ sin F 4 

Show that tho period equation is 8 log* 2 = tin, n = 1, 2, . . . . 

(13) A heavy uniform chain of length l hangs freely from 
ono end, and performs small lateral vibrations. Show that 
the normal vibrations are given by tho expression 

y = A J o(2pV{ x lv}) cos (/>*+<), 

where J 0 represents Bessel’s function (§ 7) of order zero, x 
being measured from the lower end. 

Deduce that tho period equation is J Q (2p^{lig}) = 0. 


[Answers : 

1. 150 cms./sec. ; 2. 

r ( r 4- 4)rrX 

5. y = zb r cos = — 

a 


IFaWsin* Icct, ±Fa*k* cos* feet ; 3. 1/1 G; 

. (r+J)77C/ 

s»n - , b r = ( — 1 ) r 4<i s /( r -f 1 ) 1 rr 4 c ; 


0. 8pa®/15; 7. y=a r cos ^cos ; 8. energy in rth 

. , 27 c*a*p ttt 

normal m °do = sin 3 - ; sum = 3c*a*p/4l ; 11. 2 n/p 

whore c, tan (/>//<,) = -r, tan (/»//. ,). c l * = F/ Pl , c t * = F/p,.] 
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WAVES IN MEMBRANES 

§ 27. The vibrations of a plane membrane stretched to a 
uniform tension T may be discussed in a manner very 
similar to that which we have used in Chapter II for 
strings. When wo say that the tension is T we mean that 
if a line of unit length is drawn in the surface of the 
membrane, then the material on one side of this line 


T 6* 



exerts a force T on the material on the other side, and this 
force is perpendicular to the line we have drawn. Let us 
consider the vibrations of such a membrane ; wo shall 
suppose that its thickness may bo neglected. If its 
equilibrium position is taken as the xy plane, then we are 
concerned with displacements z(xy) perpendicular to this 
plane. Consider a small rectangular element AliCD 
(fig. 7) of sides 8x, By. When this is vibrating the forces 

48 
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on it are (a) two forces T8x peqtendicular to AB and CD, 
and (6) two forces TSy perpendicular to AD and BC. 
These four forces act in the four tangent planes through 
the edges of the element. An argument precisely simjiar 
to that used in Chapter II, § 13, shows that the forces (a) 

3 2 z 

give a resultant T8x . — 8y perpendicular to the plate. 


Similarly, the forces (6) reduce to a force T8y . — 8z. Let 

8x m 

the mass of the plate be p per unit area : then, neglecting 
gravity, its equation of motion is 



8.rSy+T Sy 




This may be put in the standard form 


— A — - - — 
8x- + 'ey- c'-L *' 



where 


c a = Tip 



1 bus we have reduced our problem to the solution of the 
standard equation of wave motion, and shown that the 
velocity of waves along such membranes is c = \/(T/p). 


§ 28. Let us apply these equations to a discussion of the 
transverse vibrations of a rectangular membrane ABCD 
(fig. 8) of sides a and 6. Take AB and AD as axes of x 
and y. Then we have to solve (I) subject to certain 
boundary conditions. These are that 2 = 0 at the boundary 
-I the membrane, for all /. With our problem this means 
(hat 2 = 0 when x = i), x = a, y = 0, y = b, independent 
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of the time. The most suitable solution of the equation 
of wave motion is that of § 8, equation (29). It is 


z 


cos cos cos o i 

. px . qy .ret , p 2 +q 
sin sin sin 



If z is to vanish at a: = 0, y = 0, we shall have to reject 
the cosines in the first two factors. Further, if z vanishes 



Flo. 8 


at x = a, then sin pa = 0, so that p = mir/a, and similarly 
q = n7r/6, m and n being positive integers. Thus the 
normal modes of vibration may be written 


. mnx . nny . /qv 

Z = A 8U1 Sin — - cos ( rct+c ), . . W) 

a o 


where 


r 2 = (m 2 /a 2 -fn 2 /<r)7r 2 . 


Wc may call this the (m, n) normal mode. Its frequency 
is rc/277, i.e. 


V((S ♦ SKI ■ 


V \ • • ' 

The fundamental vibration is the (1,1) mode, for which the 


frequency ie Th ° ovortonC3 (4) a "' 

not related in any simple numerical way to the funda- 
mental, and for this reason the sound of a vibrating p ate, 
in which as a rule several modes are excited together, is 
much less musical to the car than a stiing, where tie 
harmonics are all simply related te the fundamenta . 
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In the (m, n) mode of (3) there are nodal lines x = 0, 

a/m, 2 ajm, a, and y = 0, b/n, 2b/n, ... b. On opposite 

sides of any nodal line the displacement has opposite sign. 
A few normal modes are shown in fig. 9, in which the 
shaded parts are displaced oppositely to the unshaded. 



The complete solution is the sum of any number of 
terms such as (3). with the constants chosen to give any 
assigned shape when t = 0. The method of choosing 
these constants is very similar to that of § 19 , except that 
thero are now two variables x and y instead of one, and 
consequently wo have double integrations corresponding 
to (27). 

According to (4) the frequencies of vibration depend 
upon the two variables in and n. As a result it may 
happen that there are several different modes having the 
same frequency. Thus, for a square plate, the (4, 7), 
(7, 4). (1, 8) and (8, 1) modes have the same frequency; 
and for a plate for which a = 36, the (3, 3) and (9, 1) 
modes have the same frequency. When we have two or 
more modes with the same frequency, we call it a 
cto onerate case. It is clear that any linear combination 
hose modes gives another vibration with the same 
. .jjency. 
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§ 29. We can introduce normal coordinates as in the 
case of a vibrating string (cf. § 22). According to (3) 


the full expression for z is 



. . mrrx . 

wry 

. (6) 

z = ZA mn cos (rct+e r ) Bin — Bin 

m, n 

b 

We write this 



. tmtx . ntry 
z = Z4> mn sin — sin — , 

in, n 

• 

. (6) 

where (f> mn are the normal coordinates. 

The 

kinetic 

energy is 




• 

. (?) 


and this is easily shown to bo 


T = X 5 pab fmn 

m.n “ 


( 8 ) 


The potential energy may be calculated in a manner 
similar to § 14. Referring to fig. 7 wo see that in the 
displacement to the bent position, the two tensions Toy 
havo done work T By . (arc AB—hx). As in §14, this 


reduces to approximately - T 8xSy. The other two 
tensions TSa: have done work jj 


The total 


potential energy is therefore 


( 9 ) 


In the case of the rectangular membrane this reduces to 

V = Z \ pabcVti. . ■ ■ ( 10 ) 

m,n° 

It will be seen that T and V are both expressed in the form 
of Chapter II, equation (38), typical of normal coordinates 
in mechanical problems. 
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§ 30. With a circular membrane such as a drum of 
radius a, we have to use plane polar coordinates r,9 
instead of Cartesians, and the solution of equation (1), 
apart from an arbitrary amplitude, is given in § 8, 
equation (35a). It is 

C08 

z = J m (nr) s;n md cos net. 


We have neglected the Y m {nr) term since this is not finite 
at r = 0. If we choose the origin of 9 properly, this normal 
mode may be written 

2 = J m [nr) cos md cos net. . . (11) 


If 2 is to be single-valued, m must be a positive integer. 
The boundary condition at r = a is that for all values of 
9 and /, .7 m (?w) cos m9 cos net equals zero, i e., J m (na) = 0. 
For any assigned value of m this equation has an infinite 
number of real roots, each one of which determines a 
corresponding value of n. These roots may bo found 
from tables of Bessel functions. If we call them n m , „ 
n m. a> ••• n m. • • •> then the frequency of (11) is nc/ 2n, 
i o. cn m and we may call it the (m, k) mode. The 

allowed values of m are 0, 1 , 2, ... and of k are 1 , 2, 3 

There are nodal lines which consist of circles and radii 
vectores. Fig. 10 shows a few of these modes of vibration, 
shaded parts being displaced in an opposite direction to 
unshaded. 


The nodal lines obtained in figs. 9 and 10 are known 
as Chladni’s figures. A full solution of a vibrating mem- 
brane is obtained by superposing any number of these 
normal modes, and if nodal lines exist at all, they will 
not usually be of the simplo patterns shown in these 
figures. As in the case of the rectangular membrane so 
also in the case of the circular membrane, the overtones 
bear no simple numerical relation to the fundamental 
frequency, and thus the sound of a drum is not very 
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musical. A vibrating bell, however, is of very similar 
type, but it can be shown * that some of the more important 
overtones bear a simple numerical relation to the funda- 
mental ; this would explain the pleasant sound of a well- 
constructed bell. But it is a little difficult to see why 
the ear so readily rejects some of the other overtones 



Fio. 10 

whose frequencies arc not simply related to the fundamental. 
A possible explanation f is that the mode of striking may 
bo in some degree unfavourable to these discordant 
frequencies. In any ease, wo can easily understand \\ ly 
a bell whose shape differs slightly from the conventional, 
will usually sound unpleasant. 

• Soe Slater and Frank. Introduction to Theoretical Physics, 

93 f Lamb, Dynamical Theory oj Sound (Arnold), 1010. p. 165. 
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§ 31. Examples 

(1) Find two normal modes which are degenerate (§ 28) 
for a rectangular membrane of sides 6 and 3. 

(2) Obtain expressions for the kinetic and potential 
energies of a vibrating circular membrane. Perform the in- 
tegrations over the ^-coordinate for the case of the normal mode 

z = A J m (nr) cos md cos net. 

(3) A rectangular drum is 10 cm. x20 cm. It is stretched 
to a tension of 5 kgm., and its mass is 20 gm. What is the 
fundamental frequency ? 

(4) A square membrane bounded by x = 0, a and y = 0, a 

is distorted into the shape z = A sin ^^.sin and then 

a a 

released. What is the resulting motion ? 

(5) A rectangular membrane of sides a and b is stretched 
unevenly so that the tension in the x direction is T x and in 
the y direction is Tj. Show that the equation of motion is 

T 8 2 z eh bh OL 

I I gji + 1 * - a — p -^ 5 - Show that this can bo brought into 

tho standard form by changing to new variables x/ \Zy” 

III nnd honco find tho normal modes. 

(6) Show that tho number of normal modes for the 
rectangular membrane of § 28 whoso frequency is less than N is 

approximately equal to tho area of a quadrant of tho ellipse 
x 3 y 2 4 p 

^ Jji = j Hence show that tho number is roughly 
rrpabN'/T. 

[Answers: 1. (2, 0) and (0, 1): in general (2m, n) and 
7,, > ; 2 - T = Jjrp/iVM 1 sin 2 

V = ^p;M* cos* [n*{ r )}* -f m*{ J m (;ir)}*// *] r dr, 
which becomes, after integration by parts 

V = l n P n * ci A* cos* net Ju (Jjnr)}* rdr; 3. 175-1 per seo. ; 

4. z =A sin (2:7 x/«) sin (3 ny/a) cos ( % /13ir ct/a ) ; 

5 ■ z = A sin (mnxju) sin [nzry/b) cos npt, 

p'p = m : T 1 /a ! +» ! T I /6 1 .] 
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LONGITUDINAL WAVES IN BARS AND SPRINGS 

$ 32. The vibrations which we have so far considered have 
all been transverse, so that the displacement has been 
perpendicular to the direction of wavo propagation. We 
must now consider longitudinal waves, in which the 
displacement is in the same direction as the wave. Sup- 
pose that AD (fig. 11) is a bar of uniform section and 



mass p per unit length. The passage of a longitudinal 
wave along the bar will be represented by the vibrations 
of each element along the rod, instead of perpendicular 
to it. Consider a small element PQ of length 8x, such 
that AP = x, and let us calculate the forces on this element, 
and hence its equation of motion, when it is dis- 
placed to a new position P'Q ' . If the displacement of 
P to P’ is £, then that of Q to Q' will bo so that 

P'Q' =hx-VH- We must first evaluate the tension at 
P\ We can do this by imagining 8x to shrink to zero. 

61 
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Then tlio infinitesimally small clement around P' will be 
in a stato of tension T where, by Hooke's Law, 


T/v = A . 


extension 
orig. length 


w . 8x+8£-8.r 
A Lim r 

Sx-0 bx 




Returning to the element P'Q', we see that its mass is the 

st a ' 

Therefore 


same as that of PQ, i e. p8x, and its acceleration is 



by (1). 


Thus the equation of motion for these longitudinal waves 
reduces to the usual equation of wave motion 

d ■ ' vhcre c ’ = Xlp ■ ■ (2) 


The velocity of waves along a rod is therefore \/(A Ip), 
a result similar in form to the velocity of transverse 
oscillations of a string. 

The full solution of (2) is soon found if we know the 
boundary conditions. 

(i) At a free end the tension must vanish, and thus, 


from (1), 


c {* 

--= 0, but the displacement will not, in general, 


vanish as well. 

(ii) At a fixed end the displacement $ must vanish, 
but. ti-.e tension will not, in general, vanish also. 
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§33. If wo are interested in the free vibrations of a 
bar of length l, wo shall use stationary type solutions of 
(2) as in § 8, equation (27). Thus 

£ = (a cos px+b sin px) cos {cpt +«}. 

If we take the origin at one end, then by (i) dg/dx lias to 
vanish at a; = 0 and x = l. This means that b = 0, and 
sin pi = 0. i.e. pi = nn, where n = 1, 2, ... . The free 
modes are therefore described by the functions 

71TTX In VCl | ,o\ 

g = a„ cos -j- cos j — + €„ j . • W 

This normal mode has frequency nc/21, so that the funda- 
mental frequency is e/2 1, and the harmonics are simply 
related to it. There are nodes in (3) at the points x — ll-n, 
3i/2n, 5l/2n, .... (2n — l)l/2n ; and there are antinodes (§ 6) 
at x = 0, 2I/2n, 4l/2n .... I . From (1) it follows that 
these positions are interchanged for the tension, nodes of 
motion being antinodes of tension and vice versa. We 
shall meet this phenomenon again in Chapter VI. 


§ 34 The case of a rod rigidly clamped at its two ends 
is similarly solved. The boundary conditions are now 
that g = 0 at x = 0, and at x = l The appropriate 

solution of (2) is thus 

g = a n sin 009 | ’ ‘ ^ 

This solution has the same form as that found in Chapter II, 
§ 19, for the transverse vibrations of a string. 


§35. We may introduce normal coordinates for these 
vibrations, just as in §§ 22 and 29. Taking, for 
example, the ease of § 34, we should write 


. £ , • nnx 

g = L<Pn -J- . 


n - 1 




where 
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The kinetic energy of the element PQ is |p3x . so 
that the total kinetic energy is 

•* . 

\p?dx = ... (6) 

0 n 

The potential energy stored up in P'Q' is approximately 
equal to one-half of the tension multiplied by the increase 

1 d£ 

in length : i.e. - A — . d£. Thus the total potential energy is 

M CX 


j: 


j> ' 


( 7 ) 


§ 36. The results of §§ 33, 34 for longitudinal vibrations 
of a bar need slight revision if the bar is initially in a state 
of tension. We shall discuss the vibrations of a bar of 
natural length l 0 stretched to a length l, so that its equili- 
brium tension T 0 is 

T - A l ~ to 


I. 


. . . . ( 8 ) 

Referring to fig. 1 1 , we see that Sx now represents the length 
of P'Q' when in the stretched, non-vibrating state. The 
completely unstretched length is then fore not 8x but 

-°8x, so that the tension at P' is not given by (1), but 


by the modified relation 

Sx+S£— y 8x 

T/v = A Lira — 

8^0 io Sj; 

= T °+ JYx' U3ing ( 8 ) ' * 

I he mass of PQ is p 0 (/ 0 /)Sx where p 0 refers to the unstretched 
bar, so the equation of motion is 

p 0 (W* ~ = T q - -T r =— Sx 

cl- * dx 

A l v-P n r 
— j~dx from (9). 
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We have again arrived at the standard equation of wave 
motion 


d 2 { _ \_d_H 

fa 2 ~ c 2 ct- ' 



( 10 ) 


It follows that c = (l// 0 )r 0 , where c Q is the velocity under 
no permanent tension. Appropriate solutions of (10) are 
soon seen to bo 


. . mrx 

£ = a n sin — cos 


7177 ct 
~ 


+ e L ,»=1,2 ... (ID 


The fundamental frequency is c/2/, wliich, from (10), can 
bo written c 0 /2/ 0 . Thus with a given bar, the frequency is 
independent of the amount of stretching. 

The normal mode (11) has nodes where x = 0, //«, 
2Z/n, . . . Z. A complete solution of (10) is obtained by 
superposition of separate solutions of type (11). 


§ 37. Wo shall conclude this chapter with a discussion 
of the vibrations of a spring suspended from its top end 
and carrying a load M at its bottom end. When wo 
neglect the mass of the spring it is easy to show that 
the lower mass M (fig. 12) executes Simple Harmonic 
Motion in a vertical lino. Let us, however, consider the 
possible vibrations when we allow for the mass in <>t the 
spring. Put m = pi, where p is the unstretched mass per 
unit length and Z is the unstretchcd length. Wc may 
consider the spring in three stages. In stage (a) we have 
the unstretchcd spring of length Z. The element II o 
length 8a: is at a distance x from the top point A. In stage 

(6) wo have the equilibrium position when the spring is 
stretched duo to its own weight and the load at the bottom. 
The element PP' is now displaced to QQ' . P § ** displaced 
a distance X downwards and P' a distance X -J-8X. Lastly, 
in stage (c) we suppose that the spring is vibrating am 
the element QQ' is displaced to RIP. The displacement s 
of Q and Q' from their equilibrium positions are £ and £ . 
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Tho new length RR' is therefore 8x+8X +8£. The mass of 
the element is tho same as the mass of PP', viz. p8x, and 
is of courso the same in all three stages. 

We are now in a position to determine the equation 
of motion of RR'. The forces acting on it are its weight 


ABC 
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downwards and the two tensions at R and R'. Tho 
tension T/ { may he found from Hooke’s Law, by assuming 
that 8x is made infinitesimally small. Then, as in § 36, 




extension 
or g. length 


A Lim + 

5x~0 5x 



( 12 ) 
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So tho equation of motion of RR' is 
d 2 £ 

pSx = resultant force downwards 

= ffpS* +Tjr ~T r , 

= gpte + ^ 8x. 


Dividing by fix and using (12), this becomes 
d 2 $ , A (d 2 X 8H\ 

ct* 9 + p\dx 2 + dx-)' 

This last equation must be satisfied by £ = 0, since this is 
merely the position of equilibrium (b). So 


„ , A d 2 X 

0 = g + ~pTx i 


By subtraction we discover once more the standard equation 
of wave motion 


= «, = A = A J . . as) 

dx 2 c 2 cl 1 p m 

This result is very similar to that of § 36. However, 
before wo can solve (13) we must discuss tho boundary 
conditions. There are two of these. Firstly, when x = 0, 
we must have £ = 0 for all t. Secondly, when x = l, 
(i.e. the position of the mass M) wo must satisfy the law 
of motion 


Mg - [TJx-j- 



Using (12), tlxis becomes 

= _ A [dX dj] 

[dt 2 \ z ., 9 M [dx + dx\ x .,' 
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As before, this equation must bo satisfied by | = 0, since 
this is just the equilibrium stage ( b ). Thus 




So, by subtraction we obtain the final form of the second 
boundary condition 



The appropriate solution of (13) is 

£ = a sin px cos {yci-fc}. . . (15) 

This gives £ = 0 when x = 0, and therefore satisfies the 
first boundary condition. It also satisfies the other 
boundary condition (14) if 

pi tan pi = m/M. . . . (16) 

By plotting the curves y = tan x, y = (m/J/)/x, wo see 
that there are solutions of (16) giving values of pi in the 
ranges 0 to tt/2, tt to 377/2, .... The solutions become 
progressively nearer to nv as n increases. 

We are generally interested in the fundamental, or 
lowest, frequency, since this represents the natural vibra- 
tions of M at the end of the spring. The harmonics 
represent standing waves in the spring itself, and may bo 
excited by gently stroking the spring downwards when in 
stage (b). If m/M is small, the lowest root of (16) is 
small ; writing pi = z, we may expand tan z and get 

z(z-fz 3 /3+...) = m/M. 

Approximately 

z 3 (l -f z*/3) = m/M. 

Wo may put z 2 in the term in brackets equal to the first 
order approximation z 2 = mJM, and then we find for the 
second order approximation 


o 


m AI 
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The period of the lowest frequency in (15) is 2irjpc, i.e., 
27 tI/cz. Using the fact that c 2 = A / m, this becomes 



If the mass of the spring m had been 


neglected we should havo obtained the result 27t\/(IMIX). 
It thus appears that the effect of the mass of the spring 
is equivalent, in a close approximation, to adding a mass 
one-third as great to the bottom of the spring. 


§ Examples 

(1) Find the velocity of longitudinal waves along a bar 
whose mass is 2*25 gms. per cm. und for which the modulus 
is 9-0 . 10 0 dynes. 

(2) Two semi-infinite bars are joined to form an infinite 
rod. Ihcir moduli aro A! and A* and the densities arc Pi 
and p 2 . Investigate the reflection coefficient (see § 16 ) and 
the phase change on reflection, when harmonic waves in the 
first medium meet the join of the bars. 

(3) Investigate tho normal modes of a bar rigidly fastened 
at one end and free to move longitudinally at tho other. 

(4) A uniform bar of length l is hanging freely from one 
end. Show that tho frequencies of tho norma) longitudinal 
vibrations aro (n + J) c/2 t, whore c is the velocity of longi- 
tudinal waves in tho bar. 

(5) Tho modulus of a spring is 7-2 . 10 s dynes. Its moss 
is 10 gms. and its unstrctchcd length is 12 cms. A mass 

gins, is hanging on tho lowest point, and the top point is 
tixed. Calculate to an accuracy of 1 por cont. tho periods of 
the lowest two vibrations. 

(0) Investigate tho vertical vibrations of a spring of un- 
Btretched length 21 and mass 2m, supported at its top end 
and carrying loads M at tho mid point and tho bottom. 


[Answers: 1. 2 km. per sec.; 2. Jt = ( 


V ' 1 rl' VV *2 Pi 


^V(^\pi)+\ / (\ 2 p 2 )/ * 
|(r+J)nC' , ) 

\ j Hr ( ; 5. 1-690 secs., 

i i 


3- f = A, sin cos 

0-252 « !C s. ; 0. Period = 2n/L where fc*-3fccot,d+ cot* „/ = 1. 
k = Minim.) 
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WAVES IN LIQUIDS* 

§39. In this chapter we shall discuss wave motion in 
liquids. We shall assume that the liquid is incompressible, 
with constant density p. This condition is very nearly 
satisfied by most liquids, and the case of a compressible 
Huid is dealt with in Chapter VI. We shall further assume 
that the motion is irrotational. This is equivalent to 
neglecting viscosity and assuming that all the motions 
have started from rest due to the influence of natural 
forces such as wind, gravity, or pressure of certain bound- 
aries. If the motion is irrotational, we may assume 
the existence of a velocity potential if) if we desire it. 
It will be convenient to summarise the formula which 
we shall need in this work. 

(i) If the vector u | with components («, v, w) J 
represents the velocity of any part of the fluid, then from 
the definition of (f> 

u = -V<£ = -grad <}>, . . (1) 

so that in particular u = v = —dS/dy, 

W = — C<f)j?Z . 

(ii) On a fixed boundary the velocity has no normal 

• Beloro reading this chapter tho student is advised to read 
Rutherford s l eclor Methods, Chapter VI, from which several 
results will be quoted. 

t Using Clarendon type for vectors. 

} Many writers use (UjU^u-) for tho velocity components. We 
shall find (u, v, tv) more convenient for our purposos. It is necessary, 
however, to distinguish u, which is a vector representing tho velocity 
and u, which is just the x component of the volocity. 

00 
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component, and hence if djdv denotes differentiation along 
the normal 

d<f>ldv = 0. . • • • (2) 

(iii) Since no liquid will be supposed to be created or 
annihilated, the equation of continuity must express the 
conservation of mass ; it is 


du dv dw 


(3) 


(4) 


Combining (1) and (3), we obtain Laplace’s equation 

o. . 

V ? dx i + dy 2 + dz 2 

(iv) If H(x, y, z, t) is any property of a particle of the 

dH 

fluid, such as its velocity, pressure or density, then — 

is the variation of if at a particular point in space, and 
DH 

— is the variation of U when we keep to the same particle 
n fi 

of fluid, jjj- is known as the total differential coefficient, 
and it can be shown * that 


l.O. 


DH 

Dt 

DH 

Dt 


dH dH , dH , dll 

= ^ u — 4- v b w — • 

dt dx dy dz 



(v) If the external forces acting on unit mass of liquid 
can bo represented by a vector F, then the equation of 
motion of the liquid may bo expressed in vector form 


Du 

Dt 



1 

-VP- 

P 


• Soo Rutherford. § 
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In Cartesian form this is 




with two similar equations for v and w. 

(vi) An important integral of the equations of motion 
can be found in cases where the external force F has a 
potential V, so that F = — yK. The integral in question 
is known as Bernoulli’s Equation : 

i + ltf+r-% = c, . . ( 7 ) 

P 6 Cl 


where C is an arbitrary function of the time. Now 
according to (1), addition of a function of t to <f> does 
not affect the velocity distribution given by ^ ; it is often 

dcf> , 

convenient, therefore, to absorb C into the term — and 

at 

(7) can then be written 

-+|u ! |F-r = const. . . (8) 

p 2 at 


A particular illustration of (8) which we shall require later 
occurs at the surface of water waves ; here the pressure 
must equal the atmospheric pressure and is hence constant. 
Thus at the surface of the waves (sometimes called the 
free surface) 

1 Pi 

5 u 2 -f V — — = constant. . . (9) 


§ 40. Wo may divide the types of wave motion in 
liquids into two groups ; the one group has been called 
tidal waves/ and arises when the wavelength of the 
oscillations is much greater than the depth of the liquid. 
Another name for these waves is long waves in shallow 
water. With waves of this type tho vertical acceleration 

• Lamb, Hydrodynamics, Chapter VIII. 
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of the liquid is neglected in comparison with the horizontal 
acceleration, and we shall be able to show that liquid 
originally in a vertical plane remains in a vertical plane 
throughout the vibrations ; thus each plane of liquid 
moves as a whole. The second group may be called 
surface waves, and in these the disturbance does not 
extend far below the surface. The vertical acceleration 
is no longer negligible and the wavelength is much less 
than the depth of the liquid. To this group belong most 
wind waves and surface tension waves. We shall consider 
the two types separately, though it will be recognised that 
Tidal Waves represent an approximation and the resu • 
for these waves may often bo obtained from the formulae 
of Surface Waves by introducing certain restrictions. 

TIDAL WAVES 

8 41 Wo shall deal with Tidal Waves first. Here wo 
assume that the vertical accelerations may bo neglected 
One important result follows immediately. If we draw 
the z axis vertically upwards (as wo shall continue to do 
throughout this chapter), then the equation of motion 
the z direction as given by (6), is 

Dw 1 d? 

Dt ~ 9 P dz ‘ 

Du j . . 

We are to neglect — and thus 

8 1 = -gp, i.e. p = -gp- + constant. 
dz 

Let us take our xrj plane in the undisturbed free surface, 
and write £(*, y, 0 for the elevation of the water ab 
the point (x, y, 0). Then, if the atmospheric pressure is 

£ wo milt Ll v = Vo when z = £• So tho equation 

for the pressure becomes 

p = P'+gpit-z)- ■ • • (,0) 
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We can put this value of p into the two equations of 
horizontal motion, and we obtain 


Du _ _ SC Dv 
Dt ~ 9 dx’ Dt 




The right-hand sides of these equations are independent 
of z, and we deduce therefore that in this typo of motion 
the horizontal acceleration is the same at all depths. 
Consequently, as we stated earlier without proof, on still 
water the velocity does not vary with the depth, and 
the liquid moves as a whole, in such a way that particles 
originally in a vertical plane, remain so, although this 
vertical plane may move as a whole. 


§42. Let us now apply the results of the last section 
to discuss tidal waves along a straight horizontal channel 
whose depth is constant, but whoso cross-section A varies 



from place to place. Wo shall suppose that the waves 
move in the x direction only (extension to two dimensions 
will come later). Consider the liquid in a small volume 
(fig. 13) bounded by the vertical planes x, x+dx at P 
and Q. The liquid in the vertical plane through P is all 
moving with the same horizontal velocity u(x) independent 
of the depth. We can suppose that A varies sufficiently 
slowly for us to neglect motion in the y direction. Wo 
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have two equations with which to obtain the details of 
the motion. The first is (11) and may be written 


8u Su iu 
— + « — T»r 
8t dx cz 



Since u is independent of z, — = 0. Further, since we 
shall suppose that the velocity of any element of fluid 
is small, we may neglect u £ which is of the second order, 
and rewrite this equation 

. . . ( 12 ) 

it 9 dx 

The second equation is the equation of continuity. Equa- 
tion (3) is not convenient for this problem, but a suitable 
equation can be found by considering the volume of liquid 
between the planes at P and Q, in fig. 13. Let 6(x) be the 
breadth of the water surfaco at P. Then the area of the 
plane P which is covered with water is [A +6£]p therefore 
the amount of liquid flowing into the volume per unit 
time is [{A+bQu] P . Similarly, the amount flowing out 
per unit time at Q is [(A +&£)«]< j- The difference between 
these is compensated by the rate at which the level is 
rising inside the volume, and thus 

+*£)«!« = bdr Tt 


Therefore 


-!-{L4 +&{)«} 


6 


H 

it' 


Since blu is of the second order of small quantities, wo may 
neglect this term and the equation of continuity becomes 
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Eliminating u between 


(12) and (13) gives us the equation 


b 


cH 
it 2 




In the case in which A is constant, this reduces to the 
standard form 


, 2 

dx 2 e-ct 2 ’ 


Agjb 



This is the familiar equation of wave motion in one 
dimension, and we deduce that waves travel with velocity 
y/(Agjb). If the cross-section of the channel is rectangular, 
so that A = bh,h being the depth, 

c = y/igh) .... (16) 

With an unlimited channel, there are no boundary 
conditions involving x, and to our degree of approximation 
waves with any profile will travel in either direction. 
With a limited channel, there will be boundary conditions. 
Thus, if the ends are vertical, u = 0 at each of them. 

Wo may apply this to a rectangular basin of length l, 
whoso two ends are at x = 0, l. Possible solutions of (15) 
are given in § 8, equation (27). They are 

£ = (a cos p.r+j3 sin px) cos (cpi+e). 

Then, using (13) and also the fact that *4 = bh, we find 


5u 

dx 


C ~ (a cos px+B sin px) sin (cpt+e). 
h 


u = - (a sin px—fi cos px) sin (cpt+e). 


The boundary conditions u = 0 at x = 0, /, imply that 
j6 — 0, and sin pi = 0. So 
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It will be noticed that nodes of u and £ do not occur at 
the same points. 

The vertical velocity may be found from the general 
form of the equation of continuity (3). Applied to our 
case, this is 

t + !? = o. 

dx dz 

Now u is independent of z and w = 0 on the bottom of 
the liquid where 2 = —h. Consequently, on integrating 

8u —irrafi 

w =- (z + h) 8-x = ~uT 

We may use this last equation to deduce under what 
conditions our original assumption that the vertical 
acceleration could be neglected, is valid, hor similarly to 

Die viv 

(12), the vertical acceleration — is effectively — • »- c * 


(z+A)cos y sinj-y +«r|- ( 19 > 


7T 2 r 1 c'-a r 


(z + k) CoS 


rvx [ rvet 

T 008 \T 


+ <r 


The maximum value of this is and may be com- 

pared with the maximum ho. izontal acceleration nrca, • 
The ratio of the two is rrr h/l, ie. 2*A/A since, from . ( <) 
A = 2 1/r. We have therefore confirmed the ccmd.t on 
which we stated as typical of these long waves, m th at 
the vertical acceleration may be neglected if the wavdeng 
is much greater than the depth of water. 

§43. We shall new remove the restriction i»P«l “ 
the last section to waves in one dimension. Ut ns use 
the same axes as before and consider the rate of fl 
liqui d into a vertical bounded by 

surface, ^FG// is the bottom of the liquid, ^ 
the moving surface at height £ (*, y) above ABCD. 
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rate of flow into the prism across the face PEHS is 
[u(h-\-t,)dy] x , and the rate of flow out across RQFO is 
[u(h-\-{,)dy] x+<tx . The net result from these two planes is 



B 
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a gain — — {u(h+£)}dxdy. Similarly, from the other two 

c .r 

vertical planes there is a gain — — {v(h-\-t,)}dzdy. The 

Cx J 

total gain is balanced by the rising of the level inside the 
prism, and thus 

— T^'Mh + i^drdy - ~{v(h+^)}dxdy = ^ . dxdy. 

As in § 42, we may neglect terms such as and r£ and 
write the above equation of continuity 

d(hu) c[hv ) _ cr, 

ox + cy ct * * 


( 20 ) 
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We have to combine this equation with the two equations 
of motion (11), which yield, after neglecting square terms 

in the velocities 


du _ — = -n — 

ft ~ ~ 9 Jx ’ dt ~ 9 dy 



Eliminating u and v gives us the standard equation 



If h is constant (tank of constant depth) this becomes 


m . (23) 

a.r* + cy- 0 - ct* 


This is the usual equation of wave motion in two dimensions 
and shows that the velocity is y/{gh). If we are concerned 
with waves in one dimension, so that £ is independent of 

y (as in § 42) we put 0 = ° and retrieve (15). 

We have therefore to solve the equation of wave 
motion subject to the boundary conditions 
(i) w = 0 at z = — h, 

Q = 0 at a boundary parallel to the y axis, and 


— = 0 at a boundary parallel to the x axis, 

/ V d 

(iii) C ^~ = 0 at any fixed boundary, where — denotes 

differentiation along the normal to the boundary. Tins 
latter condition, of which (ii) is a particular case, can be 
seen as follows. If Ix+my = 1 is the fixed boundary, 
then the component of the velocity perpendicular to tins 
line has to vanish. That is, lu+mv = 0. By different .at.ng 
partially with respect to t and using (21), the condition (in) 

is obtained. 
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§ 44. We shall apply these formulae to two cases ; first, 
a rectangular tank, and, second, a circular one, both of 
constant depth. 

Rectangular tank . — Let the sides be x = 0, a and y = 0, 6. 
Then a suitable solution of (23) satisfying all the boundary 
conditions (i) and (ii) would be 

C = A cos — — cos cos [rvci-\-e), . (24) 

a b 


where p = 0, 1, 2 ... , q = 0, 1, 2, ... , and r 2 = p 2 /a 2 -\-q 2 lb 2 . 


This solution closely resembles that fora vibrating membrane 
in Chapter III, §28, and the nodal lines are of the same 
general type. The student will recognise how closely the 
solution (24) resembles a “ choppy sea.” 

Circular tank . — If the centre of the tank is origin and 
its radius is a, then the boundary condition (iii) reduces to 

^ = 0 at r = a. Suitable solutions of (23) in polar 

coordinates have been given in Chapter I, equation (35a). 
We have 


£ = A cos md J m [nr) cos (ertf+e) . (25) 


We have rejected the }' m solution since it is infinite at 
r = 0, and we have chosen the zero of 9 so that there is 
no term in sin mO. This expression satisfies all the condi- 
tions except the boundary condition (iii) at r = a. This 
requires that J m '(na) = 0. For a given value of rn (which 
must be integral) this condition determines an infinite 
number of values of n, whose magnitudes may be found 
from tables of Bessel Functions. The nodal lines are 
concentric circles and radii from the origin, very similar 
to those in fig. 10 for a vibrating membrane. The period 
of this motion is 2rr cn. 


§ 45. It is possible to determine the actual paths of 
individual particles in many of these problems. Thus, 
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referring to the rectangular tank of § 42, the velocities 
u and w are given by (18) and (19). We see that 

w —irr(z+h) r-nx 

- = — cot — . 

u i 1 

This quantity is independent of the time and thus any 
particle of the liquid executes simple harmonic motion 
along a line whose slope is given by the above value of 
w/u. For particles at a fixed depth, this direction changes 
from purely horizontal beneath the nodes to purely vertical 
beneath the antinodes. 

§ 46. We shall conclude our discussion of tidal waves 
by applying the method of reduction to a steady wave, 
already described in § 25, to the case of waves in a channel 
of constant cross-section A and breadth of water-line b. 
This is the problem of § 42 with A constant. U‘t c be the 
velocity of propagation of a wave profile. Then super- 
impose a velocity -c on the whole system, so that the 
wave profile becomes stationary and the liquid flows under 
it with mean velocity c. The actual velocity at any point 
will differ from c since the cross-sectional area of the liquid 
is not constant. This area is A +&£, and varies with £. 
Let the velocity be c+6 at sections where the elevation 
is £. Since no liquid is piling up, the volume of liquid 
crossing any plane perpendicular to the direction of flow 
ia constant, i.e. 

(A +6£)(c+0) = constant = Ac. . (26) 

We have still to use the fact that the pressure at the free 
surface is always atmospheric. In Bernouilli's equation 
at the free surface (9) we may put dftdl = 0 since the 
motion is now steady motion ; also V = <j£, at the free 
surface. So, neglecting squares of the vertical velocity, 

this gives 

i(c+0)Hi/S = 
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E limin ating 6 between this equation and (26), we have 

A~c 2 


U+ 60 * 


+ 2 gC = e 2 , 


i.e. 


2^ = * 


1- 


,4 2 


M +&£) 


5 } = *«*£ 


2 ^+ 6 $ 

U +«)*]■ 


Whence 


2?U+60 2 

6 2 , 1 + 6 £ 



If £ is small, so that we may neglect £ compared with A/b, 
then this equation gives the same result as (16), viz. 
c 2 = gA/b. We can, however, deduce more than this 
simple result. For if £>0, the right-hand side of (27) is 
greater than gA/b, and if £<0, it is less than gA/b. Thus 
an elevation travels slightly faster than a depression and 
so it is impossible for a long wave to be propagated 
without change of shape. Further, since the tops of waves 
travel faster than the troughs, we have an explanation of 
why waves break on the sea-shore when thoy reach shallow 
water. 


SURFACE WAVE8 

§ 47 . We now consider Surface Waves, in which tho 
restriction is removed that the wavelength is much greater 
than the depth. In these waves the disturbance is only 
appreciable over a finite depth of tho liquid. We shall 
solve tliis problem by means of the velocity potential cf>. 

must satisfy Laplace’s equation (4) and at any fixed 
boundary dtftldv = 0, by (2). There are, however, two 
other conditions imposed on <f> at the free surface. Tho 
first arises from I» rmulli’s equation (9). If the velocity 
is so small that u 2 may be neglected, and if tho only 
external forces are the external pressure and gravity, we 
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may put u« = 0 and V = gi to this equation, which 
becomes 

( 28 ) 


4 Act 


free surface 

The second condition can be seen as follows. A particle 
of fluid originally on the free surface will remamjo ahcajs 
Now the equation of the free surface, where z - 4(*. V. l > 
may be written 

0 =/(x, y, Z, o = £(•*. y> l )~ z ■ 

Consequently, / is a function which is always zero /or a 
partick on (ho free eurface. We may therefore use (5) 
with H put equal to /, and we find 

0 _S+.2 + ,?*-» 

u Dt 8t ax cy 

M f mot e £ _ 1 1 ( 8 i\ = - - — on the surface. 

Now from (28) - - - dl g & 


Thus - C is a small quantity of order of magnitude not 

8X dl ?£ 

greater than u ; consequently u - and v-, being of 

order of magnitude not greater than u*. may be neglected. 
We are left with the new boundary condition 

= w = . • • < 29) 
’ct az 

Combining (28) and (29) we obtain an alternative relation 

= 0 . • • ( 3 °) 

ct 2 ™ dz 

We summarise the conditions satisfied by as follows . 

(i) Laplace’s equation vV = 0 m thc ,iqU1<1 ' jjj 

(ii) dtf>ldv = 0 on a fixed boundary ..-(-) 
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(iii) £ = - on the free surface 

0 cl 

ei ?<£ 

(iv) — = — — - on the free surface 

cl cz 


( 28 ) 


(29) 


(v) + g = 0 on the free surface . . (30) 

8t“ 8z 


Only two of the last three conditions are independent. 


§ 48. Let us apply these equations to the case of a 
liquid of depth h in an infinitely long rectangular tank, 
supposing that the motion takes place along the length 
of the tank, which we take as the x direction. The axes 
of x and y lie, as usual, in the undisturbed free surface. 
Condition (i) above gives an equation which may be 
solved by the method of separation of variables (see § 7), 
and if wo want our solution to represent a progressive 
wave with velocity c, a suitable form of the solution would 
be 

<f> = (.4e m *4- Be~ mz ) cos m{x— ct). 

A, B,m and c are to be determined from the other condi- 
tions (ii)-(v). At the bottom of the tank (ii) gives 8<f>ldz= 0, 
i.e. Ae~ mh —Be, mK = 0. So Ae~ mh = Bn™* = \C , say, and 
hence 

4> = C cosh 7 »(z+A) cos tn(x-cl). . (31) 

Condition (v) applies at the free surface whero, if the 
disturbance is not too large, we may put z = 0 ; after 
some reduction it becomes 


c * = (< 7 >) tanh mh. 


Since m = -n,\, where X is the wavelength, we can write 
this 



( 32 ) 
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Condition (iii) gives us the appropriate form of £ ; it is 

__ rn j l 8 j n m(x — ct). 

9 

This expression becomes more convenient if wo write a 
for the amplitude of £ ; i.e., a = ^ cosh mh. Then 

. . . . (33) 


£ = a sin m(x— ct), 
ga cos 


(34) 


me cosh m/i 


If the water is very deep so that tanli (2 tt«/A) 1, then 

(32) becomes c> = jA/2ir, and if it is very shallow so 
that tanh (2nh/X) = 2tt^/A, we retrieve the formula of § 

for long waves in shallow water, viz. c 2 — >jh. 

We have seen in Chapter I that stationary waves result 
from superposition of two opposite progressive harmonic 
waves. Thus we could have stationary waves analogous 
to (33) and (34) defined by 

/og\ 

r = a sin mx cos met, • 1 ' 




ga cosh m(zj h) ^ sin mcl 


(30) 


me cosh mh 


We could uso these last two equations to discuss stationary 
waves in a rectangular tank of finite length. 

§ 49. Wo shall now discuss surface waves in two dimen- 
sions, considering two cases in particular. 

rtedamukr Lk.- With a rectangular tank bounded 
by the planes r=0 ,..nd,= 0, b it is eas.ly verified 
that all the conditions of § 47 are satisfied by 

£ = A cos — cos cos rcl , 

gA « ah r(H h ) ^ pwx c(js ?rry ^ ^ _ 

^ rc cosh rh <* & 
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where 


p = 1 , 2 , ; q = 1 , 2 , .... ; r ! = r 5 (p 2 /a 2 +? 2 / 6 2 ) and 

c 2 = ( 0 /r) tanh rh. ... (37) 

Circular tank . — Suppose that the tank is of radius a 
and depth h. Then choosing the centre as origin and 
using cylindrical polar coordinates r, 9, z, Laplace’s equation 
(cf. Chapter I, § 7) becomes 


H , ie l, i &<f> n 

dr 2 ^ r dr ^ r 2 cd 2 ' Sz* ~ 


(38) 


A suitable solution can be found from Chapter I, equation 
(35a), which gives us a solution of the similar equation 

4. 1 + 1 'ti _ 1 _ o 

dr 2 t rcr + -* U 


r 2 c0 2 


c 2 ci- 


in the form 



, . COS COS 

r. < nr) sin ml) sin 


Tn this equation let us make a change of variablo, writing 
ct = iz, where i 2 = —1. We then get Laplace’s equation 
(3S) and its solutions are therefore 



m 


C08 


cosh 


r. (nr > sin mB ainh nz ' m = °' >- 2 


In our problem we must discard the Y solution as Yjr) 
is infinite when r = 0. So, choosing our zero of 6 suitably, 
wo can write 

<}> = Jm(nr) cos mO (.4 cosh nz+B sink nz). 

At the bottom of the tank condition (ii) gives, as in §48, 
A sinh nh = B cosh » h. so that 

<f> = C J m (nr) cos md cosh n(z+h). 

I he constants m and n are not independent, since we 
have to satisfy the boundary condition at r = a. This 
gives J m '(na) = 0, so that for any selected m, n is restricted 
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to have one of a certain set of values, determined from 
the roots of the above equation. The function C above 
will involve the time, and in fact if we are interested in 
waves whose frequency is /, we shall try C cc sin 2nft. 
Putting C = D sin 2i rft, where D is now a constant 
independent of r, d, z or t, we have 

(ft = DJ m (nr) cos md cosh n(z+h) sin 2v ft. (39) 

The boundary condition § 47 (iii) now enables us to find £ ; 
it is 

_ ZttDJ cQgjnfl cosh nh cos 277 ft . (40) 

Q 

The remaining boundary condition § 47 (iv) gives us the 
period equation ; it is 

— 4t7 2 / 2 D J m (nr) cos md cosh nh sin 2nft 

+gnD J m (nr) cos md sinh nh sin 2-nft = 0. 

i_ 0i 4 tt 2 / 2 = gn tanh nh. . . . (41) 

For waves with a selected value of m (which must be 
integral) n is found and hence, from (41) / is found. We 
conclude that only certain frequencies aro allowed. Apart 
from an arbitrary multiplying constant, the nature of the 
waves is now completely determined. 

§ 50. In § 48 we discussed the progressive wavo motion 
in an infinite straight channel. It is possible to determine 
from (34) the actual paths of the particles of fluid in this 
motion. For if X, Z denote the displacements of a particle 
whose mean position is (.r, z) we have 

cVi ga cosh m(z + A) . , ,, 

X = — — = — sin m{x — cl), 

'Ox c cosh mh 

2 = _ = _ ?? ■i'd' ’’■(*+ *) me 

dz c cosh mh 
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in which we have neglected terms of the second order of 
small quantities. Thus 

ga cosh m(z+h) 

X = 5 r — 7 cos m(x—ct), 

me- cosh mh 

ga 8 inh m(z-\-k) . 

Z = 5 i — r - sin m[x—cl). 

me- cosh mh 

Eliminating t, we find for the required path 

X - Z 2 o 2 a 2 

-L — y /42) 

cosh 2 m(z+A) Binli 2 m(;-f-A) m’-c* cosh- mh' 

These paths are ellipses in a vertical plane with a constant 
distance (2 ga/mc-) sech mh between their foci. A similar 
discussion could be given for the other types of wave 
motion which wo have solved in other paragraphs. 

§ 51. The Kinetic and Potential energies of these waves 
aro easily determined. Thus, if we measure the P.E. 
relative to the undisturbed state, then, since £(*, y) is the 
elevation, the mass of liquid standing above a base dA 
in the xy plane is p£ dA. Its centre of mass is at a height 
K. antl thus the total P.E. is 

| \gppdA ( 43 ) 

the integral being taken over the undisturbed area of 
surface. Likewise the K.E. of a small element is t pu 2 dr, 
dr being the element of volume of the liquid, so that the 
total K.E. is 



the integral being taken over the whole liquid, which may, 
within our approximation, be taken to be the undisturbed 
volume. 

With the progressive waves of § 48, £ and <f> are given 
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by (33) and (34), and a simple integration shows that the 
K.E. and P.E. in one wavelength (27r/m) are equal, and 
per unit width of stream, have the value 

\gpa 2 A (45) 


In evaluating (44) it is often convenient to use Green's 
Theorem in the form * 



^ dS. 

cv 


The latter integral is taken over the surfaco 5 which 
bounds the original volume, and d/d v represents differen- 
tiation along the outward normal to this volume. Since 
d<f>/dv = 0 on a fixed boundary, some of the contributions 
to T will generally vanish. Also, on the free surface, if 
£ is small, we may put d<f>l8z instead of 8<f>ldv. 


§ 52. We shall next calculate the rate at which energy 
is transmitted in one of these surface waves. We can 



illustrate the method by considering the problem discussed 
in §48 i.e. progressive waves in a rectangular tank ol 
depth h. Let AA' (fig. 15) be an imaginary plane fixed 
in the liquid perpendicular to the direction of wave 
• Soo Rulhorford, Chapter VI, p. 6G (ii). 
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propagation. We shall calculate the rate at which the 
liquid on the left of A A' is doing work upon the liquid on 
the right. This will represent the rate at which the energy 
is being transmitted. Suppose that the tank is of unit 
width and consider that part of AA' which lies between 
the two lines z, z+dz (shown as PQ in the figure). At 
all points of this area the pressure is p, and the velocity 
is u. The rate at which work is being done is therefore 

f° 

pudz. Thus the total rate is I pudz. W’e use Bernoulli’s 

equation (8) to give us p ; since u 3 may be neglected, 
and V = gz, therefore 


e<t> 

v = Po+p— -gpz. 


Now, according to (1) u = -d^jdx and 


irum 


, gi cosh m(z+h) 

9 — ~ ; — i — cos m(x— cl), 

me cosh mh 


Putting these various values in the required integral we 
obtain 


sin m(.r— cl) 



o 

-/> 


g'i coshmfc-l h) 
c cosh mh 


(Po-9p')d: 


i • < f° cg"-a- cosh 

4- sin* m(x — ct) _ 

J -h c a. 


2 m(2+A) 

slrmA 


dz. 


This expression fluctuates with the time, and we are 
concerned with its mean value. The moan value of 
sm m(x-cl) is z.-ro, and of sin 2 rnfx-c/) is $. Thus the 
mean rate at which work is being done is 

pg-a- n 1*0 

-ST & c\\ 2 mh cosh 2 m(z-\-h)dz. 

~ c J -h 

After some reduction this becomes 


\gpa 2 c (I + -2mh coscch 2 mh). 
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Tn terms of tho wavelength A = 2 t r/m, this is 

1 f 4rrA . 4 t:A| 

- gpn 2 c|l+ — cosech — j . . (46) 

Now from (45) we see that the total energy with a stream 
of unit width is \qpa~ per unit length, thus the velocity 
of energy flow is 

c ( . 4-Trh , 4 -h I 

— cosech — V . . (47) 

2\ A A j 

We shall see in a later chapter that this velocity is an 
important quantity known as the Group V elocity. 


§53. In the preceding paragraphs we have assumed 
that surface tension could bo neglected. However, with 
short waves this is not satisfactory and wo must now 
investigate the effect of allowing for it. When we say 
that the surface tension is T, wo mean that if a line of 
unit length is drawn in the surface of tho liquid, then 
tho liquid on one side of tills lino exerts a pull on tho 
liquid on tho other side, of magnitude T. Thus the effect 
of Surface Tension is similar to that of a membrane 
everywhere stretched to a tension T (as hi Chapter IN, 
§27), placed on tho surface of the liquid. We showed in 
Chapter III that when the membrane was bent there was 
a downward force per unit area approximately equal to 

_t 4 . Thus in fig. 16, the pressure just inside 

\dx 2 dy 2 } 

tho liquid does not equal tho atmospheric pressure p 0 , 
but rather 



Tho reader who is familiar with hydrostatics will 
recognise that tho excess pressure inside a stretched film 
(as in a soap bubble) is Kl/^ + l/i?*). where R x and IL 
are the radii of curvature in any pair of perpem uuar 


F 
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planes through the normal to the surface. We may put 
= —c'-^JBxr and R 2 = —B 2 t,!By 2 to the first order of 
small quantities, and then (48) follows immediately. 



Thus, instead of being p = p 0 at the free surface of 

the liquid, the correct condition is that p-f- t(°— 4- — \ 

} dx 2 By 2 ) 

is constant and equal to p 0 . We may combine this with 
Bernoulli's equation (9), in which we neglect u 2 and put 
V = gz. Then the new boundary condition which replaces 
§ 47 (iii) is 

d(f> r . t [c 2 { , a*n 

dt p\Bx 2 + By 2 ] ~ 0 • (49 > 

We still have the boundary condition § 47 (iv) holding, 
sinco this is not affected by any sudden change in pressure 
at the surface. By combining (29) and (49) we find the 
new condition that replaces § 47 (v). It is 

n <4> T/0* ,d 2 \e<f> 

T? + g - p \e*+ = 0 • < 80 ' 

We may collect these formula- together ; thus, with surface 
tension 


(i) y 2 <£ = 0 in the body of the liquid 

(ii) 8<f> Bv = 0 on all fixed boundaries 

it 


(mi 3 - rf + IfS . 


P [ ?x ' 


0H 


(4) 

( 2 ) 


on the free surface 

(49) 

(iv) d^dt = -e<f>ldz on the free surface . . (29) 

tv\ ^ Tf a 2 an a<* 
n sr- +g ez + ep \ £ = 0 °" the fr “ 

J) surface . (50) 

Only two of the last three equations are independent. 



WAVES IN LIQUIDS 


83 


§54. Waves of the kind in which surface tension is 
important are known as capillary waves. We shall 
discuss one case which will illustrate the conditions (i)-(v). 
Let us consider progressive type waves on an unlimited 
sheet of water of depth h, assuming that the motion takes 
place exclusively in the direction of x. Ihen, by analogy 
with (31) wo shall try 

<f> = C cosh m(z+h) cos m(x—cl). • • (51) 

This satisfies (i) and (ii) (iv) gives the form of £, which is 
£ = (C/c) sinh ink sin m(z- cl). . • (52) 

We have only one more condition to satisfy ; if wo choose 
(v) this gives 

-fflV-C cosh mh cos m{x-cl)+mCg sinh mh cos m(x-ct) 


-f - m 3 C sinh mh cos m(x—cl) — 0, 

P 

c* = (j/m+Tm/p) tanh mh. . . . (53) 


This equation is really the modified version of (32) when 
allowance is made for the surface tension , it T — 0, it 

reduces to (32). . 

When h is large, tanh mh = 1, and if we write in = -tt/a, 

wo have 


g X 2t7T 

c'- = *- + -y~- 
2 TT Xp 



The curve of c against A is shown in fig. 17, from which 
it can be seen that there is a minimum velocity which 
occurs when A* = 4w*T fgp. Waves shorter than tins, in 
which surface tension is dominant, are called ripples, an 
it is seen that for any velocity greater than the minimum 
there are two possible types of progressive wave, oneo 
which is a ripple. The minimum velocity is (igJlp) . 
and if, as in water, T = 75, p = 100 and g = 931 e.g.s 
units, this critical velocity is about 23 cms. per sec., an 
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the critical wavelength is about 17 cms. Curves of c against 
A for other values of the depth h are very similar to fig. 17. 



§ 55. Examples 

(1) Find the Potential and Kinetic energies for tidal 
waves in a tank of length l, using the notation of § 42. 

(2) Find the velocity of any particle of liquid in the 
problem of tidal waves in a circular tank of radius a (§ 44). 
Show that when m = 0 in (25), particles originall)’ on a vertical 
cylinder of radius r coaxial with the tank, remain on a coaxial 
cylinder whose radius fluctuates ; find an expression for the 
ainplitudo of oscillation of this radius in terms of r. 

(3) Tidal waves are occurring in a square tank of depth h 
and side a. Find the normal modes, and calculate the Kinetic 
and Potential energies for each of them. Show that when 
more than one such mode is present, the total energy is just 
the sum of the separate energies of each normal mode. 

(4) What aro the paths of the particles of the fluid in the 
preceding question ? 

(5) A channel of unit width is of depth h, where h = hx, 
k being a constant. Show that tidal waves are possible with 
frequency p/2ir, for wh 

C — -1.7o(axb2) C03 pi' 
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where a* = 4 p'/kg, and J 0 is Bessel’s function of order zero. 
It is known that the distance between successive zeros o 
JJx) tends to tt when i is largo. Hence show that the wave- 
length of these stationary waves increases with increasing 
values of x (This is the problem of a shelving beach.) 

(6) At the end of a shallow tank, we have * = 0, and the 

depth of water h is h = V*"- Also tho br( ' a,lth J f . tho ta ‘* 
b is given by b = b„x\ Show that tidal waves of frequency 

pl2n are possible, for which 

( = Ax u J g (rx*) cos pt 9 

where . 

« - l-m. n» - P'I9K r = a/s. 2« = l-2m-» and q = | «/- |. 

Use the fact that JJx) satisfies tho equation 

dx* xdx \ x ) 

(7) Prove directly from the conditions (i)*(v) in § 47 
without using the results of § 48 that tho velocity of surface 
waves in a rectangular channel of infinite depth is 

(8) Find the patlis of particles of fluid in the case of surface 
waves on an infinitely deep circular tank of radius a. 

(9) A tank of depth h is in tho form of a sector of a circle 
of radius a and angle 72°. What are tho allowed normal modes 

for surface waves ? , . , „r 

(10) If X, Y, Z denotes tho displacement of aportic 

fluid from its mean position x, y, z in a rectangu ar “ 1 ^ 

sides a and b when surface waves given by equa'io ) 
occurring, prove that tho path of the particle is the straight lino 

JL cot — X = — cot — Y = -coth r(z + /i) Z. 

pir a qrr b r 

(11) Show that in surface waves on a cylindrical tank of 
radius a and depth h, tho energies given by the norm 

(39) are 

v = 2 * D -£l coshViA cos*2 nft f° J m Hnr) r dr, and 

g Jo 

T - nnpD 1 sin*2ir ft cosh nh sinh nh J r dr ‘ 

2 Jo 
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Lso tho fact that the total energy must be independent of the 
time to deduco from this that the period equation is 

4t r 2 /* = gn tanh nh. 

(12) Show that when we use cylindrical polar coordinates 
to describe the capillary waves of § 53, the pressure condition 
at the free surface § 53 (iii) is 

cl JS> p\dr 3+ rdr ' r* cv* | ' 

Use tin’s result to show that waves of this nature on a 
eircu'ar basin of infinite depth are described by 

<t> = C J„(nr) cos m6 e"* cos 2i rft, 

-nC 


i = 


-V 


J m (nr) cos tn$ sin 2nft, 


where J m ’(na) = 0 and 4, r 5 /* = gn+Tn'fp. 

(13) Show that capillary waves on a rectangular basin of 
sides a, b and depth h are given by 

CO>h mnX tint/ 

* " •' 1 ~^T ni, rh - C ° S T 003 b CM 2 " A 

y r - 4 rnirx nnu . 

t ~ cos cos ~ sin 2 rft, 

u b 

where m = 0, 1. 2. ... ; n = 0, 1. 2 ... ; r* = „*(m*/a* + n*/&*), 
and the period equation is 

4 ~'P = {'jr+Ti^lp) tanh rh. 

\ t*rify, that when n = 0, this is equivalent to the result of 
§ 54, equation (53). 

[Answers : 

~ cl ^ 1 TirCt 


(1) i r Jphl, 2 COS 


* S i Op'or’- sin 2 + <,j ; 


(-) radial vel. is -(gA/c) cos md J n '(>ir) sin (cnt -f c), trans- 
verso velocity is (gAm/cnr) sin mO J m („r ) sin (cnt + c), 
(<jAjC)J t (nr) ; (3) £ = .1 cos (pr.xja) cos (oWu) cos (met la), 
ri = >>* + V 8 * K i sin* (rjrct/a), P.E. = * gp.U* 

(8) X : Y:Z = 


COS* (rrrd/a) ; (4) £ = P ta n 7 — cot 9 — ' 

1 <7 


a a 

" r,7m ion! / l "T”' Jm( " r) tnn ^m(nr); (9) Same as in 

cqns. (39) (41) except that m = 6fc/2. where i = 0, 1, 2. . . .] 



CHAPTER VI 


SOUND WAVES 

§ 56. Throughout Chapter V we assumed that the liquid 
was incompressible. An important class of problems is 
that of waves in a compressible fluid, such as a gas. In 
this chapter wo shall discuss such waves, of which sound 
wavos are particular examples. The passage of a sound 
wave through a gas is accompanied by oscillatory motion 
of particles of the gas in the direction of wave propagation. 
Theso waves arc therefore longitudinal. Since the density 
p is not constant, but varies with the pressure p, wo require 
to know the relation between p and p. If the compressions 
and rarefactions that compose the wave succeed each other 
so slowly that the temperature remains constant (an 
isothermal change) this relation is p = kp. But normally 
this is not the caso and no flow of heat, which would be 
needed to preserve the temperature constant, is possible ; 
in such cases (adiabatic changes) 


where k and y are constants depending on the particular 
gas used. We shall use (1) when it is required, rather 
than the isothermal relation 

§67. There are several problems in the propagation 
of sound waves that can be solved without using the 
apparatus of velocity potential <£ in the form in w u< > 
wo used it in Chapter V, §§ 47-54 ; wo shall therefore 
discuss some of theso before giving the general dove opim n 
of the subject. 
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Oar first problem is that of waves along a uniform 
straight tube, or pipe, and we shall be able to solve this 
problem in a manner closely akin to that of Chapter IV, 
§32, where we discussed t lie longitudinal vibrations of a 
rod. Wo can suppose that the motion of the gas particles 
is entirely in the direction of the tube, and that the velocity 
and displacement are the same for all points of the same 
cross-section. 

Suppose for convenience that the tube is of unit cross- 
sectional area, and let us consider the motion of that 
part of the gas originally confined between parallel planes 
at P and Q a distance dx apart (fig. 18). The plane P 


p Q 


X 

1 1 



\\ 


x.t 

1 I 

I 4 <*v| 



P' Q‘ 
Fio. 18 


is distant x from some fixed origin in the tube. During 
the vibration let PQ move to P'Q', in which P is displaced 
a distance £ from its mean position, and Q a distance 
The length P'Q' is therefore rf.r+rf£. Wo shall 
find the equation of motion of the gas at P'Q'. For this 
purposo wo shall require to know its mass and the 
pressure at its two ends. Its mass is the same as the 
mass of the undisturbed element PQ, viz. P( /ix, where p 0 
is the normal average density. To get the pressure at P' 
we imagine the element dx to shrink to zero ; this gives 
the local density p, from which, by (1), we calculate the 
pressure. We have 
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if we may neglect powers of d£/dx liigher than the first. 
The quantity (p—po)lpo will often occur in this chapter ; 
it is called the condensation s. Thus 

s = — 8£jdx , p = p 0 (l+ 5 )- • • I**) 


The net force acting on the element P'Q' is Pp.—P Q .< an< ! 
hence the equation of motion is 


3 2 £ cp 

=Pr~ p Q- = ~dx ‘ 





We may rewrite (4) in the form 

3 2 £ dp dp dp 3 2 £ 
p0 W ~ dpdx~ po dpdx 2 


from (2). 


It appears then that £ satisfies the familiar equation of 
wave motion 


8 *$ 

8x 2 




This equation shows that waves of any shape will bo 
transmitted in either direction with velocity \Z(dp/dp). 
In the case of ordinary air at 0° C., using (1) as the relation 
between p and p, wo find that the velocity is c = 3.1- 
metres per sec., which agrees with experiment. Newton, 
who made this calculation originally, took the isothermal 
relation between p and p and, naturally, obtained an 

incorrect value for the velocity of sound. 

A more accurate calculation of the equation of motion 
can bo made, in which powers of 8^dx are not neglected, 
as follows. From (2) wo have the accurate result 
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So, now using (4) in which no approximations have been 
made, 



cty+'dx*’ 


yW 


M9 


*n = m i n 

dt* p 0 {i +c$dx}Y+idz 2 * * 


Equation (5) is found from (6) by neglecting d£/dx compared 
with unity. A complete solution of (6) is, however, beyond 
the scope of this book. It is easy to see that, since (6) 
is not in the standard form of a wave equation, the velocity 
of transmission depends upon the frequency, and hence 
that a wave is not, in general, transmitted without change 
of shape. 


§58. We must now discuss the boundary conditions. 
With an infinite tube, of course, there are no such condi- 
tions, but with a tube rigidly closed at x — x 0 , wo must 
have £ = 0 at x = x 0 , since at a fixed boundary the gas 
particles cannot move. 

Another common typo of boundary condition occurs 
when a tube has one or more ends open to the ntmosphore. 
If wo suppose that the waves inside the tube do 
not extend their influence to the air beyond the end 
of the tube, then at all open ends the pressure must 
have the normal atmospheric value, and thus, from 
(1) and (2), dg;dx=0. Since the waves do extend a 
little outside the tube, this last equation is not strictly 
accurate. The usual modification is to increase the 
effective length of the tube by a small end-correction 
depending on the area of the cross-section of the tube. 
We shall not, however, include such corrections in 
this book. 

To summarise : 



1 d 2 £ . 

7i 77i ' n l * ie tu bo, and c 2 = dpldp . 

C Ci 
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(ii) £ = 0 at a closed end (") 

(iii) r = $ = 0 at an open end. ... (ty 

v dx 

§ 59. We shall apply these equations to find the normal 
modes of vibration of gas in a tube of length l. These 
waves will naturally be of stationary type. 

(a) Closed at both ends x = 0, l — This problem is the 
same mathematically, as the transverse vibrations of a 
string of length l, fixed at its ends (cf. Chapter II, § 19). 
Conditions (i) and (ii) of § 58 give for the normal modes 


f = i4 r sin cos -j- +€ r j , r — 1,2, .... (9) 

(b) Closed atx = 0, open atx = l(&“ stopped tube ”).— 
Here conditions (ii) and (iii) give f = 0 at x = 0, and 

— = 0 at x = l. The normal modes are 
dx 

f = ^ r s:n|r+^yCO8|(r+^)y+«,}. r = 0 . 1 ' 2 -"- ( ,0 > 

(c) Open at both ends x = 0, l- We have to satisfy 
the boundary condition (iii) d^jdx = 0 at x = 0, l. 

the normal modes are 


rvx [ met ) . o 

£ = i4 r cos -j -008 I ~y + f rj< r — t, S, 


( 11 ) 


In each case the full solution would be the superposition 
of any number of terms of the appropriate type with 
different r. The fundamental periods in the three cases 
are 2//c, 4 l/c, and 2//c respectively. The harmonics bear 
a simple numerical relationship to the fundamental, which 
explains the pleasant sound of an organ pipe. 

§ 60. We shall now solve a more complicated problem. 
We are to find the normal modes of a tube of unit sectional 
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area, closed at one end by a rigid boundary and at the 
other by a mass M free to move along the tube. Let 



Fig. 19 


the fixed boundary be taken as x = 0, and the normal 
equilibrium position of the moveable mass be at x = l 
(fig. 19). Then we have to solve the standard equation 
of wave motion with the boundary conditions that when 
x = 0. (b) gives £ = 0, and that when x = l the excess 


pressuro inside, p—p 0 , must be responsible for the 
acceleration of the mass M. This implies that 


P—Po = ^ Jp when x = 1. 


The first condition is satisfied by the function 

f = A sin nx cos M+e) . . . (12) 

To satisfy the second condition, wo observe that 

P~Po = (dp!dp)(p—Po) = —c 2 p 0 8£ldx, from (3). 

So this condition becomes 

•§--< — i 

Using (12) this gives, after a little reduction, 

nl tan nl = lp 0 /M. 

Tho allowed values of n are the roots of this equation. 
1 here is an infinite number of them, and when M = 0, so 
that tho tube is effectively open to the air at one end, 
we obtain equation (10) ; when M = co, so that the tube 
is closed at each end, we obtain equation (9). 


§61. So far we have developed our solutions in terms 
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of i, the displacement of any particle of the gas from its 
mean position. It is possible, however, to use the method 
of the velocity potential <f>. Many of the conditions which 
(f> must satisfy are the same as in Chapter V, but a few 
of them arc changed to allow for the variation in density. 
It is convenient to gather these various formula; together 
first. 

(i) If the motion is irrotational, as we shall assume, 

u = —yj» (cf. Chapter V, equation (1) ) . (13) 

(ii) At any fixed boundary, d<f>/dv = 0 (cf. Chapter V, 

equation (2)) (14) 

(iii) The equation of Continuity (cf. Chapter V, equation 

(3)) is slightly altered, and it is * 

d £+ v .pu = 0, 

ie - t + L { ^ ) +k (pv) +^ (pw) = (i ' (15) 

(iv) The equations of motion are unchanged ; if F is 

the external force on unit mass, in vector form, 
they are 

— = F— - VP (cf. Chapter V, equation (6)) . (16) 
Dt v 

(v) In cases where the external forces have a potential 

V, we obtain Bernoulli’s equation (cf. Chapter V, 
equation (8)) 

f— + $u 2 +P— T 7 = const. . . (17) 

J P Ct 

in which we have absorbed an arbitrary function of the 
time into the term dift/dt (cf. Chapter V, equation (8)). 

§ 62. In sound waves we may neglect all external 
forces except such as occurs at boundaries, and thus wo 

• Rutherford, § 67. 
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may put V = 0 in (17). Also we may suppose that the 
velocities are small and neglect u 2 in this equation. With 
these approximations Bernoulli’s equation becomes 


I 


dp 8<f> 

= const. 

ct 


Wo can simplify the first term; for = , 

and if the variations in density are small, dp/tlp may be 
taken as constant, and equal to c 2 as in (5). Thus 
Np , r,i P 

J p - c J - = c * lo &P = ^{lo&fl+aRlog^o}. 

,C ’° j~ = c 2 s 4-const., if s is small. If we absorb this con- 
stant in <j>, then Bernoulli’s equation takes its final form 

c'a-dftdt = 0 . . . (18) 

Laplace a equation for <f> does not hold because of the 
changed equation of continuity. But if u, v, w and a 
are small, (15) can be written in a simpler form by the 
aid of (13) ; viz., 

Ptfis/dt - = 0 . 

This is effectively the same as 


£a ia 


(19) 


Now let us eliminate s between (18) and (19), 
we shall find the standard equation of wave motion 


V 2 <£ 


c 2 ct 3 


and 


( 20 ) 


This shows that c is indeed the velocity of wave propaga- 
tion, but beforo wo can use this technique for solving 
problems, wo must first obtain the boundary conditions 
for At a fixed boundary, by (ii) dftdv = 0. At an 
open end of a tube, the pressure must be atmospheric, 
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and hence 5 = 0. Thus, from (18), 

d<f>ldl = 0 . . . • • ( 21 ) 

This completes the development of the method of the 
velocity potential, and we can choose in any particular 
problem whether we solve by means of the displacement 
£ or the potential tf>. It is possible to pass from one to 
the other, since from (3) and (18) 

!!=_,= _4^r- • • < 2 -’> 

dx c* ct 

§ 63. We shall illustrate these equations by solving the 
problem of stationary waves in a tube of length l, closed 
at one end (* = 0) and open at the other (x = l). This 
is the problem already dealt with in §59 (6), and with 

d 2 <f> 1 d 2 <f> 

the same notation, we require a solution of — = - — 

subject to the conditions 

d<f>ldx = 0 at x = 0, 
c<J>ldt = 0 at x = l. 


It is easily seen that 

</> = a cos mx cos (emt -fe) 

satisfies all these conditions provided that cos ml = 0, 
i.e. ml = tt/2. 377 / 2 , (r+i>* So the normal modes 


are 

^=“.co S (r+!) jCO S ](r+i) j+4 

and from this expression all the other properties of these 
waves nmy easily be obtained. The student is advised 
to treat the problems of § 59 (a) and (c)i n a similar manner. 

§ 64. Our next application of the equations of § 62 
will be to problems where there is spherical symmetry 
about the origin. The fundamental equation of wave 
motion then becomes (see Chapter I, equation (23)) 

8 2 <f> 2 C(f> 1 

+ r c7 ~ c 2 cl 2 ’ 
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with solutions of progressive typo 

<f> = l -f(r~d)-\- l -g(r+ct). 

There are solutions of stationary type (see Chapter I, 
equation (37)) 

, , cos cos 

6 = r~ l . mr . cmt. 
sin sin 

If the gas is contained inside a fixed sphere of radius a 
then we must have (f> finite when r = 0, and df/dr = 0 
when r = a. This means that 

, A . 

<P = — sin mr cos (cmt+e), 
with the condition 

tan ma =ma . . . . (23) 

This period equation has an infinite number of roots which 

approximate to ma = (n + l/2)»r when n is largo. So for 

its higher frequencies the system behaves very like a 

uniform pipe of length a open at one end and closed at 
the other. 

This analysis would evidently equally well apply to 
describe waves in a conical pipe. 


§65. Wo shall now calculate the energy in a sound 
wave. The Kinetic energy is clearly JipuW, where 

d V is an element of volume. We may put p = Po without 

1033 of accuracy. In terms of the velocity potential this 
may be written 


:/ 0 (V<f>)-d V = - l Pv I l Po f ^ ^ dS . (24) 

“ - J cv 


The last expression follows from Green’s theorem just as 
in 1 iaptcr \ , § 51, and the surface integral is taken over 
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the boundary of the gas. There is also Potential energy 
because each small volume of gas is compressed or rarified, 
and work is stored up in the process. To calculate it, 
consider a small volume F 0 , which during the passage 
of a wave is changed to V v If Sj is the corresponding 
value of the condensation, then from (3), we have, to the 
first degree in 

Fj = P 0 (l-Sj) . . . (25) 

Further, suppose that during the process of compression, 

V and 8 are simultaneous intermediary values. Then we 
can write the work done in compressing the volume from 

fP. . 

F 0 to F, in the form — p dV. But, just as in (25), 

J P# 

V = F 0 (l — a), and hence 

dV = — F 0 ds. 

We may also write p = p Q +(dpldp)(p—p 0 ) 

= Po+ c W- 

Thus the potential energy may be written 

| ' (Po+cW) = PoVi+^VW 

= P 0 (F 0 — FJ+ fcVoJVi*- 

This is the contribution to tho P.E. which arises from the 
volume V 0 . Tho total P.E. may be found by integration. 
Tho first term will vanish in this process since it merely 
represents tho total change in volume of the gas, which 
we may suppose to be zero. We conclude, therefore, that 

the Potential Energy isj ^ c 2 p 0 s 2 dV .... (26) 

It can easily be shown that with a progressive wave 
the K.E. and P.E. are equal; this does not hold for 
stationary waves, for which their sum remains constant. 

O 
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§ 66. We conclude this chapter with a discussion of 
the propagation of waves along a pipe whose cross-sectional 
area A varies slowly along its length. Our discussion is 
similar in many respects to the analysis in § 57. 

Consider the pipe shown in fig. 20, and let us measure 
distances x along the central line. It will be approximately 
true to say that the velocity u is constant across any 
section perpendicular to the x axis. Suppose that the 
gas originally confined between the two places P, Q at 



distances x, x-f dx is displaced during the passage of a 
wave, to P Q', the displacement of P being £ and of Q 
being Consider the motion of a small prism of gas 
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such as that shaded in the figure ; its equation of motion 
may be found as in § 57, and it is 



We must therefore find the pressure in terms of This 
may be obtained from the equation of continuity, which 
expresses the fact that the mass of gas in P'Q' is the same 
as that in PQ. Thus, if p is the density, 


Po A(x) dx = p A[x+Z) ■ {dx+d£}, 

.8A 


i.e. 


Po A(x) = p]A(x)+£— 


1 + 


e J\ 

dx\’ 


Neglecting small quantities, this yields 

f af fat 

= l + Fx + ASi 


Therefore 


p = po { l - 8 £-ii£}- po \ l ~aL (M) ) 


(28) 


Eliminating p between (27) and (28) we find 



where, as usual, 




So the equation of motion is 



In the case in which A is constant this reduces to tbo former 
equation (5). An important example when A is not 
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constant is the so-called exponential horn used on the 
best acoustic gramophones ; here the tube is approximately 
symmetrical about its central line and the area varies with 
the distance according to the law A = A 0 e 2 **, where a and 
A 0 are constants. 

With this form of *4, (29) reduces to 



A solution is possible by the method of separation of 
variables (see § 7). We soon find 

£ = e icn, {B l e m ' z -' r B 2 e mtZ }, 


where m l and m 2 are given by — a±V( a2 ~ n 2 ). In most 
exponential horns n 2 is considerably larger than o 2 in the 
range of audible frequencies, so that and m a may be 
written —a±in. Thus 


£ = e - "{ B l e in{et + £ 2 e ,n(rt + *>} . 



Tho first term represents a wave going outwards and the 
second a wavo coming inwards. Wo conclude from this 
that waves can be sent outwards along tho horn with a 
velocity c which is approximately independent of the 
frequency, and with an attenuation factor e~* x which is 
also independent of the frequency. It is this double 
independence which allows good reproduction of whatever 
waves are generated at the narrow end of the horn, and 
which is responsible for this choice of shape in the best 
gramophones. Other forms of .4 will not, in general, give 
rise to the same behaviour. 


§ 67. Examples 

(1) Use tho method of § 58 to investigate sound waves in a 
closed rectangular box of sides o„ o, and a,. Show that if the 
box is large, the number of such waves for which the frequency 
is less than n is approximately equal to one-eighth of the 
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volume of the quadric x J /ai 1 +!/V a 2 2 + 2 V a 3 2 = 4 m*/c*. Hence 
show tliat this number Is approximately 4trn 3 ci 1 fi 2 a 3 /3c 3 . 

(2) Investigate the reflection and transmission of a train 
of harmonic waves in a uniform straight tube at a point 
where a smooth piston of mass M just fits into the tube and 
is free to move. 

(3) Show tliat the kinetic and potential energies of a plane 
progressive wave are equal. 

(4) Show that the kinetic and potential energies of 
stationary waves in a rectangular box have a constant sum. 

(5) Find an equation for the normal modes of a gas which 
is confined between two rigid concentric spheres of radii 
a and 6. 

(6) Show that a closer approximation to the roots of 
equation (23) Is ma = (n + J)tt — 1/{(h + J)rr}. 

(7) Find numerically the fundamental frequency of a 
conical pipe of radius 1 metre open at its wide end. 

(8) The cross-sectional area of a closed tube varies with the 
distance along its central line according to the law A = AqX n . 
Show that if its two ends are x = 0 , and x = l, then standing 
waves can exist in the tube for which the displacement is 
given by the formula 

( = x< l “ n ^ 2 J m ( 7 x/c) cos {qct + (}< 

where m = (n + l)/2 and J m (?// c ) = 

Use the fact that J m (x) satisfies the equat on 


dV 1 dJ 

dx 2 + x dx 



0 . 


[Answers: 2. reflection coefft. 11 = {l+4p 0 V- w *' 42 }“ 1 * 
transmission coefft. T = (l+MV/W)" 1 ! for dofmltlonS ° f 
R and T see § 16 ; 5. period = 2nlpc, where 
- (dbp 1 + 1 ) sin p(b —a) = p(b-a) cos p(fc-o) ; 7. 1GG per sec.] 
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ELECTRIC WAVES * 


§ 68. Before we discuss the propagation of electric waves, 
we shall summarise the most important equations that we 
shall require. These are known as Maxwell’s equations. 
Let the vectors E (components E x , E v , E z ) and H (com- 
ponents H x , H v , H z ) denote the electric and magnetic 
field strengths. These are defined | as the forces on a 
unit charge or pole respectively when placed inside a 
small needle-shaped cavity, the direction of the cavity 
being the same as the direction in which we wish to measure 
the component of E or H. We shall suppose that all our 
media aro isotropic with no ferromagnetism or permanent 
polarisation ; thus, if wo write e for the dielectric constant, 
and [x tor the permeability, then the related vectors, 
viz. the magnetic induction B and the dielectric dis- 
placement D are given by the equations B = pH, D = eE. 
further, let J (components j x , j y , j,) denote the current 
density vector, and p the charge density. Then, if we 
measure j, B and H in electromagnetic units, E and D 
in electrostatic units, writing c for the ratio between the 
two so ts of units, * Maxwell’s equations may bo summarised 
in vector form as follows : 


div D = inp (1) 

(,iv B = 0 (2) 

* Before reading this chapter, the student is advised to 
ami mrise nmsolf with tho equations of olectromagnetisin, as 
round in text hooks such ns those by Joans. Pidduck. Abraham- 
neckor or the* present writer. 

! £'‘ e ’ Coulson - Eleclricili/, Oliver and Boyd, Chapter II, VI. 

, This system is known ns tho Mixed System. If we had used 

entirdy e.s.u.. or entirely o.m.u.. tho powers of c would have been 

different. Pnrticiilur enro is required in discussing tho units in 

■ ) and (7). In this chapter c will always denote the ratio of tho 
two sot« of units. 
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1 dD 

curl H = 47rJ H — . . . . (3) 

c ct 

IdB 

curl E = — . . . . (4) 

c dt 

D = eE (5) 

B = /xH. . . (0) 

To these equations we must add the relation between J 
and E. If a is the conductivity, which is the inverse of 
the specific resistance, this relation is 

] = ctE . .... (7) 

For conductors a is large, and for insulators it is small. 

Tho above equations have been written in vector form ; 
until the student has acquired familiarity with the use of 
the vector notation and operation, he is advised to verify 
tho various calculations of this chapter, using the equations 
in Cartesian form as well as vector form. This will 60 on 
show how much simpler the vector treatment is, in nearly 
every case. If we wish to write these equations in their 
full Cartesian form, we have to remember that 

dD z dD * , 8D > 1 .u * 

div D = y . D = — - - + — ’ + — and that 

dx Gy dz 


„ (dH, cH y cU t ill, cH v 
curlH SV xH = |— 

The preceding equations then become 


_ cJIA 

<y 1 


dD x 

dD v 

dD, 


(1') 

dB, 

dB„ 

+ 

dB, 

A 

dx 

V 

d'J 

4- = 

dz 

477/3 

dx dy 

dz 

= u 

dH, 

dHy 

A 2 1 

i dD x 


dE, 

dEy 


1 

cB z 

dy 

dz 

= 4t r; z + 

c dt 


dy 

dz 


c 

ct 

dH x 

dH, 


1 dD. 


dE, 

dE, 


1 

dB y 

dz 

dx 

= ±”jv + 

V 

c dt 

■ (3') 

dz 

dx 

= 

c 

17 

dHy 

dll z 


l dD, 


cE, 

cEx 


1 

dB, 

dx 

dy 

= 477;,+ 

c dt 

< 


dx 

dy 


c 

dt 

4- 
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D x = eE x , D v = eE v , D t = eEj . . (5') 

E x = v-Hx , By = tiHy , B t = jiH t . . (O') 

jx = oE x , jy = oEy , j t = oE, . . (7') 

Equations (l)-(4) aro sometimes called Maxwell’s 
Equations and equations (5)-(7) constitutive relations. 
Simple physical bases can easily bo given for (l)-(4). 
Thus, (1) represents Gauss’ Theorem, and follows from the 
law of force between two charges ; (2) represents the fact 
that isolated magnetic poles cannot be obtained ; (3) is 
Ampfere’s Rule that the work done in carrying a unit pole 
round a closed circuit equals 4 tt times the total current 
enclosed in the circuit ; part of this current is tho conduc- 
tion current J and part is Maxwell’s displacement current 
1 dD 

— — ; (4) is Lcnz’s law of induction. 

47TC Ct 

These seven equations represent the basis of our 
subsequent work. They need to be supplemented by a 
statement of the boundary conditions that hold at a change 
of medium. If suffix n denotes the component normal to 
the boundary of the two media, and suffix s denotes tho 
component in any direction in tho boundary plane, then 
on passing from the one medium to the other 

D n , B n , E, and H, are continuous . . (8) 

In cases where there is a current sheet (i.e. a finite 
current flowing in an indefinitely thin surface layer) some 
of theso conditions need modification, but wo shall not 
discuss any such cases in tliis chapter. 

There are two other important results that wo shall 
use. First, wo may suppose that the electromagnetic 
field stores energy, and the density of this energy per unit 
volume of the medium is 

^{cE*+/*H*} . . . (9) 

Second, there is a vector, known as the Poynting 
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vector, which is concerned with the rate at which energy 
is flowing. This vector, whose magnitude and direction 
are given by 

^(ExH) ( 10 ) 

represents the amount of energy which flows in unit time 
across unit area drawn perpendicular to it. E and H are 
generally rapidly varying quantities and in such cases it is 
the mean value of (10) that has physical significance. 


§69. We shall first deal with non-conducting media, 
such as glass, so that we may put a = 0 in (7) ; wo suppose 
that the medium is homogeneous, i.e. e and p are constants. 
If, as usually happens, there is no residual charge, wo may 
also put p = 0 in (1), and with these simplifications, 
Maxwell’s equations may bo written 


div E = 0 , div H = 0, 

, p8H <cE 

cur , E= _e_, curlH = -- 



These equations load immediately to the standard equation 
of wave motion, for we know * that 

curl curl H = grad div H— y 2 H. 

Consequently, from the fourth of the equations in (11), 
wo find 

c , cE e d . 

grad div H-y 2 H = - curl — = - ^ ™ rl E 


Substituting for div H and curl E, we 
equation 



discover the standard 


• Rutherford, Vector Methods, p. 69, equation (10). 


(12) 
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Eliminating H instead of E we find the same equation for E . 



According to our discussion of this equation in Chapter I, 
this shows that waves can be propagated in such a medium, 
and that their velocity is c/y/(€fi). In free space, where 
e = n = 1, this velocity is just c. Now c, which was 
defined as the ratio of the two sets of electrical units, 
has the dimensions of a velocity, and its magnitude can 
be obtained experimentally; it is approximately 2-998 . 10 10 
cms. per sec. But it is known that the velocity of light 
in free space has exactly this same value. We are thus 
led to the conviction that light waves are electromagnetic 
in nature, a view that has subsequently received complete 
verification. X-rays, y-rays, ultra-violet waves, infra-red 
waves and wireless waves are also electromagnetic, and 
differ only in the order of magnitude of their wavelengths. 
We shall be able to show later, in § 71, that these waves 
are transverse. 

In non-conducting dielectric media, like glass, e is not 
equal to unity ; also jz depends on the frequency of the 
waves, but for light waves in the visible region we may 
put p = 1 . The velocity of light is therefore c/\/c- Now 
:n a medium whoso refractive index is K, it is known 
experimentally that the velocity of light is c/K. Hence, 
if our original assumptions are valid, e = K 2 . This is 
known as Maxwell's relation. It holds good for many 
substances, but fails because it does not take sufficiently 
detailed account of the atomic structure of the dielectrio. 
It applies better for long waves (low frequency) than for 
short waves (high frequency). 


§ 70. A somewhat different discussion of (11) can bo 
given in terms of the electric and magnetic potentials. 
Since div B = 0, it follows that wo can write 

B — /zH = curl A, . (14) 
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where A is a vector yet to be determined. 'I his equation 
does not define A completely, since if i/r is any scalar, 
curl (A+grad ip) = curl A. Thus A is undefined to the 
extent of addition of the gradient of any scalar, and we 
may accordingly impose one further condition upon it. 

l 

If B = curl A, and curl E = - - — , it follows, by 

C Cl 


elimination of B, that 


Integrating, 


i idA 

curl Ed — 

I c Cl 


= 0 . 


IdA . , 

E+ - — = - grad <f>, 
c ct 


where <f> is an arbitrary scalar function. 


i.e. 



, , ldA 

*~lu 



In cases where there is no variation with the time this 
becomes E = — grad <f>, showing that <f> is the analogue 

of the electrostatic potential. 

Eliminating H from the relations /iH = curl A, 

curl H = - and using (15) to eliminate E, we find 
c dt 


ue d(f> €fid- A 

grad div A— y 2 A = — grad — — ^ • 


i.e. 


fiia 2 A , H 

v A = # a^ + grad l dlvA+ 7 « 


I.ct us now introduce the extra allowed condition upon A, 

and write a . 

divA+^^-0 . ■ • (18) 

c ct 

Then A satisfies the standard equation of wave motion 

€/ x a 2 A 
V A “ c* dt- 


( 17 ) 
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Further, taking the divergence of (15), we obtain, by (16) 
0 = div E = -vV- i a - div A = -vV+ ^ ^ • 


Tlius <f> also satisfies the standard equation 




A similar analysis can be carried through when p and J 
are not put equal to zero, and we find 


B = curl A . 

. (14') 

_ . . 10A 

E = ^grad 4,--- . 

• (15') 

0„divA+iL <a i . 

c cl 

. (16') 


. (17') 

_ «/* ^ i~P 

V V c 2 et 2 e 

. (18) 


<f> and A are known as the electric potential and magnetic 

or vector potential respectively. It is open to our 

choice whether we solve problems in terms of A and <f>, 

or of E and H. The relations (14')-(18') enable us to pass 

from the one system to the other. The boundary condi- 

tions for tf> and A may easily be obtained from (8), but 

since wo shall always adopt the E, H typo of solution, 

" uch is usually the simpler, there is no need to write 
them down here. 

There is one other general deduction that can be made 
lore. wo use (3), (5) and (7) we can write, for 
homogeneous media, 

curl H = 4 ttoE + - — . 

C 8t 
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Taking the divergence of each side, and noting, from (1), 
that div E = 4 t rp/e, we find 

c dp 

—£ rap = 0. 

c cl 

Thus, on integration, 

p = p Q e~ 'M, where 9 = el±-noc . . (19) 

9 is called the time of relaxation. It follows from (19) 
that any original distribution of charge decays exponentially 
at a rate quito independent of any other electromagnetic 
disturbances that may be taking place simultaneously, 
and it justifies us in putting p = 0 in most of our 
problems. With metals such as copper, 9 is of the order 
of 10 _1S secs., and is beyond measurement ; but with 
dielectrics such as water 9 is large enough to be deter- 
mined experimentally. Equation (19) only applies to the 
chargo at an internal point in a medium ; charges at the 
boundary of a conductor or insulator do not obey this 
equation at all. 

§71. We next discuss plane waves in a uniform non- 
conducting medium, and 6how that they are of transverse 
typo, E and H being perpendicular to the direction of 
propagation. Let us consider plane waves travelling with 
velocity V in a direction l, m, n. Then E and H must bo 
functions of a new variable 

u = Ix+my+nz — Vt . . • (20) 

When wo say that a vector such as E is a function of u, 
wo mean that each of its three components separately 
is a function of u, though the three functions need not 
be the same. Consider the fourth equation of (11). I* s 
x-component (see (3')) is 

BH, ML, _ € dE, 

dy dz c dt 



110 


WAVES 


If dashes denote differentiation with respect to u, this is 


mHz —nH y ' = - - E,'. 

c 


Integrating with respect to u, this becomes 


mH x —vH y = E x , 

c 


in which we have put the constant of integration equal 
to zero, since we are concerned with fluctuating fields 
whoso mean value is zero. There are two similar equations 
to the above, for E y and E„ and we may write them as 
one vector equation. If wo let n denote the unit vector 
in the direction of propagation, so that n = (l, m, n), we 
havo 

eV 

nxH E (21) 

c 


Exactly similar treatment is possible for the third equation 
of (11) ; wc get 

llV 

n x E = — H . . . . (22) 

c 


Equation (21) shows that E is perpendicular to n and H, 
and (22) shows that H is perpendicular to n and E. In 
other words, both E and H are perpendicular to the direc- 
tion of propagation, so that the waves are transverse, and 
in addition, E and H are themselves perpendicular, E, 
H and n forming a right-handed set of axes. If we 
eliminate H from (21) and (22) and use the fact that 

n x (n x E] = (n . E)n-(n . n)E = -E, 

since n is perpendicular to E and n is a unit vector, we 
discover that P = c 2 ‘€fi, showing again that the velocity 
of these waves is indeed c[\/(en). 

It is worth while writing down the particular cases of 
(21) and (22) that correspond to plane harmonic waves 
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in the direction of the z axis, and wit h the E vector in the 
x or y directions. The solutions are 


E t = 0 

H x = -v/(e/p)ae IJ *‘-* /r) 


U 

V 

5 

a 

11 

a 

H v = 0 

(23) 

E t = 0 

H z = 0. 


E x = be ipV ~ tlV) 

H z = 0 


E v = 0 

11 v = +V(e//x)5e'>'‘-" F » 

(24) 

E, =0 

//, = 0. 



In accordance with § 10, a and b may bo complex, the 
arguments giving the two phases. It is the general 
convention * to call the plane containing H and n the 
plane of polarisation. Thus (23) is a wave polarised 
in the xz plane, and (24) a wave polarised in the yz plane. 
By the principle of superposition (§ 6) we may superpose 
solutions of types (23) and (24). If the two phases are 
different, we obtain elliptically polarised light, in which 
the end-point of the vector E describes an ellipse in the 
xy plane. If the phases aro the same, wo obtain plane 
polarised light, polarised in the plane y/x = —b/a. If 
the phases differ by tt/2, and the amplitudes are equal, 
wo obtain circularly polarised light, which, in real form, 
may be written 

E x = a cos p(t—z/V) H z = — V( € lp) a s ‘ n ) 

E v = a sin p(t-z/V) //„ = + V(«//*) a 308 P( l ~ z ! V ) 
E t = 0 //, = 0. (25) 

The end-points of the vectors E and H each describe 
circles in the xy plane. 

In all three cases (23)-(25), when wo are dealing with 
free space (c = p = 1) the magnitudes of E and H are 
equal. 

§ 72. By the use of (10) wo can easily write down the 
rate at which energy is transmitted in these waves. Ihus, 


• To which, unfortunately, not all writers conform. 
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with (25) the Poynting Vector is simply (0, 0, — 

\ 47 r 

This vector is in the direction of the positive z-axis, showing 
that energy is propagated with the waves. According to 
(9), the total energy per unit volume is 



s' fE,+ ' JP > 



From these two expressions we can deduce the velocity 
with which the energy flows ; for this velocity is merely 
the ratio of the total flow across unit area in unit time 
divided by the energy per unit volume. This is c/\/ (cp), 
so that the energy flows with the same velocity as the 
wave. This does not hold with all types of wave 
motion ; an exception has already occurred in liquids 
(§ 52). 

When wo calculate the Poynting Vector for the waves 

£ 

(23) and (24), wo must remember that — ExH is not a 

4tt 

linear function and consequently (see § 10) we must choose 
either the real or the imaginary parts of E and H. Taking, 
for example, the real part of (23), the Poynting Vector 
lies in the z direction, with magnitude 

V?. . L a’cos^-i), 


This is a fluctuating quantity whoso mean value with 
respect to tho time is — yj~. The energy density, from 
• € ® a 

(9)> is cos 2 p(/— z/V), with a corresponding mean value 

ea 2 /87T. Once again the velocity of transmission of energy 
. ca J /e ea 2 

18 87tV/I 8tt = w h‘ c h is the same as the 

wave velocitv. 
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§ 73. We shall next discuss the reflection and refraction 
of plane harmonic light waves. This reflection will be 
supposed to take place at a plane surface separating two 
non-conducting dielectric media whose refractive indices are 
K x and K 2 . Since we may put = /i 2 = 1 , the velocities 
in the two media are c/K v cjK t . In tig. 21 let Oz be the 


x 



direction of the common normal to the two media, and lot 
AO, OB, OC be the directions of the incident, reflected 
and refracted (or transmitted) waves. Wo have not yet 
shown that these all lio in a plane ; lot us suppose that 
they make angles 0, n—6' and $ with tho z axis OA being 
in the plane of the paper, and let us take the piano ol 
incidence (i.e. the plane containing OA and Oz) to be the 
12 plane. The y axis is then perpendicular to the piano 
of tho paper. 

H 
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Since the angle of incidence is 0, then as in (20), each 
of the three components of E and H will be proportional to 

g}P{ct -KtizslnO+zcosQ))' 

Let the reflected and transmitted rays move in directions 
(l v m v n,) and (l 2 ,m 2 , n 2 ) respectively so that n, = —cos 9', 
n 2 = cos <J>. Then the corresponding components of E 
and H for these rays will be proportional to 

gip{rf an( j g<p{rt -K,dtZ+ m,v + n,z)) 

Thus, considering the E x components, wo may write the 
incident, reflected and transmitted values 

A e >p{'i-K l (x sin 8+z cos &)) , ^^{pfet-KtOp+m,!/ + »»,*)) 

A n C. ip ' ri ~ K tU,z+m, v + n^)). 

These functions all satisfy the standard equation of wave 
motion and they have the same frequency, a condition 
which is necessary from the very nature of the problem. 

We shall first show that the reflected and transmitted 
waves lie in the plane of incidence. This follows from tho 
boundary condition (8) that E x must be continuous on the 
plane 2 =0, i.e. for all x, y, t, 

Aeipfrt-K.x Bln 0) _J_ J i e«p{d-/r l a 1 x+m,y)} _ Ajjrttt-Ktoe+mMl 

This identity is only possiblo if the indices of all three 
terms are identical : i.e. 

ct-A> sin 6 = d-K^x+m#) = d-K^+m#). 

Thus sin 9 = Ay, = Ay,, 

0 = K x m x = K 2 m 2 . 

The second of these relations shows that m, = m 2 = 0, 
so that the reflected and transmitted rays OB, OC lie in 
the plane of incidence xOz. The first relation shows that 
= sin 9, i.e. that the angle of reflection 9' is equal to 
tho angle of incidence 9, and also that 

A', sin 9 = K a sin (f> . . . (26) 
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This well-known relationship between the angles of 
incidence and refraction is known as Snell’s law. 

Our discussion so far has merely concerned itself with 
directions, and wo must now pass to the amplitudes of 
the waves. There are two main cases to consider, according 
as tho incident light is polarised in the plane of incidence, 
or perpendicular to it. 

Incident light polarised in the plane of incidence . — The 
incident ray AO lias its magnetic vector in the xz plane, 
directed perpendicular to AO. To express this vector in 
terms of x, y, z it is convenient to use intermediary axes 
£, tj, £ through 0 (see fig. 22, where the directions of £ 
and £ arc shown ; t) coincides with the y axis which is 



perpendicular to the plane of the paper). £ is in the 
direction of propagation, and £ is in the plane of incidence. 
Referred to these new axes, H lies entirely in the £ direction, 
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and E in the 77 direction. We may use (23) and write 

Ef =Et=0,E v = 

H v = Hi = 0 , Hf = -Kja,!***-*#. 

Now £ = £ sin 0 -f "2 cos 9, and so it follows that : 
incident wave 

E x = 0 , H x = —K l a l cos 9 e'rt rf - ^,(* sin 0 +* eoa 0)}, 

E v = a l e ip ^ ct ~ K ^ xslD ^+ iCos ^) j H = 0 

E t = 0 , H z = K 1 a l sin 9 e ip ^ ct ~ K ^ zsl "^+ tcos ^). 

Similar analysis for the reflected and refracted waves, in 

which wo replace 9 by ir~9 and <f> in turn, enables us 
to write 

reflected wave 

E x = 0 , H z — K x b x cos 9 e'rt'i-ZAzdnO-ico&O)), 

E v = 6 1 c ,p!rt-ff ‘ (Isin ^ - * co *0)) t H = 0 
E 2 = 0 , H t = K x b x sin 9 e»*to-*i<*sin0-*co80)} . 
refracted wave 

E x = 0, H z = — ^ 2 a 2 cos^e'>{rt-x,(* 6 in^+*cos^)} ( 

E u = a 2 e>p{rt-^.(j-3in^+tcos^)) t H v = 0, 

E z = 0 , H t = K t a j sin <f> e'X<*-.ff 1 (xsin^+*co 9 ^)} i 

Wo may write the boundary conditions in the form that 
Ey,'K 2 E s , Ii x , H v and H z are continuous at z = 0. 
These six conditions reduce to two independent relations, 
which we may take to bo those due to E y and H x : 

a i+ 6 i = a s , 

~ K i a i C03 0+K x b x cos 9 = -K<a z cos f 

Thus 


a 




a. 


K x cos 9+K a cos f> ^ K x cos 0-K z cos <f> ~ 2K X cos 0 ' 

Using Snell’s law (26) in the form K x : K 2 = sin d> : sin 9, 
this gives 


a, 


61 


sin (9+cf>) -si n(9- < f>) 2 sin <f> cos 9’ 


(27) 
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Equation (27) gives the ratio of the reflected and refracted 
amplitudes. If medium 2 is denser than medium 1, K Z >K V 
so that 9><f>, and thus b 1 /a l is negative ; so there is a 
phase change of i t in the electric field when reflection takes 
place in the lighter medium. There is no phase change on 
reflection in a denser medium, nor in the refracted wave. 
The same conclusion is true for the magnetic field H, though 
it must be remembered here and later in this chapter that 
the positive directions for E and H are defined by %r)l, 
in fig. 22 and their counterparts relative to OB, OC, 
the rj direction being along the y axis throughout. 

Incident light polarised perpendicular to the plane of 
incidence. — A similar discussion can be given when the 
incident light is polarised perpendicular to the plane of 
incidence ; in this case the roles of E and H are practically 
interchanged, H y for example being the only non-vanishing 
component of H. It is not necessary to repeat the analysis 
in full. With the same notation for the amplitudes of the 
incident, reflected and refracted waves, we have 


1 V 1 u 2 

sin 29+ sin 2<f> sin 20— sin 2<f> 4cosds\n<f> ^ ^ 

It follows from (28) that the reflected ray vanishes if 
sin 26 = sin 2 <f>. Since 6 ^ tf>, this implies that 9-\-<f> = tt/2, 
and then Snell’s law gives 



K j sin 6 = K 2 sin <f> = K 2 cos 9, 

tan 9 = KJK l = y/(* 2 l*i) • • (29) 


With this angle of incidence, known as Brewster’s angle, 
there is no reflected ray. 

There will be a phase change of n in E on reflection 
when sin 2 9< sin 2 <f>. If reflection takes place in the 
lighter medium so that ATj<A 2 , this holds for 9 greater 
than Brewster’s angle ; but if K l >K 2 it holds for 9 less 
than Brewster’s angle. In all other cases there is no phase 
change on reflection. 
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In general, of course, the incident light is composed of 
waves polarised in all possible directions. Equations (27) 
and (28) show that if the original amplitudes in the two 
main directions are equal, the reflected amplitudes will 
not be equal, so that the light becomes partly polarised 
on reflection. When the angle of incidence is given by 
(29) it is completely polarised on reflection. This angle is 
therefore sometimes known as the polarising angle. 

§74. An interesting possibility arises in the discussion 

of § 73, which gives rise to the phenomenon known as 

total or internal reflection. It arises when reflection 

takes place in the denser medium so that <f>>8. If we 

suppose 8 to be steadily increased from zero, then <f> also 

increases and when sin 8 = Ay A', , <f> = tt/ 2. If 8 is 

increased beyond this critical value, <j> is imaginary. 

lhere is nothing to disturb us in this fact provided that 

we interpret the analysis of § 73 correctly, for we never 

had occasion to suppose that the coefficients were real. 

We can easily make the necessary adjustment in this 

case. Take for simplicity the case of incident light 

polarised in the plane of incidence. Then the incident 

and reflected waves are just as in our previous calculations. 

I he refracted wave has the same form also, but in the 

exponential term, A' 2 sin <f> = A', sin 8 , and is therefore 
real, whereas 

A, cos <£ = \/(Ay-Ay sin 2 <£) = % /(Ay — Ay sin 2 8), 

and is imaginary, since we are supposing that internal 
reflection is taking place and therefore A', sin 8>K 2 . We 
may therefore write A\ cos <f> = ± i q , where q is real, 
thus the refracted wave has the form 

E v = (l 2 e'rtct ' KtSlnd z± iqz) 

= a^-r<P e ip(ct^K x 8ln0*) # 

lor reasons of finiteness at infinity, we have to choose 
the negative s.g so that it appears that the wave is 
attenuated as it proceeds into the less dense medium. 
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For normal light waves it appears that the penetration 
is only a few wavelengths, and this justifies the title 
of total reflection. The decay factor is 

e -m _ e -;> v /(K 1 , 6!n , 0-/r 1 , )z 

This factor increases with the frequency so that light of 
great frequency hardly penetrates at all. In actual 
physical problems, the refractive index docs not change 
from K x to K 2 abruptly, as we have imagined ; however, 
Drude has shown that if we suppose that there is a thin 
surface layer, of thickness approximately equal to one 
atomic diameter, in which the change takes place smoothly, 
the results of this and the preceding paragraphs are hardly 
affected. 


§ 75. In our previous calculations we have assumed 
that the medium was non-conducting, so that we could 
put o = 0. When wo remove this restriction, keeping 
always to homogeneous media, equations (l)-(7) give us 

div E = 0, 
div H = 0, 

e 8E 

curl H = 4 wctE H — — , 

c 8t 

, /^H 

curl E = — • 

c cl 


Now curl curl E = grad div E— y 2 E = — y 2 E, so that 


„„ M , cH n d , „ 4 nop cE , «/* d ' E 

V 2 E = - curl — = - - curl H = — + — -rz , 

v c 81 c 81 c Cl c- Cl 2 





4ttoh c E 

c Hi 



A similar equation holds for H. Equation (30) is the 
well-known equation of telegraphy (see § 9). The first 
term on the right-hand side may be called the displacement 

1 8D 


term, since it arises from the displacement current — — 
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and the second is the conduction term, since it arises from 
the conduction current j. If we are dealing with waves 
whose frequency is pj2ir, E will be proportional to e ipt ; 
the ratio of these two terms is therefore epfimeo. Since 
e is generally of the order of unity, this means that if 
pl2ir is much greater than ca, only the displacement term 
matters (this is the case of light waves in a non-conducting 
dielectric) ; but if pfi-n is much less than cc, only the 
conduction term matters (this is the case of long waves 
in a good metallic conductor). In the intermediate region 
both terms must be retained. With most metals, if p< 10 7 
per second we can neglect the first terra, and if p>10 15 
per second we can neglect the second terra. 

Let us discuss the solutions of (30) which apply to 
plane harmonic waves propagated in the z direction, such 
that only E x and H v are non-vanishing (as in (24)). We 

may suppose that each of these components is pro- 
portional to 

eW-rt .... (31) 

where p/ 2 tt is the frequency and q is still to be determined. 
This expression satisfies the equation (30) if 



q is therefore complex, and we may write it 


where 




The " velocity ” of (31) is l/ g ; but we have seen in § 73 
that in a medium of refractive index K the velocity is c/K. 
So the effective refractive index is cq which is complex. 
Complex refractive indices occur quite frequently and are 
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associated with absorption of the waves ; for, combining 
(31) and (33) we have the result that E x and //, are 
proportional to 

e -Pp e <p(t-a*) . . . (34) 

This shows that a plane wave cannot be propagated in 
such a medium without absorption. The decay factor may 
be written e~ kt where k = pfi. k is called the absorption 
coefficient. In the case where i-nac/ep is small compared 
with unity (the case of light waves in most non-conducting 
dielectrics), k is approximately equal to 2tt ay/ip/e). Now 
the wavelength in (34) is A = 2i r/ap, so that in one wave- 
length the amplitude decays by a factor e - ^, approximately 
e-4j r’ac/«p. As we are making the assumption that cajcp 
is small, the decay is gradual, and can only be noticed 
after many wavelengths. The distance travelled before 
the amplitude is reduced to 1/e times its original value is 
l/k, which is of the same order as a. 

The velocity of propagation of (34) is 1/a, and thus 
varies with the frequency. With our usual approximation 
that co/cp is small, this velocity is 



We can show that in waves of this character E and H 
are out of phase with each other. For if, in accordance 
with (31), we write 

E x = a ««»-•*>, 

//„ = b «<**-«•>, 

then the y-componcnt of the vector relation 

curl E = — , 

c cl 

gives us the connection between a and b. It is 

cE x p. c7/„ 
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ie. gra = - 6 . . . . (36) 

c 

Thus 6/a is equal to (c/fi)q. Now q is complex and hence 
there is a phase difference between E x and U v equal to 
the argument of q. This is tan -1 (/3/a), and with the same 
approximation as in (35), this is just tan _1 (27r<7c/ejj)» which 
is effectively 2/mjclep. 

§ 76. It is interesting to discuss in more detail the 
case in which the conductivity is so great that we may 
completely neglect the displacement term in (30). Let 
us consider the case of a beam of light falling normally 
on an infinite metallic conductor bounded by the plane 
z = 0. Let us suppose (fig. 23) that the incident waves 



come from the negative direction of z, in free space, for 
which e — p — 1, and are polarisod in the yz plane. Then, 
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according to (24) they are defined by : 
incident wave 

E x = a x , H„ = e'X‘-* w . 

reflected wave 

E x = b, e«» + *W , H v = -b, e^ t+tle) . 


In the metal iteelf we may write, according to (31) and (30), 


E x = o 2 , H,= - qa t 

These values will satisfy the equation of telegraphy (30) 
in which we have neglected the displacement term, if 

“W.- _ o. a: 


8* = - -Z? i = “Vb 

pc 


where y 2 = 2-no pi pc. 

q = y(l-i). . 


Thus 


• • • 


(37) 


Inside the metal, E and H have a v\ 4 phase difference, 
since, as we have shown in (30), this phaso difference is 
merely the argument of q. 

The boundary conditions are that E x and II v are 
continuous at z = 0. This gives two equations 


Henoe 


** 1+^1 — a 2 > 




(38) 


Since q is complex, all three electric vectors have phase 
differences. The ratio It of reflected to incident energy 
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is | bja^ 2 , which reduces to 

(cy-fi) 2 + (cy) 2 
(cy+/x) 2 + (cy) 2 ' 

In the case of non-ferromagnetic metals, cy is much larger 
than /x, so that approximately 

cy 

This formula has been checked excellently by the experi- 
ments of Hagen and Rubens, using wavelengths in the 
region of 10~ 6 cms. 

It is an easy matter to generalise these results to apply 
to the case when we include both the displacement and 
conduction terms in (30). 

We can use (38) to calculate the loss of energy in the 
metal. It we consider unit area of the surface of the metal, 
the rate of arrival of energy is given by the Poynting Vector. 


T^ 8 * s I a i I 2, Similarly the rate of reflection of energy is 


— | b x | a . So the rate of dissipation is j | a x | a — | | 2 j. 


This must be the same as the Joule heat loss. In our 
units, this loss is caE 2 per unit volume per unit time. 
It we take the mean value of E x 2 in the metal, it is an 

easy matter to show that J* caE x 2 dz is indeed exactly 

equal to this rate of dissipation. 


§ 77. When the radiation falls on the metal of § 76, it 
exerts a pressure. V e may calculate this, if we use the 
experimental law that when a current J is in the presence 
ol a magnetic field H there is a force fi jxH acting on 
it In our problem, there is, in the metal, an alternating 
Geld E, and a corresponding current oE. The force on the 
current is therefore /xaE X H, and this force, being perpen- 
dicular to E and H, lies in the z direction. The force on 



ELECTRIC WAVES 


125 


the charges that compose the current is transmitted by 
them to the metal as a whole. Now both E and H are 
proportional to e~ p ^ z (see equation 37) so that the force 
falls off according to the relation e~ 2pY *. To calculate the 
total force on unit area of the metal surface, we must 
integrate /zoExH from z = 0toz = co. ExHisa 
fluctuating quantity, and so we shall have to take its mean 
value with respect to the time. The pressure is then 

i.e. (ca/ip) \a 2 \*. 

Using (38) this may be expressed in the form 

(c<j[i 2 lp) | PJlcy+^+M 2 )’ 1 . 

§ 78. There is another application of the theory of 
§ 76 which is important. Suppose that we have a straight 
wire of circular section, and a rapidly alternating c.m.f. 
is applied at its two ends. We have seen in § 76 that 
with an infinite sheet of metal the current falls off as wo 
penetrate the metal according to the law e" pyt - If py 
is small, there is little diminution as wo go down a distance 
equal to the radius of the wire, and clearly the current 
will be almost constant for all parts of any section (see, 
however, question (12) in §79). But if py is large, then 
the current will bo carried mainly near the surface of the 
wire, and it will not make a great deal of difference whether 
the metal is infinite in extent, as we supposed in § 76, or 
whether it has a cross-section in the form of a circle ; in 
this case the current density falls off approximately 
according to the law e“ pyr as we go down a distance r 
from the surface. This phenomenon is known as 
the skin effect ; it is more pronounced at very high 
frequencies. 

We could of course solve the problem of the wire quite 
rigorously, using cylindrical polar coordinates. I lie 
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formulae are rather complicated, but the result is 
essentially the same. 

§ 79. Examples 

(1) Prove the equations (17') and (18') in § 70. 

(2) Find the value of H when E a = E y = 0, and E t = 
A cos nx cos net. It is given that H = 0 when t = 0, and also 
€ = /x = 1, p = a = 0. Show that there is no mean flux of 
energy in this problem. 

(3) Prove the equation (28) in § 73 for reflection and 
refraction of light polarised perpendicular to the plane of 
incidence. 

(4) Show that the polarising angle is less than the critical 
angle for internal reflection. Calculate the two values if 
K x = 6, K t = 1. 

(5) Show that the reflection coefficient from gloss to air at 
normal incidence is the samo as from air to gloss, but that 
the two phase changes are different. 

(6) Light falls normally on the plane face which separates 
two media K x , K r Show that a fraction R of the energy is 
reflected, and T is transmitted, where 

R / *.-*» \» , 4Ayv, 

\K t +Kj ’ ‘ 

Hence prove that if light fulls normally on a slab of dielectric, 
bounded by two parallel faces, the total fraction of energy 

reflected is — -rrr* anc * transmitted is ■ . it ^ 

A,*+A,* A',*+AV 

necessary to bike account of the multiple reflections that take 
place at each boundary. 

(7) Light passes normally through the two parallel faces 
of a piece of pinto glass, for which K = 1-5. Find the fraction 
of incident energy transmitted, taking account of reflection 
at the faces. 

(8) Show that when internal reflection (§ 74) is taking 
place, there is a phaso change in the reflected beam. Evaluate 
this numerically for the case of a beam falling at an angle 
of 60° to the normal when A', = 1-6, K t = 1, the light being 
polarised in the piano of incidence. 

(9) Show if n = 1. then the reflection coeflicient with 
rnduls (§ 7b may be written in the form R = l—2/^/(ca/v), 
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where v is the frequency. If o is 1-6 . 10 7 (in our mixed units), 
calculate R for A = 10 -3 cms. and A = 10- 1 cins. 

(10) A current flows in a straight wire whose cross-section 
is a circle of radius a. The conduction current j depends 
only on r the radial distance from the centre of the wire, 
and the time t. Assuming that the displacement current 
can be neglected, prove that H is directed perpendicular to t ho 
radius vector. If j(r 9 /) and II(r 9 t) represent the magnitudes 
of j and H, prove that 


d . 8j h<j cH 

— (Hr) = 4nrj , ^ . 

Cr Cr c Ct 


(11) Use the results of question ( 10) to prove that j satisfies 
the differential equation 

47t<jH dj 

C it 


dj\ 

r dr \ dr) 


By using the formula for curl in cylindrical polars twice in 
succession show that H[ = H e (r t /)] satisfies the equation 

d*H 1 CH H infia ill 

— = [curl curl H ^ 

dr 2 ^ r dr r* 1 c it 


Use the method of separation of variables to prove that 
there is a solution of the ^-equation of the form j = 
where 

dV 1 If 

jl + " -I * ia f = 0 , a = iranp/c. 

ar z r dr 


Hence show that / is a combination of Bessel functions of 
order zero and complex argument. 

(12) If a in question (11) is small, show that an approximate 
solution of the current equation is j = A( 1 -f far 2 — Ja-/ I )t ,p< , 
where A is a constant. Hence show that the total current 
fluctuates between ±J, whore, neglecting powers of a above 
the second, J = *aM( l+a 4 a*/24). Use this result to show 
that the heat developed in unit length of the wire in unit 

time is -- (1-f a 4 a 2 /12). (Questions (10), (11) and (12) 

are the problem of the skin effect at low frequencies.) 


[Answers : 2. //, = //, = 0 , II y = —A sin nx sin net ; 
4. 9° 28', 9° 36'; 7. 12/13 of the incident energy is trans- 
mitted ; 8. 100° 20' ; 9. 0-984, 0-950.J 
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GENERAL CONSIDERATIONS 

§ 80. The speed at which waves travel in a medium is 
usually independent of the velocity of the source ; thus, if a 
pebble is tlirown into a pond with a horizontal velocity, the 
waves travel radially outwards from the centre of disturb- 
ance in the form of concentric circles, with a speed which is 
independent of the velocity of the pebble that caused them. 

When wo have a moving source, sending out waves 
continuously as it moves, the velocity of the waves is 
often unchanged,* but the wavelength and frequency, as 
noted by a stationary observer, may be altered. 

Thus, consider a source of waves moving towards an 
observer with velocity «. Then, since the source is moving, 


ntx 



(a) Waves whoa source is stationary. 

(b) Waves when source is moving. 


the waves which are between the source and the observer 
will be crowded into a smaller distance than if the source 
had been at rest. This is shown in fig. 24, where the waves 
are drawn both for a stationary and a moving source. If 
the frequency is n, then in timo t the source emits nt waves. 

• It is changed slightly when there is dispersion ; see § 83. 

128 
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If the source had been at rest, these waves would have 
Occupied a length AB. But due to its motion the source 
has covered a distance id, and hence these nl waves are 
compressed into a length A'B', where AB—A'B' = ut. 
Thus 

nl\—nt\' = ut, 

i.e. A' = A— u/n = A(l— u/c), .... (1) 

if c is the wave velocity. If the corresponding frequencies 
measured by the fixed observer are n and n', then, since 
nA = c = n'A', therefore 


n = 


TIC 

c—u 


( 2 ) 


If the source is moving towards the observer the frequency 
is increased ; if it moves away from him, the frequency is 
decreased. Tliis explains the sudden change of pitch 
noticed by a stationary observer when a inotor-car passes 
him. The actual change in tliis case is from nc/(c— u) to 
ncl(c- t-u), so that 

An = 2ncu/(c 3 —u i ). . . ■ (3) 

This phenomenon of the change of frequency when a source 
is moving is known as the Doppler effect. It applies 
equally well if tho observor is moving instead of the source, 
or if both are moving. 

For, consider tho caso of the observer moving with 
velocity v away from the source, which is supposed to bo 
at rest. Let us superimpose upon the whole motion, 
observer, source and waves, a velocity —v. We shall 
then have a situation in which the observer is at rest, 
the source has a velocity — v, and the waves travel with 
a speed c—v. Wo may apply equation (2) which will then 
give tho appropriate frequency as registered by the observer; 
if this is n', then 


n m = 


n(c-v) 


tl(C~v) 


(c-v)-(-v) 


(4) 
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To deal with the case in which both source and observer 
are moving, with velocities u and v respectively, in the 
same direction, we superimpose again a velocity —v upon 
the whole motion. Then in the new problem, the observer 
is at rest, the source has a velocity u— v, and the waves 
travel with velocity c—v. Again, we may apply (2) and 
if the frequency registered by the observer is we have 



n(c—v) n(c—v) 

(c—v) — (u—v) c—u 



These considerations are of importance in acoustic and 
optical problems ; it is not difficult to extend them to 
deal with cases in which the various velocities are not in 
the same line, but we shall not discuss such problems here. 


§ 81. Wo have shown in Chapter I, § 6 that we may 
superpose any number of separate solutions of the wave 
equation. Suppose that wo have two harmonic solutions 
(Chapter I, equation (11)) with equal amplitudes and nearly 
equal frequencies. Then the total disturbance is 


<f> = a cos (Zqa:— 71,0+a cos 2ir(k^c—n^) 

kj—Ic* 


2a cos 277 ( - ■ t> 8 x — k — - t 

— z 


cos 2tt 


x— 


rti-n, 


t 


The first cosine factor represents a wave very similar to 
the original waves, whoso frequency and wavelength are 
an average of the two initial values, and which moves 

with a velocity ~ a . This is practically the same as the 

velocity of the original waves, and is indeed exactly the 
same if njk^ = n 9 /k 2 . But the second cosine factor, which 
changes much more slowly both with respect to x and t, 
may be regarded as a varying amplitude. Thus, for the 
resultant of the two original waves, wo have a wave of 
approximately the same wavelength and frequency, but 
with an amplitude that changes both with time and distance. 
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Wo may represent this graphically, as in fig. 25. The 
outer solid profile is the curve 

o o Ui— K n,— 7i, 

y = 2a COS 277 p_-* x _ _L _- t ■. 

The other profile curve is the reflection of tliis in the x 
axis. The actual disturbance cf> lies between these two 
boundaries, cutting the axis of x at regular intervals, 
and touching alternately the upper and lower profile 
curves. If the velocities of the two component waves 
are the same, so that n 1 /k l = njk 2 , then the wave system 
shown in fig. 25 moves steadily forward without change 




Fio. 23 


of shape. Tho caso when njl 'c 2 is not equal to njk t is 
dealt with in § 83. 

Suppose that <f> refers to sound waves. Then we shall 
hear a resultant wave whoso frequency is the mean of 
tho two original frequencies, but whose intensity fluctuates 
with a frequency twice that of the solid profile curve. 
This fluctuating intensity is known as beats ; its frequency, 
which is known as the beat frequency, is just 74~n 2 , 
that is, the difference of tho component frequencies. Wo 
can detect beats very easily with a piano slightly out of 
tune, or with two equal tuning-forks on tho prongs of 
ono of which wo have put a little scaling wax to decrease 
its frequency. Determination of the beat frequency 
between a standard tuning-fork and an unknown frequency 
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is one of the best methods of determining the unknown 
frequency. Low frequency beats are unpleasant to the ear. 

§ 82. There is another phenomenon closely related to 
beats. Let us suppose that we have a harmonic wave 
<f> = A cos 2iT[nt—kx), with amplitude A and frequency n. 
Suppose further that the amplitude A is made to vary 
with the time in such a way that at x = 0, A — a-\- 
b cos 2 -rrpt. If the wave is to move with velocity c = n/k, 
(f> must be some function of ct—x. So for general a-, 
A = a+b cos 2irp(t—kx' l n). This is known as amplitude 
modulation. The result is 

(f> = {a+6 cos 27Tp(t—kx'n)}cos 2n(nt—kx) 

— (I cos 2ir(nt—kx) 

+ o{ cos ^)] + cos 1») ^ . 


'i lie effect of modulating, or varying, the amplitude, is to 
introduce two new frequencies as well as the original one ; 
these new frequencies n±p are known as combination 
tones. I n the same way we can discuss phase modulation 
and frequencymodulation. (.Sec § 89, questions (11)-(12).) 

§ 83. If the velocities of § 81 are not the same (njk y 
not equal to n 2 jk 2 ), then the profile curves in fig. 25 move 
with a speed (wj-n,)/^-*,), which is different from that 
of the more rapidly oscillating part, whoso speed is 
(Hj-f «o)/(£,-)-A-.,). other words, the individual waves 
in fig. Jo advance through the profile, gradually increasing 
and then decreasing their amplitude, as they give place to 
other succeeding waves. This explains why, on the sea- 
shore, a wave which looks very large when it is some 
distance away from the shore, gradually reduces in height 

as „! t . m< ? ves , in ’ and ,na y even disappear before it is 
sufficiently close to break. 


This situation arises whenever the velocity of the 
waves i.e. their wave velocity V, is not constant, but 
depends on the frequency. This phenomenon is known 
< tspersion, \e deduce that in a dispersive system 
the only wave profile that can he transmitted without 
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change of shape is a single harmonic wave train ; any 
other wave profile, which may be analysed into two or 
more harmonic wave trains, will change as it is propagated. 
The actual velocity of the profile curves in fig. 25 is known 
as the group velocity U. Wo see from (6) that if the 
two components are not very different, V = njk, and 

U = (n l —n 2 )l(k l —k 2 ) = dn',dk. . . (7) 


In terras of the wavelength A, we have k = 1/A, so that 

• ' ,8) 

We could equally well write this 

u*W_ r+k ir_ r _ x dv (#| 

dk dk dk d\ 


Our calculation has considered just two waves. But 
the form of equation (7) shows that wo could equally well 
consider any number of waves superposed, provided that 
for any two of them n,— n 2 and k x - k 2 were sufficiently 
small for us to take their ratio constant and equal to dn dk. 
If this condition is not satisfied we have to go to a closer 
approximation, as in § 84. 

The importance of group velocity lies in the fact that 
the energy is propagated with this velocity. We have 
already met several cases in which the wave velocity 
depends on the frequency ; we shall calculate the group 
velocity for three of them. 

Surface waves on a liquid of depth h : 

The analysis of Chapter V, equation (32) shows that the 
velocity of surface wave3 on a liquid of depth h is given by 


V 2 = 




277 


According to (9) therefore, the group velocity is V — Xd V/dX, 





4 nh . 4 -h 

— — cosech — — 

A A 
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When h is small, the two velocities are almost the same, 
but when h is large, U = F/2, so that the group velocity 
for deep sea waves is one-half of the wave velocity. 
Equation (10) is the same as the expression obtained in 
§ 52, equation (47), for the rate of transmission of energy 
in these surface waves. Thus the energy is transmitted 
with the group velocity. 

Electric waves in a dielectric medium : 

The analysis in Chapter VII, § 69, shows that the wave 
velocity in a dielectric medium is given by 

V 2 = c 2 /cfi. 

We may put p. = 1 for waves in the visible region. Now 
the dielectric constant e is not independent of the frequency, 
and so V depends on A. The group velocity follows from 
(9) ; it is 



In most regions, especially when A is long, e decreases when 
A increases so that U is less than V. For certain wave- 
lengths, however, particularly those in the neighbourhood 
of a natural frequency of the atoms of the dielectric, there is 
anomalous dispersion, and U may exceed V. When A 
is large, we have the approximate formula 

e = A +fl/A t +C , /A 4 . 

It then appears from (11) that 

u = v ^ 

A+B, AH C/A 4 * 

Electric waves in a conducting medium : 

1 ho analysis in Chapter VII, § 76, shows that thoelectrio 
vector is propagated with an exponential term 

where = ^., h lpc. Thus V‘ = I =. J!L . According 
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to (7), the group velocity is 




If we suppose that a and /x remain constant for all 
frequencies, then this reduces to 

U = 2/y = 2V . . . . (12) 


The group velocity here is actually greater than the wave 
velocity. 

§84. Wo shall now extend this discussion of group 
velocity to deal with the case of more than two component 
waves. We shall suppose that the wave profile is split 
up into an infinite number of harmonic waves of the type 

e 2 m(tr-«n (13) 


in which the wave number k has all possible values ; we 
can suppose that the wave velocity depends on the 
frequency, so that n is a function of k. If the amplitude 
of the component wave (13) is a(k) per unit range of k, 
then the full disturbance is 



a(k) . e ini{kl ~ n,) Jk 


-30 



This collection of superposed waves is known as a wave 
packet. The most interesting wave packets are those in 
wliich the amplitude is largest for a certain value of k, 
say k 0 , and is vanishingly small if k—k 0 is large. Then 
the component waves mostly resemble e 2,n,lvr and 
there are not many waves which differ greatly from this. 
Wo shall discuss in detail the case in which 

a(k)=A . . . (15) 

This is known as a Gaussian wave packet, after the 
mathematician Gauss, who used the exponential function 
(15) in many of his investigations of other problems. 
A, a and k t are, of course, constants for any one packet. 
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Let us first determine the shape of the wave profile at 
t = 0. The integral in (14) is much simplified because 
the term in n disappears. In fact, 

00 

0) = J A . g2rriXi dk. 

— X 

On account of the term the only range of k 

which contributes significantly to this integral lies around 
*o I since when k— k 0 = l/\/a this term becomes er 1 , and 
for larger values of k—k 0 it becomes rapidly smaller, this 
range of k is of order of magnitude A k = 

In order to evaluate the integral, wo use the result * 



- X> 


This enables us to integrate at once, and we find that 

0) = A yj- e --' z, ’ a e-* 1 ** . . (18) 

The term er m ** represents a harmonic wave, whose 



wavelength A - l /.- 0 , and the other factors give a varying 
amplitude A yj - e If we take the real part of (18), 


• Gillosnio. Integration, p. 88. 
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<f)(x } 0) has the general shape shown in fig. -0. The outer 
curves in this figure are the two Gaussian curves 

«=±A yj”- 

and <f>(x, 0) oscillates between them. Our wave packet 
(14) represents, at t = 0, one largo pulse containing several 
oscillations. If we define a half-width as the value of x 
that reduces the amplitude to 1/c times its maximum 
value, then the half-width of this pulse is ( \/ <7 )/t7. 

At later times, />0, wo have to integrate (14) as it 
stands. To do this wo require a detailed knowledge of n 
as a function of k. If wo expand according to laylors 
theorem, wo can write 


where 


n = 7i 0 +a(&— k 0 )+f}(k-l: 0 ) 2 l2+... 


(19) 


a = (dnjdlc) 0 , /3 = (d'-nldlc*) 

As a rule the first tw r o terms are the most important, and 
if wo neglect succeeding terms, we may integrate, lining 
(17). The result is 

+ CD 

A e-W-to'e 2 "* 111 “ r(n# + a(i ' 


<f>{x, l) = J 


— CO 


= A ^/- e- 7T ' ,x - cu, ’ ia ■ e-’ !ia,z ~ n,,) - 


( 20 ) 


When t = 0, it is seen that this does reduce to (IS), thus 
providing a check upon our calculations. The last term 
in (20) shows that the individual waves move with a 
wave velocity njk 0 , but their boundary amplitude is given 

by the first part of the expression, viz. A yj- e n(z 

Now this expression is exactly the same as in (IS), di.iun 
in fig. 20, except that it is displaced a distance at to the 
right. Wo conclude, therefore, that the group as a uhole 
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moves with velocity a = (dnjdk) 0 , but that individual 
waves within the group have the wave velocity n 0 /k 0 . 
The velocity of the group as a whole is just what we have 
previously called the group velocity (7). 

If we take one more term in (19) and integrate to 
obtain <f>(x, t) we find that <f> has the same form as in (20) 
except that a is replaced by a+nfiit. The effect of this 
is twofold ; in the first place it introduces a variable 
phase into the term and in the second place it 

changes the exponential term in the boundary amplitude 
curve to the form 

—Tr 2 a[x—at) i 
6XP O'+TT-f-r- • 

This is still a Gaussian curve, but its half-width is increased 

t0 {(v*+n*Pt*)lon'}W . . . . ( 21 ) 

Wo notice therefore that the wave packet moves with 
the wave velocity 7iJk Qt and group velocity ( dn/dk ) 0 , 
spreading out as it goes in such a way that its half- width 
at time t is given by (21). 

The importance of the group velocity lies mainly in 
the fact that in most problems where dispersion occurs, 
the group velocity is the velocity with which the energy 

is propagated. We have already met this in previous 
paragraph. 

§ 85. W e shall next give a general discussion of the 

standard equation of wave motion y 2 <j> = — ^ , in which 

c 2 8t 2 

c is constant. We shall show that the value of cf> at any 
point P (which may, without loss of generality bo taken 
to bo the origin) may be obtained from a knowledge of 

(ho values of <f>, and on any given closed surface S, 

which may or may not surround P ; the values of <f> and 
its derivatives on S have to be associated with times 
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which differ somewhat from the time at which we wish 
to determine (f>p. 

Let us analyse tf> into components with different 
frequencies ; each component itself must satisfy the 
equation of wave motion, and by the principle of super- 
position, which holds when c is constant, we can add the 
various components together to obtain the full solution. 
Let us consider first that part of <f> which is of frequency 
p ; we may write it in the form 

tfr(x, y, z) e ikct , . . . (22) 

where k = 2 ttpIc (23) 

ip is the space part of the disturbance, and it satisfies the 
Poisson equation 

(V l +W = °- • • • (24) 

This last equation may be solved by using Green’s theorem.* 
This theorem states that if ifi 1 and ip 2 are any two functions, 
and S is any closed surface, which may consist of two or 
more parts, such that iJj 1 and ip 2 have no singularities inside 

it, then 

J{0 iVVi-^iVVa)* = 


S 


Fio. 27 

The volume integral on the left-hand side is taken over 
• See Rutherford, p. 65, equation (29). 



J{ 


Ti dn ri on 


dS. (25) 
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the whole volume bounded by S, and 8/dn denotes 
differentiation along the outward normal to dS. 

In this equation ip 1 and ift« are arbitrary, so we may 

e -Ur 

put ip 1 equal to tfj, the solution of (24), and tp 2 — , 

r being measured radially from the origin P. We take 
the volume through which we integrate to be the whole 
volume contained between the given closed surface S 
(fig. 27) and a small sphere 27 around the origin. We 
have to exclude the origin because becomes infinite at 
that point. Fig. 27 is drawn for the case of P within S ; 
the analysis holds just as well if P lies outside S. 

Now it can easily bo verified that y 2 i/r 2 = —k s tfi 2 , 

so that the left-hand side of (25) becomes J* 0 2 (V 2 +^‘ 2 )0 dr, 

and this vanishes, since (y 2 +k 2 )ifi = 0 by (24). The 
right-hand side of (25) consists of two parts, representing 
integrations over S and 27. On 27 the outward normal 
is directed towards P and hence this pait of the full 
expression is 

When wo make tho radius of Z tend to 0, only one term 
remains ; it is 



where dw is an element of solid angle round P. Taking 
the limit as r tends to zero, this gives us a contribution 
—477^/.. Equation (25) may therefore be written 
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Since by definition <f> = ifi(ryz)e itrt , we can write this 
last equation in the form 


+' = 


XdS 


( 20 ) 


where 


o iklct-r) gjj 

X = - — —Jie iklet ~ 

r bn 




= A —B +C, say. 

We may rewrito X in a simpler form ; for on account of 
the time variation of cp, ip e iklet ~ T) is tho same as <p, taken, 
not at time t, but at time l—r/c. If we write this symboli- 
cally [<f>]t- r i e , then B = ^ [<f>] t _ r/e . In a similar way, 

4 = , and » whero - for 

r Y dn \t-rie cr bn Ibl },- r/e 


~ 1 m 

~ r UnJi-* a 

, m 

example, — 

L dn \ t-rlc 


means that wo evaluate btybn as a 


function of x, y, z, t and then replace t by i—rjc. Wo 
call t—r/c the retarded time. We have therefore proved 
that 


* r = rj xis ‘ 


where 


X = l\ e J + I*|¥] .,27) 

f L^ n J(-f/e d n \ r ) V 1 T,C Cr&lL^Jf-r/e 


So far wo have been dealing with waves of one definite 
frequency. But there is nothing in (27) which depends 
upon tho frequency, and hence, by summation over all 
the components for each frequency present in our complete 
wave, we obtain a result exactly the same as (27) but 
without the restriction to a single frequency. 

This theorem, which is due to Kirchhoff, is of great 
theoretical importance ; for it implies (a) that the value 
of <f) may be regarded as the sum of contributions X/ \tt 
from each element of area of S ; this may bo allied the 
law of addition of small elements, and is familiar in a 
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law of addition of small elements, and is familiar in a 
slightly different form in optics as Huygens’ Principle ; 
and (6) that the contribution of dS depends on the value 
of <f>, not at time t, but at time l—r/c. Now rjc is the 
time that a signal would take to get from dS to the point 
P, 80 that the contribution made by dS depends not on 
the present value of (f> at dS, but on its value at that 
particular previous moment when it was necessary for a 
signal to leave dS in order that it should just have arrived 
at P. This is the justification for the title of retarded 
time, and for this reason also, [<f>]t- T / e is sometimes known 
as a retarded potential. 

It is not difficult to verify that we could have obtained 
a solution exactly similar to the above, but involving 
t-\-r/c instead of t— r/c ; wo should have taken in the 

g + iitf g— iifcf 

previous work to bo — — instead of . In this way 

we should have obtained advanced potentials, [<£] (+r / e , 
and advanced times, instead of retarded potentials and 
retarded times. More generally, too, we could have 
superposed the two types of solution, but we shall not 
discuss this matter further. 

In the case in which c = co, so that signals have an 
infinite velocity, the fundamental equation reduces to 
Laplace’s equation,* y-<f> = 0, and the question of time 
variation does not arise. Our equation (27) reduces to 
the standard solution for problems of electrostatics. 

§ 86. We shall apply this theory to the case of a source 
0 sending out spherical harmonic waves, and wo shall 
take S to be a closed surface surrounding the point P 
at which we want to calculate <f>, as shown in fig. 28. 
Consider a small element of dS at Q ; the outward normal 
makes angles 0, and d with QO and PQ, and these two 
distances are r x and r. The value of <f> at Q is given by 
the form appropriate to a spherical wave (see Chapter I, 
• See Rutherford, p. 67, equation (33). 
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equation (24)) : 


Thus 


a 

<Pq = — cos m(d— r x ) 

an x ar, 


(28) 


= acosflj j-^ cos m(c/— fj)— — sin m(c/— r t ) . 
I r i r i 

Now A = 277 /m, so that if fj is much greater than A, which 



Fio. 28 


will almost always happen in practical problems, we may 
put 

d<f> ma cos 0, . 

— = sin m(cl—r,). 

dn r, 


Also 


and 


2(1) = _I 

bn\r) r a 


cos a 


d<f> amc . 

— sin micl—Ty). 

Cl fj 


The retarded values are easily found, and in fact, from (27), 

„ ma cos 9. . , 

X = ' sinm(c/— [r-frj) 

rr i 

+ j-cob 6 cob m(ct— [r+r.]) a —cos 0 sin m(c/ — [r+rj). 

r ^i cr r l 



WAVES 



We may neglect the second term on the right if r x is much 
greater than A, and so 

X = — — (cos 0+cos 9j) s\nm(d— [r+rj) . (29) 

rr i 


Combining (29) with (2G) it follows that 

. 1 f inn n 

<f>p = — — I — (cos 9 + cos dj) sin m(cl— 

477 -j rr x 

a r l 

= — — (COS0 4- cosSj smm(ct— [rd-r^ciS . (30) 

-AJ rr t 


If, instead of a spherical wave, we had had a plane 
wave coming from the direction of 0, we should write 

<t>q = a cos m(cl — r^), 

r l now being measured from some plane perpendicular to 
OQ, and (30) would be changed to 

<t>r = ~ ^ J* (cos 9 + cos 0,) sin m(cl~ [r-f rJJiS. (31) 

Wo may interpret (30) and (31) as follows. The effect 
at P is the same as if each element dS sends out a wave 

.. r . , A /cos 0+cos 0,\ 
ot amplitude — I - h/S, A bemg the amplitude 

of the incident wave at dS ; further, these waves aro 
a quarter of a period in advance of the incident wave, 
as is shown by the term -sin rj) instead of 

cos vi(ct rj). I (cos 0+ cos 0 X ) is called the inclination 
factor and if, as often happens, only small values of 0 
and 0 X occur significantly, it has the value unity. This 
interpretation ot (30) and (31) is known as Fresnel’s 
principle. 

The presence of this inclination factor removes a 
difficulty which was inherent in Huygens’ principle ; this 
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principle is usually stated in the form that each element 
of a wave-front emits wavelets in all directions, and these 
combine to form the observed progressive wave-front. In 
such a statement there is notliing to show why the wave 
does not progress backwards as well as forwards, since 
the wavelets should combine equally in either direction. 
The explanation is, of course, that for points behind the 
wave-front cos 6 is negative with a value either exactly 
or approximately equal to —cos 6 V and so the inclination 
factor is small. Each wavelet is therefore propagated 
almost entirely in the forward direction. 

Now let us suppose that some screens are introduced, 
and that they cover part of the surface of S. II we assume 
that tho distribution of <f> at any point Q near the screens 
is the same as it would have been if the screens were not 
present, we have merely to integrate (30) or (31) over 
those parts of S wliich are not covered. This approxi- 
mation, which is known as St Venant’s principle, is not 
rigorously correct, for there will be distortions in the 
value of <f>Q extending over several wavclengtlis from the 
edges of each screen. It is, however, an excellent approxi- 
mation for most optical problems, where A is small ; 
indeed (30) and (31) form the basis of the whole theory 
of tho diffraction of light. With sound waves, on the other 
hand, in which A is often of tho same order of magnitude 
as tho size of tho screen, it is only roughly correct. 

§ 87. Let us illustrate this discussion with an example 
of tho analysis summarised in (31). Consider an infinite 
screen (fig. 29) which we may take to be the xij plane. 
A small part of this screen (large compared with the wave- 
length of the waves but small compared with other distances 
involved) is cut away, leaving a hole through which waves 
may pass. We suppose that a set of plane harmonic waves 
is travelling in the positive z direction, and falls on tho 
screen ; we want to find the resulting disturbance at a 
point P behind the screen. 
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In accordance with § 86 we take the surface 5 to be 
the infinite xy plane, completed by the infinite hemisphere 
on the positive side of the xy plane. We may divide the 



contributions to (31) into threo parts. The first part 
arises from the aperture, the second part arises from the 
rest of the screen, and the third part arises from the 
hemisphere. 

If the incident harmonio waves are represented by 
(f> = a cos m(ct—z) this first contribution amounts to 



(1-f- cos 0) sin m(ct—r)dS. 


We have put 0 t = 0 in this expression since the waves 
fall normally on to the xy plane. We shall only be con- 
cerned here with points P which lie behind, or nearly 
behind, the aperture, so that we may also put cos 0=1 
without loss of accuracy. This contribution is then 



sin m(ct—r)dS 



The second part, which comes from the remainder of 
the xy plane, vanishes, since no waves penetrate the 
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screen and thus there are no secondary waves starting 
there. 

The third part, from the infinite hemisphere, also 
vanishes, because the only waves that can reach this part 
of S are those that came from the aperture, and when 
these waves reach the hemisphere their inclination factor 
is zero. Thus (32) is in actual fact the only non-zero con- 
tribution and we may write 

(f>P =— sin m(ct—r)dS . . (33) 

Let P be the point (x, y, z) and consider the contribution 
to (33) that arises from a small element of the aperture at 
Q (£> fj' 0)* If OP = /, and QP = r, we have 

p = **+ y 2 + 2 2 , 

r*=(x-t) i +(y~ri) i +z * 
=P-2xf-2y V +P+T ) * . . (34) 

Let us make the assumption that the aperture is so small 
that plf 2 and 7j 2 // 2 may bo neglected. Then to this 
approximation (34) shows us that 

_ _ , xt+yr) 

~ — J e 



Again without loss of accuracy, to the approximation to 
which we are working, we may put 1/r = 1//, and then 
we obtain 

<f> P = — A sin {m(d —/)+«}, 

where 

A 2 = C 2 +5 2 , tan f = 5/(7, 

K 2 
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and 


C[x, y) = cos ^ (xf +yr j )d£dr ) , 

S(x, y) = ^J* sin (x£ + yr))d£dri 


( 35 ) 


Once we know the shape of the aperture it is an easy 
matter to evaluate these integrals. Thus, if we consider 
the case of a rectangular aperture bounded by the lines 
£ = ± a , 77 = ±[3, we soon verify that S = 0, and that 

+0 + /? 

c = cos yj {xZ+yvWv d t 

-a J -p 


4a sin pax sin pfiy 
A/ px py 



where p = 2-rr/Xf. If we are dealing with light waves, 
then the intensity is proportional to C 2 and the diffraction 
pattern thus observed in the plane 2 = / consists of a 
grill network, with zero intensity corresponding to the 
values of x and y satisfying eithor sin pax = 0, or 
sin pf3y = 0 but excluding x = 0 and y = 0. 


§ 88. The discussion of the last paragraph related to 
the case of plane waves falling normally on an aperture 
whose size, while large compared with the wavelength, 
was still small compared with the distance from the 
aperture to the screen on which the pattern was being 
observed. We might refer to this as diffraction at a 
pin-hole. But the equations (35) arise in another far 
more important way which we must now explain, and 
which is known as Fraunhofer Diffraction. 

Consider (fig. 30) a plane wave shown as A A' in the 
diagram falling normally on a convergent lens L. (L 
replaces the previous pin-hole.) This lens will convert 
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the plane wave into a spherical wave which converges 
at Z t the focus. On account of the finite size of the lens 
the focus is not perfect, and we ask the question : what 
will be the intensity observed at a point P in the focal 
plane through Z ? 



INCIDENT 

WAVEFRONT 


FlO. 30 


To answer this question it is convenient to draw the 
wavefront liOR' of a wave that has just left the lens. 
We may regard this as part of a spherical surface with 
centre Z and radius equal to the focal length /. If we 
take 0 as origin and OZ as axis of Z, then the coordinates 
(f, t), £) of any point Q on the surface satisfy the equation 


t.e. 


= 2/C 



Now by reasoning very similar to that used in § 87 we may 
argue that if P is near Z, the total effect at P is just the 
sum of separate contributions arising from all the elements 
d8 within the curved wavefront RR'. Let u8 suppose 


ISO 
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that the inclination factor (p. 144) may be put equal to 
unity, and that the amplitude at all points on RR' is a ; 
this is the same as the amplitude in the incident wave AA' . 
Thus 

<f>Q = a cos m{ct+f) . . • (38) 


Let us write (x, y, f) for the coordinates of the point P 
at which the observation is made, and put QP = r. Then 
the appropriate form of (30) is 


<f>p = — sin m(ct— r+f) dS . . (39) 


We may replace 1/r by l//in the first part of this integral : 
and note that 


r 2 = (*-£) 2 +(y-7?) a +(/-£) a 

= r-+x 2 +y*-2x£-2yr, by (37). 


In cases where this type of diffraction is important, / is 
large and x and y are small. We may therefore write 


so that effectively 


r a =P-2x£-2y 1 , 


r=f 





Combining (39) and (40) we have 



Thus, on integrating, 


<f> P = —A sin(mc/-f-€), . . . (41) 


where A and e are given by preciselj r the same formulae 
as in (35). 

This kind of analysis will apply particularly to the 
image of a star in a telescope of long focal length. The 
star is so far away that it may be regarded as giving out 
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a beam of parallel light. We have just shown therefore 
that the image of the star is not a point, but a pattern 
with maxima and minima, depending on the shape and 
size of the lens. For example, if a rectangular aperture 
(bounded by f = ±a, t) = ±0) is placed immediately 
behind the lens L, the diffraction pattern is a grill network, 
as in (36). And a circular aperture (question 9 at the end 
of the chapter) gives diffraction rings around Z. In any 
case the finite extent of the central maximum, or zone, 
will put limits to our power of resolving the light from 
two close stars. For if the geometrical images of the 
two stars lie within one another’s central zones, we shall 
experience difficulty in distinguishing whether there is 
really only one star, or two. But there is no space here 
to deal with this important matter any more closely. 

§ 89. We conclude this chapter with a discussion of 
the equation 

W = ■ ■ ■ 1421 

where p is some given function of x, y, z and t. When p = 0 
this is the standard equation of wave motion, whose solution 
was discussed in § 85. Equation (42) has already occurred 
in the propagation of electric waves when charges were 
present (Chapter VIII, equations (17') and (18)). 
may solve this equation in a manner very similar to that 
used in § 85. Thus, suppose that p(x, y, z, t) is expresse 
in the form of a Fourier series with respect to t, viz., 

p(x, y, z, t) = Za k {x, y, z)e‘ k " . • («) 

There may be a finite, or an infinite, number of different 
values of k, and instead of a summation over discrete 
values of k we could, if we desired, include also an integra- 
tion over a continuous range of values. W e shall discuss 
here the case of discrete values of k ; the student will 
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easily adapt our method of solution to deal with a 
continuum. 

Suppose that <f>(x, y, z, t) is itself analysed into com- 
ponents similar to (43), and let us write, similarly to (22), 

<£(*, y> z,t)=Z ifj k (x , y , z)e iUt . . (44) 

k 


the values of k being the same as in (43). If we substitute 
(43) and (44) into (42), and then equate coefficients of 
e ikct , we obtain an equation for ip k . It is 

(V a -r^ 2 )«Ajt = -4wz* . . (45) 


This equation maj' be solved just as in § 85. Using 
Green’s theorem as in (25), we put = ift k (x, y, z), 
e~ ikr 

V^ 2 = ~y> taking Z and S to bo the same as in fig. 27. 

With these values, it is easily seen that the left-hand side 
of (25) no longer vanishes, but has the value 

. . . (46) 

the integral being taken over the space between Z and S. 
The right-hand side may be treated exactly as in § 85, 
and gives two terms, one due to integration over Z, and 
the other to integration over S. The first of these is 

-4tt ifi k {x P , y P , 2 P ) . . . (47) 

The second may be calculated just as on p. 140. Gathering 
the various terms together, we obtain 


+ k 


<K(*e, 2/e, z P ) = J - t(T, r y ’ e~ ikT dt 
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Combining (43), (44) and (48) we can soon verify that our 
solution can be written in the form 

^,!/r,3 P ) = J [ ^'*+i|zrf.S. . (49) 


where X is defined by (27). This solution reduces to (27) 
in the case where p = 0, while it reduces to the well-known 
solution of electrostatics in the case where c = co. 

We have now obtained the required solution of (42). 
Often, however, there will be conditions imposed by the 
physical nature of our problem that allow us to simplify 
(49). Thus, if p(x, y, z, t) is finite in extent, and has only 
had non-zero values for a finite time t>t 0 , we can make 


X = 0 by taking S to be the sphere at infinity. This 
follows because X is measured at the retarded time t — r c, 
and if r is large enough, we shall have t — r'c<t 0 , so that 
[<£],_„„ and its derivatives will be identically zero on S. 
In such a case we have the simple result 

, , » [[p]t-rit 

<{>{xp, Vr< z r ) = I — — 

the integration being taken over the whole of space. 
Retarded potentials calculated in this way are \ci\ 
important in the Classical theory of electrons. 


dr, 


(50) 


§ 89. Examples 

(1) An observer who is at rest notices that the frequency 
of a car appears to drop from 272 to 25G per second ns 
the car passes him. Show that the speed of the car is 
approximately 23 m.p.h. How fast must ho travel in ilio 
direction of the car for the apparent frequency to rise to 
280 per second, and what would it drop to in that case T 

(2) Show that in the Doppler effect, when the source and 
obsorvor are not moving in the same direction, the formu x 
of § 80 are valid to give the various changos in frequency, 
provided that u and v donote, not the actual velocities, but 
the components of the two velocities along the direction in 
which the waves roach the observer. 
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(3) The amplitude A of a harmonic wave A cog 2; r(n$— fcr) 
is modulated so that A = a +6 cos 2irpt+c cos 3 2ttpL Show 
that combination tones of frequencies n±p 9 n±2p appear, 
and calculate their partial amplitudes. 

(4) The dielectric constant of a certain gas varies with the 
wavelength according to the law e = A+B/X 2 —C\ t 9 where 
A y B and C are constants. Show that the group velocity U 
of electromagnetic waves is given in terms of the wave velocity 
V by the formula 



A -2CA 3 


(5) In a region of anomalous dispersion (§ 83) the dielectric 

AX 1 

constant obeys the approximate law c = 1 +— — — . A more 

A — Aq 

A A 2 ( A 3 — Aq 3 ) 

accurate expression is c = 1 H — — — — , where A , B and 

(A Aq ) +BA 

A 0 are constants. Find the group velocity of electric waves 
in t heso two cases. 

(6) Calculate the group velocity for ripples on an infinitely 
deep lake. (§ 54, equation (54).) 

(7) Investigate the motion of a wavepacket (§ 84) for 
which the amplitude a is given in terms of the wave number k 
by the relation 

a(k) = 1 if |fc— & 0 | <k x 
= 0 otherwise. 


k 0 and being constants. Assume that only the first two 
terms o! the Taylor expansion of n in terms of k are required. 
Show t lat at time t the disturbance is 


<{>(Ty t) 


sin {2tt(s— oQfcJ 
n(x—at) 




where a = ( dn 'dk ) Q . Verify that the wavepacket moves as 
a whole with the velocity a. 

(8) Show that when dS is normal to the incident light 

(§ 86)> the inclination factor is j - ^~ CQS Plot this function 

2 

against 6, and thus show that each little element dS of a 



GENERAL CONSIDERATIONS 


1 55 


wave gives zero amplitude immediately behind the direction 
of wave motion. Using the fact that the energy is proportional 
to the square of the amplitude of <f>, show that, taken alone, 
each small element sends out 7/8 of its energy forwards in 
front of the wave, and only 1/8 backwards. 

(9) A plane wave falls normally on a small circular 
aperture of radius b. Discuss the pattern observed at a large 
distance / behind the aperture. Show that with the formula! 
of § 87, if the incident wave is <f> =a cos m(d—z), then 
5 = 0, and if P is the point (x, 0, /), then 



co 9 p( d f where p = 2itx/A/, 


= I " n ^ cos (pb cos 6 ) sin ; 0 dO. 

A/ Jo 

Expand cos ( pb cos 0) in a power series in cos 6, and hence 
show that 

where k = pb\2 = nbx/Xf. Since the system is symmetrical 
around the z axis, this gives the disturbance at any point 
in the plane z =/. It can be shown that the infinite senes 
is in fact a Bessel function of order unity. It gives rise 
diffraction rings of diminishing intensity for large values of *. 

(10) The total charge q on a conducting sphere of radius « 
is made to vary so that q = 4*a*a. where a •= 0 f " r '<°- an 
a = a 0 Sin pt for *>0. Show that if « = /* - . <1™ » 

the electric potential * at a distance R from the centre of the 

sphere is given by 


ct<R—a, 

R—a<ct<R+a, 


R+a<cl, <t> — 


pR r ' \ c 

iiracoo " n ' 11 

pR c 


. pa . ( R\ 

sin — sin py— c y 


(11) The wave represented by ^ - A cos 2n(nl kx + () 
suffers phase modulation in which « = a+6 cos 2n; . . 
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and p are constants, and 6* may be neglected. Show that in 
addition to the wave of given frequency n and amplitude A, 
combination tones appear, with frequency n ±p, and amplitude 
nAb . 

(12) The wave represented by <f> = A cos 2;r(nJ — Axr + c) 
suffers frequency modulation in which n = n 0 +o cos 2 npt. 
a, n 0 and p are constants, a* may be neglected and p<<n 0 . 
Show that in addition to a wave of frequency n 0 and amplitude 
A, there are four combination tones of frequency n 0 ±p±a 
and amplitude $4. It may be assumed that at is small. 


[answers : 1 .c/34, where c = velocity of sound, 249 per 

_ _.r. A Ao*A* 1 

second; 3. a +c/2, 6/2, c/4 ; 5. 1 A *_ Aq i )( A *_ a o *+A A») J ’ 

V f ! + ^AW-W’-Vm whera y . wave velocity ; 

V + «{( a* — a„*) s + u a*}* r 

2r,J 1 

f; = ic + v’J 



INDEX 


Tho Numbers refer to pages 


Absorption coefficient, 121 
Adiabatic, 87 
Advanced potentials, 142 
times, 142 
Ampere’s Rule, 104 
Amplitude, 2 
modulation, 132 
Amplitudes, partial, 33 
Anomalous dispersion, 134 
Antinodes, 6 

Bars and springs, longitudinal 
waves in, 61-59 
Basins, tides in, 66, 70, 74-76 
Boats, 131 
frequency, 131 
Boll, vibrations of, 49 
Bernoulli’s equation, 62, 72, 82, 
93, 94 

Boundary conditions, 1, 27, 30, 
38, 62, 69, 73, 82, 90, 94 
Brewster’s angle, 117 

Capillary waves, 81-84 
Chladni’s figures, 48 
Circularly polarised light, 1 1 1 
Combination tones, 132 
Compressible fluid, 87 
Condensation, 89 
Conductivity, 103 
Conical pipe, sound waves in, 96 
Constant of separation, 9 
Constitutive relations, 104 
Coordinates, normal, 36, 47, 63 

D’Alombort, 7 
Damping, 15, 39 
Decay, modulus of, 15 
Degenerative vibrations, 46 
Dielectric displacement, 102 
Diffraction of light, 145 
Diffraction at a pin hole, 145, 148 
theory, Fraunhofer, 148 
Dispersion, 132 
anomalous, 134 

Displacement current, Maxwell’s 
104, 119, 122 


Doppler effect, 129 
Drude, 119 

Electric and magnetic field 
strengths, 102 
waves, 102-127 

Elliptically polarised light, 1 1 1 
End correction, 90 
Energy, kinetic, 23, 33, 47, 54, 
78, 96 
loss of, 124 

potential, 24, 33, 47, 54, <8, 96 
rote transmitted, 79, 111, 153 
Equation of telegraphy, 15, 39, 
119 

wave motion, 1-20, 5 
wave motion, complex solu- 
tions, 16-17 
Exponential horn, 100 

Field strengths, electric and 

magnetic, 102 

Fraunhofer diffraction theory, 
148 

Free surface, 62, 82 
Frequency, 3 
modulation, 132 
Fresnel’s principle, 144 
Fundamental, 35, 49, 58, 91 


Gaussian wave packet, 135 
Gauss’ thoorem, 104 
General considerations, 1 1.8- 1 ol 


ound note, 35 

oup velocity, 81, 133, 135, 


io n lio 


Hagen and Rubens, 124 
Half-width, 137 
Harmonic wave, 2, 16-17 
Horn, 96, 100 

Huygens' Principle, 142, 144 

Inclination factor, 144 
Incompressible liquid, 60 
Index, rcfractivo, 106, 119 
Induction, magnotic, 102 


167 



158 


INDEX 


Intensity, 148 

Internal or total reflection, 118 
Isothermal, 87 

Joule heat, 124 

Kinetic energy in bars, 54 
liquids, 78 
membranes, 47 
sound, 96 
strings, 23, 33 
Kirchhoff, 141 
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polarised light. 111 
nodal, 6 
wave, 4 

Polarisation, plane of, 111 
Polarising anglo, 118 
Potential, advanced, 142 
electric, 106, 108 
energy in bars, 54 
energy in liquids, 78 
energy in membranes, 47 
energy in sound, 96 
energy in strings, 24, 33 
magnetic or vector, 106, 108 
retarded, 142, 153 
velocity, 60, 72, 87, 92 
Poynting vector, 104, 112, 124 
Pressure, radiation, 124 
Principle of superposition, 6, 
130, 135, 138 
Profile, wave, 2 

Progressive waves, 6, 13, 23, 26, 
28, 30, 40, 66, 71, 74, 77, 
83, 95, 100, 109-125 

Reduction to a steady wavo, 40, 
71 

Reflection coefficient, 27, 117, 
123 

of light waves, 113 
total or internal, 1 1 8 
Refraction of light waves, 113 
Refractive index, 106, 118 
complex, 120 
Relaxation, time of, 109 
Resistance, specific, 103 
Rotarded potential, 142, 153 
time, 141 
Ripples, 83 

Screen, 145, 146 
Separation constant, 9 
Skin effect, 125 
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Snell's law, 115 
Sound, velocity of, 89 
waves, 87-101 
Source, moving, 128 
Springs and bars, longitudinal 
waves in, 61-59 
vibration of, 55 

Stationary waves, 0, 32, 38, 45, 
48, 53, 75, 95 
Strings, normal modes, 31 
waves on, 21-42 
St. Venant's Principle, 145 
Superposition, principle of, 6, 
130, 135, 138 
Surface, free, 62, 82 
tension, 63, 81 
waves in liquids, 63, 72-81 

Telegraphy, equation of, 15, 39, 
119 

Tidal waves, 62, 63-72 
Time of relaxation, 109 
Tone, 35 

combination, 132 
Total or internal reflection, 118 
Transmission coefficient, 28, 117, 
123 

Transverse waves, 21, 44, 109 

Vector, Poynting, 104, 112, 124 
Volocity, group, 81, 133, 135, 
137, 138 
of light, 106 


Velocity, of sound, 89 
potential, 60, 72, 87, 93 
wave, 132 

Vibrations, degenerate, 46 

Wave, capillary, 81-84 
electric. 102-127 
harmonic, 2, 16-17 
in bars and springs, 51-59 
in liquids, G0-86 
in mombranes, 43-50 
long, in shallow water, 62 
longitudinal, 21, 61-59, 87 
motion, equation of, 1-20, 5 
number, 3 
on strings, 21-42 
packet, 135 
plane, 4 
profile, 2 

progressive, 6, 13, 23, 26, 28, 
30, 40, 66, 71, 74, 77, 83, 
95, 100, 109-125 
reduction to a steady, 40, 71 
sound, 87-101 

stationary, 6, 32, 38, 45, 48, 
53, 75, 95 
surface, 63, 72-81 
tidal, 62, 63-72 
trans verso, 21, 44, 109 
velocity, 132 
Wave-front, 4 
Wavolongth, 3 



PRINTED IN GREAT BRITAIN BY OLIVER AND BOYD LTD.. 


EDINBURGH 
































